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ABSTRACT 


As  the  inductively  coupled  plasma  (ICP)  becomes  in¬ 
creasingly  more  popular  as  a  spectrochemical  source,  the  need 
to  understand  analyte  emission  characteristics  at  a  fundamental 
level  becomes  more  urgent.  A  detailed  knowledge  of  the 
emission  characteristics  as  a  function  of  spatial  position  is 
an  important  part  of  this  characterization.  A  self-scanning 
linear  photodiode  array  measurement  system  has  been  used  to 
obtain  vertical  and  horizontal  spatial  profiles  of  emission 
intensity . 

Vertical  spatial  profiles  are  normally  bell-shaped,  with 
the  position  of  maximum  intensity  above  the  load  coil  varying, 
depending  on  ICP  instrumental  parameters  and  analyte  excitation 
characteristics.  This  dependence  has  been  investigated.  It 
has  been  found  that  the  spatial  structure  of  many  atom  lines 
can  be  predicted  theoretically.  This  involves  the  cal¬ 
culation  of  emission  intensity  for  a  particular  atom  line  as 
a  function  of  the  excitation  and  ionization  potential  of  that 
line,  and  source  "temperature"  gradients.  The  spatial  dep¬ 
endence  of  ion  emission  on  these  parameters  has  also  been 
investigated . 

The  photodiode  array  provides  an  almost  ideal  sensor 
for  the  acquisition  of  lateral  intensity  data  (horizontal 
profiles)  exactly  in  a  format  for  input  to  an  Abel  inversion 
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computer  program.  This  program  converts  spatially  integrated, 
lateral  profiles  to  spatially  resolved  radial  profiles.  A 
group  of  lateral  profiles  taken  at  different  heights  above 
the  load  coil  can  provide  a  complete  3  dimensional  radial  map 
of  analyte  emission  in  the  ICP.  These  detailed  maps  provide 
information  about  the  excitation  environment  experienced  by 
analyte  species. 

A  spatial  study  of  the  effect  of  an  excess  of  an  easily 
ionizable  element  (EIE) ,  such  as  Na,  has  been  undertaken, 
using  both  vertical  and  radial  emission  profiles.  This 
study  has  clarified  the  observations  of  other  ICP  spectro- 
scopists  studying  this  complex  system.  The  spatial  data  has 
provided  sufficient  information  to  postulate  several  possible 
mechanisms  for  the  effect. 
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CHAPTER  I 


INTRODUCTION 

In  the  evolution  of  any  analytical  technique  there 
are,  in  general,  two  complementary  activities  which  con¬ 
tribute  to  its  development.  These  are  an  investigation 
of  the  fundamental  properties  or  characteristics  of  the 
technique,  and  application  of  the  technique  to  routine  and 
specialized  problem  solving.  The  success  of  a  particular 
method  must  be  evaluated  in  terms  of  the  latter  activity, 
the  "ideal"  being  the  analysis  of  all  substances  at  all 
concentration  levels,  in  any  sample  matrix  without  inter¬ 
ferences.  The  value  of  the  former  activity  is  in  sug¬ 
gesting  ways  in  achieving  this  "ideal",  through  an  under¬ 
standing  of  the  mechanisms  of  operation  of  the  method.  It 
is  often  the  case,  that  the  fundamental  and  practical 
pursuits  do  not  progress  at  the  same  rate;  such  is  the 
case  with  Inductively  Coupled  Plasma-Optical  Emission 
Spectroscopy  (ICP-OES) .  As  a  method  for  multielement 
atomic  analyses,  it  has  been  used  for  a  wide  variety  of 
analysis  and  is  becoming  firmly  established  in  the  analytical 
chemists  arsenal  as  a  result  of  its  early  successes.  An 
understanding  of  its  fundamental  characteristics,  in 
particular  the  mechanism  of  energy  transfer  to  analyte 
species,  has  been  slower  in  coming,  in  large  part  due  to 
the  complex  nature  of  the  "plasma  state".  Of  course  the 
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practical  application  of  plasma  spectroscopy  can  and  will 
progress  quite  nicely  on  empirical  data  but  the  mechanistic 
understanding  is  important  for  the  technique  to  reach  its 
full  potential.  What  is  required  is  a  description  of  the 
interactions  between  analyte  species  injected  into  the 
plasma  gas  and  the  plasma  gas  itself,  as  well  as  inter¬ 
analyte  interactions  in  the  plasma  gas  environment.  Such 
a  description  should  be  capable  of  rationalizing  the 
experimental  data  available  about  analytical  plasmas. 

Spatial  studies  are  considered  an  important  aspect 
in  the  characterization  of  emission  sources  [1-3] .  The 
acquisition  of  spatial  information  involves  the  observation 
of  some  measurable  quantity,  such  as  emission  intensity 
from  an  analyte  line,  as  a  function  of  a  distance  co¬ 
ordinate  related  to  the  source.  Scheeline  and  Walters  [2] 
have  summarized  the  importance  of  such  spatial  studies; 

"Spatial  resolution  of  the  numerous  properties 
of  spectral  light  sources  may  give  insight  into 
the  mechanisms  of  operation  of  these  sources 
and  thence  improved  analytical  procedures." 

The  rationalization  for  this  statement  is  the  fact  that 

the  spatial  distribution  of  some  measurable  physical 

property  is  a  reflection  of  the  environment  which  caused 

the  appearance  of  that  property. 

The  general  theme  of  this  thesis  is,  as  suggested  by 

the  title,  a  study  of  the  spatial  distribution  of  analyte 

species  in  a  Radio  Frequency  Inductively  Coupled  Argon 
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Plasma  (RF-ICAP) .  The  investigations  carried  out  are 
fundamental  in  nature  and  the  data  presented  are  used  to 
rationalize  or  clarify  the  observations  of  other  re¬ 
searchers,  and  to  support  the  author's  own  speculations 
about  physical  processes  in  ICAP's.  Before  proceeding,  a 
brief  introduction  to  the  ICAP  is  warranted. 

What  is  a  Plasma? 

Martin  Umen,  in  his  book,  "An  Introduction  to  Plasma 
Physics"  [4]  defines  a  plasma  as  a  "conglomeration  of 
mobile  positively  and  negatively  charged  particles,  there 
being  approximately  the  same  density  of  positive  charges 
as  negative".  In  argon,  the  ionization  process: 

Ar°  t  Ar+  +  e"  (1) 

creates  mobile  argon  ions  and  electrons  which  are  the  pre¬ 
dominant  charged  particles  at  low  temperatures  (<15,000  K) . 

The  analytical  RF-ICAP  is  generally  a  weakly  ionized  plasma, 
the  degree  of  ionization  at  10,000  K  being  9%  at  1.1  atm 
[5]  at  equilibrium. 

There  are  a  variety  of  methods  used  for  producing  plasmas 
for  analytical  spectroscopy  [6].  These  include:  D.C.  plasma 
jets,  microwave  induced  plasmas,  capacitively  coupled  plasmas, 
and  RF-ICAP' s.  The  type  of  plasma  under  study  in  this  work 


is  the  RF-ICAP  or  ICP  for  short. 
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The  Inductively  Coupled  Plasma 


The  first  mention  in  the  literature  of  an  inductively 
coupled  plasma  was  by  Babat  [7]  in  1947.  The  area  lay 
fallow  until  in  1961  Reed  [8]  published  an  article  de¬ 
scribing  an  ICP  operating  at  atmospheric  pressure  sup¬ 
ported  by  flowing  argon  gas.  The  plasma  torch  consisted 
of  a  quartz  tube,  open  at  one  end  supplied  with  flowing 
argon  gas  at  the  other.  A  five  turn  copper  coil  was 
wrapped  around  the  quartz  tube  near  its  outlet  and  RF 
power  to  a  maximum  of  10  KW  at  4  MHz  was  delivered  to  the 
coil.  The  plasma  was  initiated  by  inserting  a  graphite  rod 
into  the  RF  field  which  when  heated  inductively  initiated  break¬ 
down  of  the  argon  gas  in  its  vicinity  after  which  the 
plasma  became  self-supporting.  Reed  also  introduced  the 
idea  of  vortex  stabilization  whereby  the  argon  gas  is 
introduced  into  the  torch  tangentially  causing  it  to  flow 
in  a  spiral  motion  up  along  the  walls  of  the  quartz 
tube.  This  centres  the  plasma,  recirculates  some  of  the 
plasma  gas  between  the  coils,  and  cools  the  quartz  tube 
to  prevent  melting.  This  proved  to  be  an  important  con¬ 
cept  and  has  been  widely  employed  for  analytical  plasmas. 

It  was  quickly  recognized  that  the  ICP  could  be  a 
useful  tool  for  spectroscopic  analysis  and  two  groups 
independently  published  their  application  of  the  Reed 
type  torch.  Greenfield  et  al.  [9]  described  a  torch 
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which  is  very  close  to  that  utilized  in  present  com¬ 
mercially  available  plasma  systems.  It  consisted  of  two 
concentric  glass  tubes  with  tangential  gas  inlets  on  each 
and  a  central  gas  tube  for  injecting  sample  into  the  plasma. 
Emission  from  the  elements  injected  was  observed  in  the 
tailflame  region  of  the  plasma.  This  paper  was  followed 
by  a  description  by  Wendt  and  Fassel  [10]  of  a  laminar 
flow  ICP  torch.  However  Fassel' s  group  eventually  con¬ 
verted  to  tangential  gas  inlets  favouring  Reed's  approach 
of  vortex  stabilization  [11,12].  From  these  beginnings 
the  use  of  plasmas,  the  ICP  in  particular,  as  a  source 
for  spectrochemical  analysis  has  blossomed  to  the  extent 
that  a  journal  has  been  initiated  which  is  devoted  entirely 
to  the  source  [13] . 

Some  popular  areas  of  application  for  the  ICP  have 
been  in  the  analysis  of  rocks  [14]  and  soils  [15,16],  in 
the  metals  industry  [17]  for  the  analysis  of  trace  el¬ 
ements  in  steels,  trace  elements  in  biological  materials 
[16,18]  and  natural  waters  [19],  and  in  forensic  science 
[20] .  The  ICP  has  also  proven  to  be  useful  for  the 
analysis  of  elements  in  oils  [21]  and  organic  compounds 
[22] .  For  good  historical  and  application  reviews  of 
plasmas  as  spectrochemical  sources  the  reader  is  referred 
to  Greenfield  et  al.  [23-25],  Barnes  [26],  and  Dahlquist 


and  Knoll  [16] . 
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Generating  and  Sustaining  the  Plasma 

A  schematic  diagram  outlining  the  main  features  of 
the  ICP  discharge  is  provided  in  Figure  1.  This  general 
approach  has  become  one  of  the  most  popular  and  is  that 
used  for  all  the  measurements  described  in  this  thesis. 

The  plasma  torch  itself  is  described  in  detail  in  Chapter 
II.  It  generally  consists  of  two  concentric  glass  tubes 
with  a  smaller  aerosol  tube  up  the  centre.  Argon  gas  is 
introduced  tangentially  at  the  base  of  the  torch  and 
spirals  up  through  the  load  coil.  Power  is  fed  into  the 
coil  from  an  RF  generator  operating  at  27.18  MHz  in  the 
range  of  1  to  3  Kw.  A  tesla  coil  is  used  to  "seed"  the 
incoming  gas  with  electrons  by  ionizing  a  portion  of  the 
argon.  The  electrons  produced  are  accelerated  in  the 
magnetic  field  produced  as  a  result  of  the  RF  power  ap¬ 
plied  to  the  coil  and  these  accelerated  electrons  collide 
with  neutral  argon  atoms  causing  further  ionization.  This 
can  be  represented  as: 

Ar  +  e  t-  Ar+  +  2e  (2) 

At  some  point  in  time  the  rate  of  production  of  argon 
atoms  by  the  forward  process  in  equation  2  will  be  bal¬ 
anced  by  the  reverse  process,  and  by  the  two  body  process: 

Ar+  +  e  +  Ar*  +  hv  (3) 
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Figure  1.  ”ain  features  of  the  ICP  discharge. 


such  that  an  equilibrium  is  established  in  which  the  rate 
of  production  of  argon  ions  is  balanced  by  the  rate  of 
argon  ion  consumption.  The  plasma  is  thus  sustained  by 
the  continuous  application  of  RF  power.  The  plasma  forms 
a  toroidal  shape  at  a  frequency  of  27  MHz  which  fac¬ 
ilitates  the  introduction  of  sample  into  the  base  of  the 
hot  gases  [11] . 

Normally  sample  is  introduced  into  the  plasma  in  the 
form  of  an  aerosol  gas,  although  powdered  samples  [27,28] 
or  solid  samples  [29-31]  may  be  introduced  using  appro¬ 
priate  techniques.  The  aerosol  gas  "punches  a  hole" 
through  the  base  of  the  plasma  forming  an  analyte  channel. 
The  processes  of  desolvation,  molecular  dissociation, 
ionization  equilibration,  and  excitation  take  place 
during  the  analyte's  progress  up  the  channel.  The  process 
of  excitation  and  subsequent  radiative  emission  cre¬ 
ates  a  "plume"  of  analyte  emission  within  the  analyte 
channel  approximately  30  mm  high  in  the  area  above  the  load 
coil.  It  is  the  delineation  of  the  spatial  structure  of  thi 
plume  as  a  function  of  analyte  species  and  sample 
matrix  that  is  one  of  the  major  goals  of  the  study  con¬ 
tained  in  this  thesis. 

Plasma  Physical  Properties  and  Spectral  Features 

The  physical  environment  of  the  RF  plasma  is  not 
characterized  by  combustion  processes  such  as  those  found 


in  flames  and  very  few  analogies  may  be  drawn  between  the 
two.  Rather  the  energy  contained  in,  and  dissipated  by, 
the  plasma  is  the  result  of  energy  transferred  from  the 
power  applied  to  the  load  coil  to  the  argon  species  present. 
The  energy  thus  absorbed  is  dissipated  in  three  main  ways 
[8] :  by  conduction  to  the  walls  of  the  quartz  tube,  by 
flow  of  hot  gases  away  from  the  load  coil,  and  by  the  em¬ 
ission  of  radiation. 

Radiative  emission  from  the  plasma  is  the  result  of 
two  main  processes  and  yields  the  characteristic  spectrum 
of  an  argon  plasma.  An  argon  line  spectrum  is  produced 
by  the  emission  of  light  as  a  result  of  transitions  from 
upper  excited  states  of  argon  neutral  atoms  (Ar*(l))  to 
lower  excited  states  (Ar*(2))  or  to  ground  state  argon 
atom  (Ar°)  through  the  processes: 

Ar*(l)  +  Ar*(2)  +  hv  (4) 

Ar*  (1)  ->  Ar°  +  hv  (5) 

where  v  is  the  transition  frequency.  In  addition  a  spec¬ 
tral  continuum  or  background  is  produced  as  a  result  of 
the  radiative  recombination  process  given  in  equation  (3). 
This  recombination  of  a  free  electron  produces  a  continuum 
which  is  characteristic  of  argon  [32],  Emission  from  argon 
ions  has  not  been  observed  in  the  ICP .  The  entrainment  of 
air  into  the  tailflame  produces  emission  from  NO  molecular 
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species  in  the  spectral  region  200-300  nm  [33] .  Mol¬ 
ecular  emission  from  OH  is  observed  in  the  region  around 
306  nm,  and  hydrogen  atomic  lines  are  observed  when  aqueous 
aerosols  are  introduced  into  the  plasma.  Some  spectra  of 
the  argon  ICP  have  been  published  demonstrating  these 
features  [33,34] . 

The  continuum  and  argon  line  intensity  is  greatest 
in  the  spatial  region  within  the  load  coils  and  drops 
quickly  with  increasing  distance  away  from  the  load  coils. 
The  iso-intensity  plot  of  the  emission  from  Ar  I  at  415.8 
nm  provided  in  Figure  2  demonstrates  this.  This  figure 
is  a  plot  of  the  radial  (from  the  centre  of  the  plasma  to 
the  edge)  intensity  distribution  of  Ar  I  emission;  the 
numbers  on  the  lines  representing  the  relative  intensities. 
It  can  be  seen  that  the  intensity  is  significantly  reduced 
in  the  centre  of  the  plasma  as  a  result  of  the  aerosol 
channel.  It  can  also  be  noticed  that  the  argon  line  and 
continuum  intensity  will  be  relatively  low  in  the  region 
15-20  mm  above  the  load  coil  where  analytical  observations 
are  usually  made.  This  is  good  from  the  standpoint  of 
maximum  signal  to  background  ratios. 

One  of  the  most  striking  features  of  the  ICP  and  one 
which  makes  it  particularly  suitable  for  emission  spectro¬ 
scopy  is  its  relatively  high  temperature.  Values  for  the 
measured  temperature  vary  depending  on  the  measurement 
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Figure  2.  Iso-intensity  profile  of  emission  from  Arl  415.8  nm  (heavy  lines  mark 
torch  edge) . 


method  employed  and  the  type  of  temperature  measured,  but 
most  workers  report  values  in  the  range  of  4000  K 
to  8000  K  [7,35-38].  The  practical  outcomes  of  these  high 
temperatures  are  several.  Firstly,  the  high  tem¬ 
peratures  mean  that  there  is  sufficient  energy  available 
to  significantly  populate  upper  levels  of  atoms  leading 
to  an  analyte  emission  spectrum  rich  in  analytical  lines, 
some  with  high  excitation  energies.  Even  elements  which 
have  been  difficult  to  analyze  by  emission  spectroscopy 
due  to  high  excitation  potentials,  such  as  sulphur,  iodine, 
phosphorus,  arsenic  and  selenium  have  been  successfully 
determined  using  the  ICP  [39-40].  Of  course  one  of  the  com¬ 
plications  associated  with  an  emission  spectrum  which  is 
rich  in  analysis  lines  is  the  possibility  of  spectral 
overlaps  [16]  which  may  turn  out  to  be  one  of  the  most 
significant  problems  associated  with  simultaneous  multi¬ 
element  analysis  with  ICP's.  Even  the  high  resolution 
afforded  by  echelle  spectrometers  does  not  completely 
overcome  these  problems  [41] .  Many  of  the  most  sen¬ 
sitive  analysis  lines  are  ionic  emission  lines,  partially 
due  to  the  fact  that  at  these  temperatures  many  species 
are  highly  ionized,  but  also  resulting  from  an  "ion  line 
advantage".  The  possible  source  of  this  advantage  will 
be  discussed  in  a  subsequent  section  but  basically  it 
means  the  population  of  many  emitting  ion  energy  levels 
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are  one  to  three  orders  of  magnitude  higher  than  that 
predicted  for  the  observed  excitation  temperature  [42]. 

The  analyte  channel  is  a  narrow  emitting  core  sur¬ 
rounded  by  hot  optically  thin  gases.  This  structure  leads 
to  very  low  levels  of  self  absorption  [43] ,  which  in  turn 
provides  analytical  curves  with  a  very  large  linear  work¬ 
ing  range;  about  4  to  5  orders  of  magnitude  [44].  From 
an  analytical  standpoint  this  means  that  samples  of  widely 
varying  composition  can  be  analyzed  without  the  need  for 
dilution . 

A  further  benefit  of  the  high  temperatures  is  that 
many  of  the  classical  vaporization  interferences  associated 
with  emission  and  absorption  flame  spectroscopy  are  re¬ 
duced  or  eliminated.  For  example,  the  suppression  of 
emission  due  to  the  presence  of  phosphate  or  aluminium 
concomitants  has  been  shown  to  be  negligible  in  the  ICP 
[45-48].  Presumably  the  plasma  is  sufficiently  energetic 
to  dissociate  any  compounds  which  may  form  before  or 
after  the  nebulization  stage.  It  has  been  reported  that 
the  effect  of  the  presence  of  easily  ionizable  elements 
(EIE's)  such  as  Na  may  be  reduced  to  a  negligible  level 
under  suitable  operating  conditions  [47].  There  is, 
however,  much  confusion  in  the  literature  about  this  effect 
and  the  exact  nature  of  the  influence  of  EIE's  on  analyte 
emission.  The  classical  interpretation  in  flames  of  a 
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shift  in  the  ionization  equilibrium  between  analyte  ion 
and  neutral  atom  species  does  not  seem  to  apply  in  the 
ICP.  This  is  rationalized  on  the  basis  that  the  already 
high  density  of  electrons  in  the  plasma  is  not  changed  by 
addition  of  EIE's.  Nevertheless,  there  is  an  effect  on 
net  analyte  emission  intensity  due  to  these  elements 
(Na,  Li,  K,  Rb,  Cs)  and  this  area  has  been  singled  out  as 
one  in  which  a  detailed  spatial  study  is  warranted.  This 
study,  presented  in  a  later  chapter,  has  led  to  a  clar¬ 
ification  of  the  "effect",  showing  that  the  variability  in 
the  influence  of  EIE's  reported  by  different  workers  is 
probably  spatial  in  origin. 

Plasma  Temperature 

Temperature  has  been  mentioned  several  times  through¬ 
out  the  preceding  section  on  plasma  properties.  The 
interpretation  of  what  is  meant  by  temperature  and  its 
measurement  is  the  subject  of  a  massive  body  of  plasma 
literature  and  will  be  considered  only  briefly  in  this 
thesis.  A  discussion  of  temperature  is  integrally  bound 
with  a  discussion  of  equilibrium  in  the  ICP. 

Thermodynamic  equilibrium  is  a  state  in  which  de¬ 
tailed  energy  balance  applies,  that  is,  "every  process  of 
energy  transfer  is  perfectly  balanced  by  the  reverse 
process"  [21] .  This  is  not  achievable  in  an  emission 
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source  since  the  loss  of  energy  from  the  system  by  em¬ 


ission  of  radiation  is  not  compensated  for.  For  this  reason 
the  concept  of  local  thermodynamic  equilibrium  (LTE)  is 
introduced,  whereby  all  energy  transfer  processes  except 
emission  of  radiation  are  in  thermodynamic  equilibrium. 

Under  the  condition  that  LTE  applies, a  single  system 
temperature  defines  the  following  distributions. 

(1)  The  velocity  of  particles  should  be  given  by  Maxwell's 
equation , 

f  (v)  =  (»-  Pm-)  eXp  (-mv^/2kT  )  4ttv^  dv 
27rkT  c  g 


(2) 


where , 
f  (v) 

m 

k 


the  probability  function  that  a  particle 
has  velocity  in  the  region  v  and  v  +  dv. 

mass  of  the  particle 

Boltzmann  constant 


Tg  =  gas  kinetic  temperature. 

For  the  velocity  distribution  of  electrons,  this 
temperature  is  called  T  ,  the  electron  temperature. 
The  population  of  excited  states  should  be  given  by 
the  Boltzmann  equation. 


n  g 

—  =  (— )  exp  [  —  ( E  -E  ) /kT  ] 
n  g  ^  g  p  exc 

P  P 


where , 

n^  and  n^  =  population  of  excited  states  g  and  p 
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g^  and  =  degeneracies  of  states  g  and  p 

E  and  E  =  excitation  energies  of  states  g  and  p 
9  P 

T  =  excitation  temperature 

0  X  c 

(3)  The  distribution  of  molecules  and  their  dissociation 
products  should  follow  the  mass  action  law  of  Guldberg 
and  Waage. 

(4)  The  distribution  of  atoms  and  ionization  products 

should  be  described  by  the  Saha  equation;  for 

M?  t  M+  +  e” 

D  3 

the  equilibrium  constant  is: 
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where  is  the  Saha  equilibrium  constant  expressed 

in  terms  of  number  density  (n)  of  species.  The 
equilibrium  constant  may  also  be  expressed  in  terms 
of  partial  pressures, 
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where  p .  .  , 
iD 


the  partial  pressures  of  ion, 
electron  and  atom  species  of 
element  j ; 


m  =  electron  mass; 


h  =  Planck's  constant; 

Z. Z  .  =  ion  and  atom  partition  function; 
ID  aD 
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=  ionization  energy  of  the  atom  ; 

T\  =  ionization  temperature . 

These  definitions  lead  to  the  possibility  of  several 
different  kinds  of  temperatures  which  are  commonly  specified 
for  a  particular  source.  These  are  the  electron  temp¬ 
erature  (T  )  defining  free  electron  energies  (velocities), 
gas  dynamic  temperature  (T  )  defining  the  gas  particle 
energies  (velocities) ,  excitation  temperature  (Texc)  de- 
fining  atom  excited  state  populations,  rotational  temp¬ 
erature  (T^.  )  defining  populations  of  excited  rotational 

states  of  molecules,  and  ion  temperature  (T^)  defining 
the  distribution  of  atom  and  ion  populations.  Each  type 
of  temperature  may  be  measured  by  spectroscopic  means  by 
observing  the  particular  system  with  which  it  is  as¬ 
sociated.  For  a  system  which  is  in  LTE  the  condition 

T  =  T  =  t  =  t  =  t 
g  e  exc  rot  i 

will  apply.  The  ICP  is  generally  considered  not  to  be 

in  LTE  [35-38,42,49,50].  The  results  of  Kornblum  and 

deGalen  [35]  obtained  for  an  ICP  operating  at  50  MHz 

indicated  that  T  <  T  <  T .  at  both  low  and  high 

g  exc  i 

aerosol  flow  rates. 

The  temperatures  which  have  been  reported  in  the 
literature  vary  widely  but  a  summary  of  representative 
values  is  provided  in  Table  X.  It  should  be  borne  in  mind 
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TABLE  I 


SURVEY  OF  PLASMA  TEMPERATURES 


Worker 

T  (K) 

exc 

T±  (K) 

T  ^  (K) 

rot 

Alder  et  al.  [38] 

5230-8510 

6320-8440 

- 

Kornblum,  deGalan  [35] 

2800-8600 

8000-9000 

1700-3800 

Kalnicky,  et  al.  [36] 

4100-6100 

- 

- 

Jarosz,  Mermet  [50] 

5500-6600 

9200 

- 

Alder,  Mermet  [51] 

5000-7000 

- 

6100 

Jarosz,  et  al.  [37] 

4800 

6700-7400 

- 

that  the  ICP  is  not  spatially  homogeneous  so  that  dif¬ 
ferent  regions  will  exhibit  different  values  for  temp¬ 
erature.  Also  there  is  significant  variation  in  the 
operating  parameters  used  in  the  different  studies  in¬ 
cluding  RF  power  level  and  frequency,  aerosol  flow  rate, 
and  spatial  region  observed.  For  these  reasons  it  is 
probably  not  valid  to  intercompare  the  various  values 
given  in  Table  I,  and  the  publications  should  be  con¬ 
sulted  for  specific  measurement  details.  It  will,  however, 
give  the  reader  some  feel  for  the  relative  magnitude  of 
the  various  temperatures.  Specific  methods  for  the  meas¬ 
urement  and  interpretation  of  source  temperatures  have 
not  been  included  in  this  discussion  although  results  of 
the  measurement  of  excitation  temperatures  will  be  included 
in  a  latter  chapter.  Reif  et  al.  [52-55]  and  Boumans  [56] 
discuss  these  matters  in  some  detail. 

Excitation  Mechanisms  in  the  ICP 

Elementary  particles  interact  by  virtue  of  grav¬ 
itational,  nuclear,  and  electromagnetic  forces.  The  sum 
total  of  these  interactions  is  termed  a  collision.  These 
collisions  may  be  either  elastic  in  which  the  total  kin¬ 
etic  energy  of  the  colliding  particles  is  conserved,  or 
inelastic,  in  which  the  total  energy  is  conserved  but 
kinetic  energy  is  not.  In  the  latter  case  kinetic  energy 
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may  be  converted  to  excitation  energy,  eventually  being 
lost  radiatively.  This  type  of  collision  is  termed  a 
"collision  of  the  first  kind".  An  example  in  which 
kinetic  energy  of  colliding  partners  A  and  B  is  converted 
to  excitation  energy  in  B  is  provided  by: 

A  +  B  ->  A  +  B*  (6) 

The  reverse  of  this  process  in  which  stored  excitation 
energy  is  partially  converted  to  kinetic  energy  of  the 
colliding  partners  is  termed  a  "collision  of  the  second 
kind".  Both  of  these  types  of  inelastic  collisions  are 
likely  to  be  significant  in  the  ICP. 

The  very  existence  of  the  plasma  state  in  the  argon 
ICP  implies  the  presence  of  three  elementary  particles. 
These  are  argon  neutral  atoms  (Ar) ,  argon  ions  (Ar+) 
and  electrons  (e  ) .  Initially  the  mode  of  production 
of  argon  ions  is  probably  due  to  impact  with  "fast" 
electrons  which  have  been  accelerated  in  the  magnetic 
field  produced  in  the  region  of  the  load  coil.  The 
ionization  energy  of  Ar  is  15.75  ev  so  the  electron  must 
possess  at  least  this  energy  to  ionize  the  argon  atom. 

Any  energy  in  excess  of  15.75  ev  is  carried  away  by  "slow" 
electrons  involved  in  the  collision.  This  process  is 
illustrated  in  Equation  7. 

Ar  +  e,.  Z  Ar+  +  e  +  e  (7) 

f  s  s 
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where  ef  is  a  fast  electron  and  e  is  a  slow  electron, 

fast  and  slow  being  terms  which  describe  the  relative 

amount  of  kinetic  energy  possessed  by  the  electrons.  The 

collision  of  fast  electrons  with  energy  greater  than  about 

* 

11.5  ev  produces  excited  states  of  the  argon  atom  (Ar  ). 
This  process  is  illustrated  in  Equation  8.  The  radiative 

_  *  _ 

Ar  +  e,  ^  Ar  +  e  (8) 

f  s 


decay  of  argon  excited  in  this  manner,  and  argon  excited 

by  ion-electron  recombination  (Equation  3)  gives  rise  to 

the  Ar  line  spectrum. 

3  3 

The  P 2  and  states  of  the  argon  atom  are 

"metastable"  levels  [57]  and  as  such  have  relatively  long 
lifetimes  compared  to  radiating  states.  The  relatively 
long  lifetimes  of  these  atoms  means  that  they  can  play  a 
role  in  the  collision  processes  in  the  ICP .  Thus  pro¬ 
duction  of  argon  ions  in  the  less  energetic  regions  of  the 

plasma  could  be  the  result  of  the  collision  of  argon 

* 

metastable  species  (Ar^)  with  relatively  slow  electrons 
[42].  The  energies  of  the  metastable  levels  are  11.55 
and  11.71  ev.  Thus  the  electrons  need  only  possess  energy 
in  the  range  of  4  ev  to  4.2  ev  to  initiate  this  process. 
This  process  is  illustrated  in  Equation  9. 


e 
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Ar+  + 


-> 


(  9) 
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In  this  way  the  plasma  state  can  exist  above  the  energy 
addition  region.  Decay  of  the  plasma  state  is  the 
result  of  the  radiative  recombination  and  radiative  de¬ 
excitation  processes  shown  in  Equations  10,  11,  and  12. 

Ar+  +  e  Z  Ar*  +  hv  (cont)  (10) 

s 

Ar+  +  e  +  e  Z  Ar*  +  e^  (11) 

s  s  f  ' 

Ar*  +  Ar °  +  hv  (line)  (12) 

When  analyte  species  are  introduced  into  the  plasma 
gas  these  argon  species  and  electrons  give  up  at  least 
part  of  their  energy  to  these  analyte  atoms.  Analyte  atoms 
(X)  may  be  excited  by  electron  impact  [58] . 

X  +  e,  Z  X*  +  e  (13) 

f  ^  s 

where  the  energy  of  the  fast  electron  must  be  greater  than 
the  excitation  potential  of  X*.  The  probability  of  the 
forward  process  is  greatest  when  the  energy  of  the  fast 
electron  just  exceeds  the  excitation  potential  of  X*. 
Analyte  ions  may  also  be  produced  by  electron-atom 
collisions  [58] , 

X  +  e~  Z  X+  +  e”  +  e~  (14) 

f  s  s 

where  the  energy  of  the  fast  electron  must  exceed  the 
ionization  potential  of  the  analyte  atom.  More  energetic 
electrons  can  produce  excited  state  ions 
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X  +  t  X+*  +  e  +  e  (15) 

f  s  s 

Analyte  ions  can  recombine  with  slow  electrons  into 
excited  states  [59],  in  two  body: 

X+  +  e  X  +  hv  (cont)  (16) 

or  three  body  recombination, 

X+  +  e”  +  e~  t  X*  +  e~  (17) 

s  s  f 

resulting  in  the  production  of  excited  state  analyte 
atom  species. 

Alternate  mechanisms  for  the  generation  of  emitting 
species  may  involve  argon  species.  Metastable  argon  has 
been  implicated  in  ionization  and  excitation  process 
through  Penning  ionization  reactions  [59-61] .  In  a 
scheme  in  which  this  type  of  reaction  is  present,  analyte 
atoms  are  ionized,  and/or  ionized  and  excited  depending  on 
the  energies  involved  with  excess  energy  being  carried 
away  in  the  form  of  kinetic  energy  of  the  free  electron. 
Thus  the  reaction  [42] 

*  +  - 

X+Ar+X+Ar+e  (18) 

m 

will  produce  analyte  ions  for  any  species  whose  ionization 
potential  is  less  than  the  metastable  values  of  11.55  and 
11.71  ev.  Similarly  the  reaction 

*  +  * 

X  +  Ar  +  X 
m 


+  Ar  +  e 


(19) 
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will  produce  excited  analyte  ions.  Brockaert  et  al . 

[62]  have  provided  indirect  evidence  that  this  reaction 
may  be  important  in  the  ICP.  They  have  shown  that  em¬ 
ission  from  yttrium  ion  levels  whose  sum  of  ionization  and 
excitation  energies  is  in  the  range  of  10-12  eV  are  most 
affected  by  the  addition  of  an  argon  metastable  quencher 
(H^BO^,  or  Na^B^O^  as  a  concomitant). 

Direct  excitation  of  atoms  or  ions  is  also  possible 
through  collisions  of  the  second  kind  with  argon  metastable 
species.  Thus  atoms  can  be  excited  by, 

*  * 

Ar+X+Ar+X  (20) 

m 

and  analogously  for  ions 

*  +  *4- 

Ar+X+Ar+X  (21) 

m 

The  energy  restrictions  on  these  two  reactions  are 
quite  critical,  in  that  the  excitation  energy  and  argon 
metastable  energy  should  match  closely  for  highest  pro¬ 
bability  of  energy  transfer.  Beenekker  [63]  observed  for 
argon  and  helium  microwave  induced  plasmas, 


"The  ions  of  analyte  atoms  produced  by  Penning 
ionization  are  excited  by  metastable  species 
of  the  plasma  support  gas,  giving  rise  to  the 
ion  line  emission  observed.  When  this  process 
is  not  possible  for  reasons  of  energy,  these 
ions  recombine  with  an  electron  to  yield  ex¬ 
cited  neutral  atoms  giving  rise  to  the  observed 
atom  line  emission." 
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Beenakker  also  suggested  the  possibility  of  charge  trans¬ 
fer  reactions  of  the  type; 

Ar+  +  X  +  Ar  +  X+*  (22) 

which  also  has  the  energy  constraint  that  the  analyte 
ionization  plus  excitation  energy  should  closely  match 
the  argon  ionization  energy.  This  type  of  charge  transfer 
reaction  is  considered  to  be  important  in  the  energy 
region  15-16  ev  [64],  the  ionization  energy  of  the  argon 
atom  being  15.75  ev. 

Jarosz  and  Mermet  [50]  also  note  the  presence  of  some 

lines  of  sulphur  ions  for  which  the  excitation  energy  is 

in  the  region  of  30  ev  and  suggests  that  argon  ion  meta- 
+  * 

stables  (Arm  )  play  a  role  through  the  mechanism; 

+  *  +* 

Ar+X+Ar+X  (23) 

m 

+  *  . 

the  Ar^  species  having  energies  m  the  range  32.13  to 
33.45  ev. 

It  is  unlikely  that  any  one  of  these  mechanisms  is 
the  sole  means  of  excitation  in  the  ICP.  Instead,  depending 
on  the  energy  requirements  of  the  analyte  species,  the 
mechanism  which  leads  to  the  most  favourable  probability 
of  energy  transfer  must  be  considered.  Additionally  sev¬ 
eral  of  the  mechanisms  may  be  acting  simultaneously  to 
produce  excited  states,  further  complicating  the  inter- 
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pretation  of  plasma  phenomena.  A  study  of  the  spatial 
structure  of  emission  in  the  plasma  will  certainly  aid 
in  delineating  which  of  the  mechanisms  should  be  con¬ 
sidered  most  important  in  the  analytical  ICP. 

Spatially  Resolved  Measurements 

The  well  known  adage,  "a  picture  is  worth  a  thousand 
words",  is  particularly  applicable  to  scientific  in¬ 
vestigation.  However,  often  it  is  not  sufficient  to 
simply  present  a  photograph  or  an  artistic  sketch  of  an 
emission  source,  because  many  of  the  physical  parameters 
which  characterize  the  source  cannot  be  directly  observed, 
or  often  a  parameter  of  interest  is  obscured  by 
another.  In  these  cases  the  investigator  must  use 
scientific  instrumentation  to  isolate  the  parameter  of 
interest  and  devise  some  system  whereby  the  spatial 
variation  of  this  parameter  may  be  mapped.  In  this  way 
the  physical  dependence  of  one  parameter  upon  another  can 
be  accessed  by  a  comparison  of  the  spatial  similarities. 

An  example  of  this  is  found  in  a  study  of  free-atom 
formation  processes  in  flames  undertaken  by  Cowley  et  al. 
[65].  They  were  able  to  correlate  the  spatial  character¬ 
istics  of  metal  oxide  and  free  atom  species  with  the 
spatial  characteristics  of  flame  produced  molecular  species 
such  as  CN,  CH,  C9  and  OH  to  account  for  the  favourable 


. 


characteristics  of  the  oxy-acetylene  flame. 

In  a  similar  manner,  Walters  [66]  was  able  to  cor¬ 
relate  the  spatial  dependence  of  emission  intensity  from 
metals  in  a  spark  discharge  with  the  energy  characteristic 
of  the  line  being  observed  [66]  ,  the  less  ionized,  less 
excited  species  being  observed  further  from  the  spark 
channel.  On  the  other  hand  continuum  intensities  peaked 
within  the  channel.  Armed  with  this  information  he  was 
able  to  postulate  possible  excitation  mechanisms  in  the 
spark  source.  Also  he  suggested  making  use  of  this 
spatial  separation  to  improve  signal  to  noise,  and  signal 
to  background  ratios  by  observing  emission  from  off  the 
channel  axis. 

Most  workers  in  the  field  of  plasma  spectroscopy 
recognize  the  utility  of  spatial  profiling  and  have  made 
use  of  the  technique  to  augment  their  studies.  The 
spatial  distribution  of  temperature  and  electron  number 
densities  in  the  ICP  have  been  studied  by  a  number  of 
workers  [3,35-37,67],  leading  to  the  suggestion  of  a 
non-LTE  plasma.  Also  the  spatial  effect  of  plasma  op¬ 
erating  conditions  on  emission  intensities  has  been 
investigated  [68-71].  These  studies  show  a  spatial 
structure  which  is  highly  dependent  upon  RF  power  and 
aerosol  flow  rate.  Also  the  spatial  effects  of  inter¬ 
ferences  have  been  studied  in  microwave  plasmas  [72-74], 
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DC  plasmas  [75]  and  ICPs  [48,71,76,77].  The  results  of 
these  studies  will  be  discussed  in  greater  detail  in  a 
subsequent  chapter. 

The  usual  modes  of  acquisition  of  spatial  information 
from  the  ICP  are  vertical,  that  is  parallel  to  the  plasma 
axis,  or  horizontal,  that  is  at  right  angles  to  the  plasma 
axis.  The  spatial  information  thus  gathered  is  the 
integral  of  emission  from  all  depths  from  within  the 
source.  From  a  fundamental  standpoint  it  is  desirable  to 
deconvolute  these  lateral  emission  intensities  such  that 
radial  distributions  can  be  obtained.  The  radial  in¬ 
formation  provides  emission  intensity  information  as  a 
function  of  distance  from  the  central  axis  of  the  dis¬ 
charge.  The  transformation  from  lateral  to  radial  in¬ 
formation  is  usually  performed  using  a  mathematical  techni¬ 
que  called  an  "Abel  inversion"  [78]  .  There  are  a  variety 
of  methods  for  implementing  this  for  horizontal  profiles 
and  it  has  been  used  extensively  for  the  study  of  spectro¬ 
scopic  sources.  Much  of  the  early  work  in  this  area  was 
in  the  domain  of  the  plasma  physicist  [79,80],  but  more 
recently  the  method  has  been  applied  to  spark  discharges 
[1,81,82]  and  ICPs  [3,36,37,67,76,83]. 

Traditionally  the  acquisition  of  detailed  spatial 
profiles  has  involved  the  translation  of  monochromator 
and  detection  system,  a  data  point  being  acquired 


' 
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for  each  translation  increment.  Such  a  process  is  time 
consuming  and  may  lead  to  inaccuracies  due  to  source 
changes  during  the  measurement  process.  Franklin  et  al . 

[85]  introduced  the  idea  of  using  photodiode  arrays  mounted 
vertically  in  the  exit  focal  plane  of  the  monochromator 
to  dispense  with  the  need  to  translate  the  measurement 
system  to  obtain  detailed  spatial  profiles.  This  system 
has  been  used  successfully  in  several  studies  involving 
measurement  of  emission  intensity  information  from  the 
ICP.  Edmonds  and  Horlick  [69]  investigated  the  effects 
of  aerosol  flow  rate  and  power  on  spatial  profiles  of  emis¬ 
sion,  and  Horlick  and  Blades  [77]  and  Kawaguiki  et  al.  [71] 
studied  the  effect  of  EIE's  on  analyte  emission  inten¬ 
sities.  Blades  and  Horlick  [83]  have  also  described  a 
photodiode  array  profiling  system  for  implementing  Abel 
inversions  on  emission  from  an  ICP. 

This  approach  has  been  used  to  acquire  all  of  the 
spatial  data  presented  in  this  thesis.  The  measurement 
system  will  be  described  in  detail  in  the  following 
chapter.  The  data  gathered  will  be  discussed  in  sub¬ 
sequent  chapters. 
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CHAPTER  II 


EXPERIMENTAL 


Plasma  Generation 

The  instrumentation  used  for  generating  and  sustaining 
the  plasma  is  outlined  in  Figure  3.  The  system  consists  of 
a  high  power  radio  frequency  (RF)  generator,  an  automatic 
tuning  and  impedance  matching  network,  a  plasma  torch  as¬ 
sembly,  a  nebulizer  and  spray  chamber  for  sample  introduct¬ 
ion,  and  a  source  of  argon  gas  with  associated  gas  pressure 
and  gas  flow  controllers. 

The  RF  generator  (Model  HFS  3000D) ,  automatic  power 
control  (Model  ADC5-3),  automatic  matching  network  (Model 
AMN-2500E) ,  and  plasma  torch  assembly  (Model  PT-2500),  were 
purchased  from  Plasma-Therm  Inc.,  Kresson.  N.J. 

The  RF  generator  is  of  the  crystal-controlled  type 
operating  at  27.12  MHz  with  variable  power  output  from  0  to 
3.0  kw.  The  tuning  network  allows  the  plasma  to  be  ignited 
at  the  proper  impedance  point,  with  automatic  returning  by 
two  variable  vacuum  capacitors,  when  the  plasma  is  es¬ 
tablished.  The  three  turn  RF  coil  is  constructed  from 
1/3"  copper  tubing  and  is  connected  to  water  cooled  copper 
blocks  inside  the  plasma  housing.  The  coil  is  wound  with 
a  26  mm  inside  diameter  and  1  mm  spacing  between  turns. 
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The  distance  between  the  torch  and  coil  is  3  mm. 
continuously  fed  through  the  load  coil  to  keep  it  cool 
during  plasma  operation. 

Plasma  Torch 

The  plasma  torches  used  in  this  study  were  constructed 
at  the  Department  of  Chemistry  glass  blowing  shop  at  the 
University  of  Alberta.  The  torch  consists  of  three  con¬ 
centric  quartz  glass  tubes  fused  to  form  a  single  unit  of 
the  type  described  by  Fassel  et  al.  [11].  A  diagram  of 
the  torch  and  its  dimensions  is  given  in  Figure  4.  The 
outer  quartz  tube,  20  mm  od  and  18  mm  id,  serves  to  confine 
the  argon  gas  during  passage  through  the  load  coils.  The 
middle  quartz  tube  is  tulip  shaped,  12  mm  id  and  14  mm  od 
at  its  exit  and  serves  to  direct  a  low  flow  of  argon  gas  at 
the  base  of  the  plasma  once  it  is  established.  This  po¬ 
sitions  the  plasma  so  that  its  base  is  1  to  2  mm  above  the 
inner  quartz  tube.  The  argon  gas  is  introduced  tangentially 
at  the  base  of  each  tube.  The  inner  or  aerosol  tube  is 
5  mm  od,  3  mm  id,  constricted  to  1.5  mm  id  at  its  exit. 

This  tube  traverses  the  length  of  the  torch  and  is  joined 
to  a  5/12  ball  joint  at  the  torch  base.  This  ball  joint 
couples  with  the  aerosol  spray  chamber  in  the  manner  shown 
in  Figure  5.  The  torch  is  positioned  such  that  the  exit 
of  the  aerosol  tube  is  2-3  mm  below  the  bottom  of  the  load 
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Figure 


IC torch  (Fi.^enoicn^  in  ir) 
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Figure  5. 


ICF  torch,  enray  chamber  assembly. 
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coil . 

Aerosol  Generation 

One  advantage  that  concentric  type  glass  nebulizers 
[86]  have  over  crossflow  nebulizers  [84]  is  that  the  tip 
alignment  remains  fixed  for  the  lifetime  of  the  nebulizer, 
whereas  the  alignment  of  a  crossflow  may  change  during  hookup 
or  cleaning.  The  operation  of  concentric  type  nebulizers 
is,  however,  more  sensitive  to  solutions  containing  high 
salt  concentrations  or  suspended  particles,  such  as  dust. 

If  the  operating  reliability  of  crossflows  under  high  salt 
conditions  is  required,  it  would  be  advantagous  to  have  one 
with  "fixed  tip  alignment",  which  could  be  withdrawn, 
cleaned,  and  re-inserted  into  the  carrier  with  no  change 
in  operating  characteristics. 

For  the  purpose  of  the  studies  carried  out  in  this 
thesis,  such  a  nebulizer  was  constructed,  tested,  and  used 
successfully  for  all  of  the  experiments  described. 

Nebulizer  construction.  Figure  6  is  a  schematic  diagram 
of  the  fixed  tip  crossflow  nebulizer.  The  crossflow  tips, 

A  and  B  are  made  of  1  mm  id,  2  mm  od  quartz  capillaries, 

7-8  cm  long,  which  have  been  drawn  out  in  a  flame;  tip  A  at 
a  right  angle.  The  two  capillaries  are  inserted  into  a 
small  Teflon  jig  with  two  holes  drilled  in  it  to  fix 
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Figure  6.  Schematic  diagram  of  fixed  tip  crossflow  nebulizer 
(dimensions  in  mm) . 
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Figure  7. 


Nebulizer  and  Teflon  carrying  block  (dimensions 
in  mm) . 
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the  centre-to-centre  separation  at  3  mm.  The  tips  are  then 
"ground  down"  using  a  diamond  saw  until  "proper"  tip  over¬ 
lap  and  separation  are  achieved.  Tip  A  is  fire  polished  to 
close  its  tip  diameter  so  that  the  gas  flow  rate  through  it 
is  less  than  1  1pm  at  20  psig  pressure.  While  still  in  the 
jig  the  capillaries  A  and  B  are  connected  to  argon  gas  and 
sample  tubes  respectively  and  the  tips  aligned  to  give  max¬ 
imum  aerosol  production.  They  are  then  bonded  together 
using  E-63  Epoxy  (Techkits,  Demarust,  New  Jersey).  The  tips 
are  subsequently  sealed  into  an  8  mm  od  glass  tube  38  mm 
long  using  RTV  sealant  (Silastic  732  RTV  Adhesive/Sealant , 

Dow  Corning,  Midland,  Mich.)  with  4  to  5  mm  of  the  tips  ex¬ 
posed.  The  nebulizer  may  then  be  inserted  into  a  standard 
Teflon  carrying  block  with  an  8  mm  hole  in  the  side  as  shown 
in  Figure  7. 

Nebulizer  performance.  The  nebulizer  performance  was  tested 
using  a  5  kw  Argon  ICP  obtained  from  Plasma  Therm  Inc.  con¬ 
sisting  of  a  27.12  MHz  RF  generator  (Model  HFS-5000  D) 
with  automatic  power  control  box  (ATPAC) ,  automatic  matching 
network  (AMN  for  HFS-5000) ,  and  plasma  torch  assembly  (PT-1) . 

The  measurement  system  consisted  of  a  Heath  1/3  meter 
monochromator  (GCA/MacPherson  EV-700)  coupled  to  a  Heath 
P.M.T.  module  (EV-701-30)  with  1P28  photomultiplier  tube. 
Emission  intensity  measurements  were  taken  using  a  micro¬ 
processor-based  system,  constructed  concurrent  to  this  study 
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for  the  purpose  of  making  such  measurements.  Details  of  the 
data  acquisition  system  are  included  in  Appendix  A. 

These  results  describe  one  of  five  nebulizers  the 
author  has  constructed  using  the  procedure  outlined.  The 
efficiency  of  the  fixed  tip  crossflow  nebulizer  was  cal¬ 
culated  by  aspirating  a  known  amount  of  distilled  water 
through  the  nebulizer  with  spray  chamber  and  torch  attached 
and  collecting  the  drain-off  liquid.  The  efficiency  meas¬ 
ured  was  3.2%  for  a  sample  uptake  rate  of  8  ml  per  minute, 
when  run  at  0.85  1pm  aerosol  flow  rate  corresponding  to  a 
pressure  of  20  psi.  The  fluid  uptake  rate  of  the  nebulizer 
can  be  controlled  at  assembly  time  by  adjusting  the  size  of 
sample  tip  B.  In  general,  the  larger  the  tip  diameter,  the 
higher  the  sample  uptake  rate  (this  also  degrades  efficiency 
of  the  nebulizer  due  to  leakage  of  the  sample  at  the  tip 
edges) .  One  of  the  results  of  a  sample  uptake  rate  of  8  ml 
per  minute  has  been  that  the  nebulizer  can  be  effectively 
used  for  high  salt  solutions.  To  show  this,  a  1  ppm  Ca 
solution  containing  2%  NaCl  was  nebulized  continuously  for 
1/2  hour.  The  emission  intensity  of  Ca(II)  measured  at 
393.3  nm  varied  no  more  than  6.8%  from  the  average  over 
this  period  of  time.  The  variation  in  signal  appeared 
random,  and  the  intensity  at  the  end  of  the  1/2  hour 
period  was  within  1%  of  the  intensity  at  the  beginning. 

Analytical  performance  of  the  nebulizer  was  compared 
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to  a  conventional  crossflow  nebulizer  of  the  type  described 
by  Kniseley  [84],  and  a  concentric  glass  nebulizer  (Meinhard 
Model  TR-30-A1),  by  aspirating  a  1.0  ppm  solution  of  Ca  and 
integrating  the  emission  intensity  at  393.3  nm.  The  average, 
standard  deviation,  and  signal  to  noise  ratio  were  cal¬ 
culated  on  16  consecutive  1  second  integrations.  The  results 
are  given  in  Table  II. 

The  emission  intensity  was  optimized  for  the  particular 
nebulizer  by  adjusting  the  aerosol  flow  rate  but  otherwise 
all  conditions  were  the  same.  The  relative  intensity  and 
its  standard  deviation  are  shown  in  column  1  and  the  signal 
to  noise  ratio  and  its  standard  deviation  in  column  2. 

These  figures  were  arrived  at  by  averaging  10  sets  of  16, 

1  second  integrations. 

The  fixed  tip  crossflow  nebulizer  constructed  by  the 
author  compares  favourably  with  the  other  two  and  can  be 
used  in  any  application  where  the  other  two  might  be  used. 

The  nebulizer  can  be  constructed  quickly  and  easily  and  once 
made  retains  its  operating  characteristics  for  its  lifetime. 

Spatial  Profiles  With  A  Photodiode  Array 

The  measurement  system  employed  for  emission  profile 
measurements  consists  of  a  Heath  (GCA/MacPherson  EV-700) 

0.35  meter  monochromator  with  a  256  element  self-scanning 
photodiode  array  (Reticon  Model  RL-256/128).  The  photo- 
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TABLE  II 


COMPARISON  OF  CONVENTIONAL  CROSSFLOW,  MEINHARD, 
AND  FIXED  TIP  CROSSFLOW  NEBULIZERS 


Nebulizer  Type 

Relative  Intensity  (a) 

S/N  (a) 

Fixed  Tip 

65.4  (1.52) 

95.5  (25.8) 

Conventional 

70. 0 (1. 16) 

98.2  (13.9) 

Meinhard  Concentric 

47.1 (1.29) 

79.4  (26.0) 
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diode  array  is  mounted  vertically  in  the  exit  focal  plane 
of  the  monochromator.  The  concept  of  spectral  source 
profiling  using  a  photodiode  array  has  been  described  and 
implemented  by  Franklin  et  al.  [85],  and  Edmonds  and 
Horlick  [69].  The  method  is  outlined  in  Figure  8.  An 
image  of  the  source  (the  plasma  in  this  case)  is  formed  on 
the  entrance  slit  of  the  monochromator  using  a  10  cm  focal 
length  bi-convex  quartz  lens.  The  entrance  slit  passes  a 
thin  vertical  slice  of  the  source  image  into  the  monochro¬ 
mator.  The  grating  disperses  the  source  light  into  its 
spectral  components  which  fall  on  the  exit  focal  plane  of 
the  monochromator.  If  one  of  the  spectral  components  is 
allowed  to  fall  on  the  vertical  photodiode  array,  a  readout 
of  the  array  provides  a  quantitative  measurement,  at  the 
wavelength  selected, of  the  image  formed  at  the  entrance 
slit  of  the  monochromator.  The  dimension  of  the  image 
formed  is  altered  by  adjusting  the  distances  between  the 
source  and  lens,  and  the  lens  and  monochromator,  according 
to  the  magnification  formulae  for  bi-convex  lenses. 

Computer  Coupled  Photodiode  Array 

All  studies  were  carried  out  using  a  Reticon  RL-256/128, 
256  element  photodiode  array.  The  individual  photodiodes 
are  25.4  ym  wide  and  are  placed  on  25.4  ym  centres;  thus  the 
light  sensitive  area  of  the  array  is  6.5  mm  in  length,  and 
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25.4  ym  wide.  The  array  and  associated  electronics  are 
housed  in  a  small  metal  cabinet.  The  front  face  plate  of 
the  monochromator  has  been  removed  and  the  cabinet  mounted 
in  its  place.  The  exit  folding  mirror  must  be  removed  to 
allow  spectral  lines  to  fall  on  the  array.  The  array,  when 
mounted  vertically  in  the  exit  focal  plane,  in  a  sense, 
acts  as  the  exit  slit  of  the  monochromator.  The  charac¬ 
teristics  of  photodiode  arrays  and  some  examples  of  their 
implementation  have  been  reviewed  [87-89]. 

A  block  diagram  of  the  computer  controlled  photodiode 
array  readout  system  is  provided  in  Figure  9.  The  array 
provides  two  control  lines,  a  start,  and  a  clock,  which  are 
used  for  timing  and  sequencing  array  signal  acquisition. 

The  start  pulse  signals  the  computer  that  the  array  has 
completed  its  integration  and  initiates  readout  of  in¬ 
formation  from  the  array.  The  time  between  start  pulses 
controls  the  integration  time  and  is  switch  selectable  in 
a  binary  sequence  from  0.0128  to  26.43  sec  when  operating 
at  50  MHz  clock  rate.  The  clock  controls  the  rate  at  which 
data  is  read  from  the  array  and  is  also  used  to  control 
digitization  and  storage  of  the  array  output  signal.  The 
clock  rate  normally  used  is  50  KHz;  thus  a  single  scan  of 
the  array  can  be  acquired  and  stored  in  about  5  msec.  The 
array  is  cooled  to  -10°C  using  a  Peltier  type  cooler  mounted 
on  its  back  face.  The  monochromator  is  purged  with  nitrogen 
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igure  9 


Block  diagram  of  computer  controlled  photodiode 
array  profiling  system. 


gas  to  prevent  icing  of  the  array. 

The  measurement  system  consists  of  a  notch  filter 
(PAR  Model  210  selective  amplifier)  which  is  set  at  1/4  the 
clock  rate  to  remove  the  fixed  pattern  signal  associated 
with  clocking  out  the  array  information.  The  output  from 
this  is  fed  into  a  series  of  operational  amplifiers  which 
serve  to  amplify,  filter,  and  provide  a  voltage  offset  for 
the  array  signal.  A  schematic  diagram  of  this  circuit  is 
provided  in  Figure  10. 

Analog  outputs  from  the  array  are  digitized  using  a 
high-speed  12  bit  analog  to  digital  converter  (Model  600-8E, 
ADAC  Corporation,  Woburn,  Mass.)  which  interfaces  to  the 
PDP8/e  computer  mainframe  and  operates  under  software 
control.  The  computer  is  a  Digital  Equipment  Corporation  PDP8/ 
computer  with  16 k of  core  memory  operating  OS/8  from  a  DEC 
tape-based  system.  The  FORTRAN/machine  language  program 
used  to  acquire  and  manipulate  data  from  the  diode  array  is 
provided  in  Appendix  B.  The  acquired  signal  may  be  stored 
on  DEC  tape,  displayed  on  an  oscilloscope,  plotted  on  an 
XY  recorder  (Hewlett-Packard  Model  7045A) ,  or  output  to  a 
Decwriter . 

The  plasma  torch  assembly,  lens,  and  monochromator- 
photodiode  combination  are  mounted  on  an  instrument  rail 
bed  analogous  to  the  system  designed  and  developed  by 
Walters  [90,91].  This  system  provides  precise,  stable 
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Figure  10.  Signal  conditioning  electronics  for  photodiode  array. 


alignment  while  allowing  flexibility  of  experimental  design. 
Chemicals 


All  solutions  were  made  up  using  reagent  grade  chem¬ 
icals  (Fisher  Scientific  Company)  and  distilled  water.  All 
dilutions  were  performed  using  laboratory  pipettes  and 
volumetric  flasks. 


CHAPTER  III 


VERTICAL  SPATIAL  FEATURES  OF  ANALYTE  EMISSION 

The  establishment  of  ICP-OES  as  a  simultaneous  multi¬ 
element  analysis  method  has  meant  finding  a  set  of  compromise 
conditions  for  operation  of  the  associated  instrumentation 
[46,47,92,93].  Apart  from  fixed  parameters  such  as  RF 
frequency,  method  of  power  coupling,  and  instrumental  dis¬ 
persion  and  detection  system,  there  are  three  major  variable 
parameters.  When  analyte  is  introduced  into  the  plasma  the 
observed  emission  intensities  are  subject  to  the  following 
parameters : 

(1)  RF  power  coupled  into  the  plasma. 

(2)  Aerosol  carrier  gas  flow. 

(3)  Observation  height  and  zone  above  the  load  coil. 

These  parameters  also  affect  the  argon  continuum  and  line 
emission  intensity. 

The  way  in  which  the  variation  of  these  parameters 
affects  net  emission  intensity  was  investigated  by  Boumans 
and  DeBoer  [46,92],  for  the  purpose  of  establishing  com¬ 
promise  conditions  for  simultaneous  multielement  analysis. 

In  these  early  investigations  they  found,  in  general,  that 
emission  intensity  "...increases  with  power ,... decreases 
with  increasing  observation  height ,... depends  to  only  a 
small  extent  on  carrier  gas  flow",  and,  "...is  not  appreciably 
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affected  by  the  addition  of  an  easily  ionizible  element." 
They  also  observed  that  spectral  line  intensities  pass 
through  a  maximum  when  the  carrier  gas  flow  rate  is  in¬ 
creased.  This  effect  was  also  observed  by  Berman  and 
McLaren  [93],  who  noted  that  as  the  aerosol  flow  rate  is 
increased,  the  line  to  background  ratio  increases,  reaches 
a  maximum,  and  then  decreases  over  the  range  0.7  to  1.0 
liters  per  minute  (1pm) .  They  also  observed  that  for  both 
Cr  atom  and  Cr  ion  emission,  net  line  intensity  passed 
through  a  maximum  as  power  was  varied  from  1.0  to  1.8  kw, 
although  atom  lines  peaked  at  a  lower  power  than  ion  lines. 

In  a  later  more  complete  study  Boumans  and  DeBoer 
[42]  observed  the  variation  of  continuum  background  in¬ 
tensity,  net  line  emission  intensity,  and  signal  to  back¬ 
ground  ratio  as  a  function  of  power  input  for  Mnl,  Bal , 

Lil,  Ball,  Cdl ,  Znl,  Mnll,  and  Cdll  species .  They  found 
that  the  net  line  signal  increased  rapidly  initially  and 
then  plateaued  for  Cdll,  Mnl I,  Znl,  and  Cdl  as  the  power 
was  varied  from  0.8  to  2.2  kw,  whereas  the  net  line  signal 
decreased  for  Lil,  Bal,  and  Mnl  over  the  same  power  range. 
Emission  intensity  for  Ball  increased  initially,  then  re¬ 
mained  unchanged  after  1.2  kw.  These  observations  were  for 
a  5  mm  high  observation  window  centered  at  15  mm  above  the 
load  coil.  They  related  this  behaviour  to  the  "norm 
temperature"  of  the  lines,  which  is  in  turn  related  to  the 
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ionization  and  excitation  potentials  of  the  species  in¬ 
volved.  Boumans  later  termed  lines  which  showed  an  increase 
in  intensity  with  an  increase  in  power  as  "hard"  lines 
and  those  which  showed  a  decrease  in  intensity  with  an  in¬ 
crease  in  power  as  "soft"  lines  [94],  A  full  treatment  and 
discussion  of  the  importance  of  "norm  temperature"  to  spatial 
structure  of  analyte  emission  will  be  given  in  the  following 
chapter . 

An  interpretation  of  the  response  of  line  intensity  to 
changes  in  RF  power  and  aerosol  flow  rate  is  provided  in 
spatial  studies  carried  out  by  Edmonds  and  Horlick  [69]  and 
Dube  and  Boulos  [68].  Both  groups  acquired  analyte  emission 
intensity  information  as  a  function  of  vertical  height  above 
the  load  coil.  Both  groups  found  that  the  emission  in¬ 
tensity  profiles  for  analyte  species  were  in  the  general 
shape  of  a  "camel's  hump".  That  is,  as  one  observes  em¬ 
ission  intensity  up  the  analyte  channel  it  initially  in¬ 
creases,  reaches  a  maximum  and  then  decreases  forming  an 
emission  intensity  peak  or  maximum.  The  position  of  this 
peak  maximum  was  found  to  shift  upward  as  aerosol  flow  rate 
was  increased.  This  suggests  that  the  emission  intensity 
maximum  observed  by  Boumans  and  DeBoer  as  a  function  of 
aerosol  flow  rate  is  the  result  of  an  upward  shift  in  the 
spatial  profile  across  the  spatial  observation  window. 

Edmonds  and  Horlick  further  found  that  analyte  spatial  em- 
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ission  profiles  fell  into  two  distinct  categories  based  on 
their  response  to  applied  power  at  a  fixed  aerosol  flow 
rate.  The  intensity  maximum  of  the  atom  lines  of  Ca,  Ba, 

Sr,  Cr ,  and  Ti  shifted  to  a  lower  position  in  the  plasma  as 
power  was  increased  from  1.25  kw  to  2  kw,  whereas  the 
intensity  maximum  of  the  ion  lines  of  Ba,  Sr,  Ti ,  and  Ca 
did  not  shift  as  power  was  varied  over  the  same  range.  They 
concluded  that  atom  lines  exhibited  one  type  of  spatial 
dependence  on  power  and  ion  lines  a  different  type  of 
spatial  dependence.  The  atom  lines  of  Cu  and  Ag ,  however, 
were  observed  to  behave  like  "ion"  lines. 

A  correlation  between  the  observations  of  Boumans  and 
DeBoer  and  Edmonds  and  Horlick  can  be  made  by  examining 
the  contents  of  Figures  11,  12,  and  13.  Figures  11  and  12 
are  similar  to  the  type  of  data  published  previously  [69]. 

The  plots  given  in  Figure  11  are  vertical  spatial  profiles 
of  Ca  emission  from  the  422.7  nm  atom  line.  The  frames 
(a),  (b)  and,  (c)  relate  to  aerosol  flow  rates  of  0.7,  0.8, 

and  0.9  1pm  respectively.  The  numbers  from  1  through  5  on 
each  frame  represent  applied  powers;  these  are  :  (1)  1.0  kw, 

(2)  1.25  kw,  (3)  1.50  kw,  (4)  1.75  kw,  and  (5)  2.0  kw.  The 
emission  intensity  is  plotted  from  0  mm  to  16  mm  above  the 
load  coil,  centered  on  the  analyte  channel.  In  agreement 
with  the  observations  of  Edmonds  and  Horlick,  as  the  RF 
power  is  increased,  the  position  of  the  maximum  in  em- 
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Figure  11.  Vertical  spatial  profiles  for  Cal  422.7  nm.  Aerosol  flow  rates;  (a)  0.7  1pm 

(b)  0.9  lpm,  (c)  0.9  1pm  RF  power;  (1)  1.0  kw,  (2)  1.25  kw ,  (3)  1.50  kw, 

(4)  1.75  kw,  and  (5)  2.0  kw. 
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Figure  12.  Vertical  spatial  priflies  for  Call  393.3  nru.  Aerosol  flew  rates;  (a)  0.7  1pm 

(h)  0.8  lnm ,  (c)  0.9  1pm,  RF  power;  (1)  1.0  kw,  (2)  1.25  kw,  (3)  1.50  kw, 

(4)  1.75  kw,  and  (5)  2.0  kw. 


ission  intensity  shifts  from  higher  to  lower  in  the  aerosol 
channel.  Similar  behaviour  has  been  observed  by  Kawaguchi 
et  al.  [71],  and  Savage  and  Hieftje  [96].  This  is  dem¬ 
onstrated  in  the  sequence  in  frame  (a)  going  from  the  plot 
labelled  1  through  to  the  plot  labelled  5.  In  addition,  the 
net  intensity  of  the  peak  increases  from  1.0  to  1.25  kw, 
and  then  decreases  from  1.25  kw  up  to  2.0  kw.  At  a  fixed 
power  (for  example  for  curve  2)  an  increase  in  aerosol  flow 
rate  shifts  the  peak  maximum  to  a  higher  position  above  the 
load  coil  and  simultaneously  decreases  the  net  emission 
intensity  of  the  peak. 

Analogous  plots  for  the  emission  from  the  ion  line  of 
Ca  at  393.3  nm  are  provided  in  Figure  12,  in  the  region  6 
to  22  mm  above  the  load  coil.  The  ion  line  peaks  at  higher 
positions  above  the  load  coil  and  the  maximum  does  not 
show  the  same  shift  with  changes  in  applied  power  that  is 
observed  for  the  atom  line.  Instead,  as  the  RF  power  is 
increased  from  1.0  to  2.0  kw,  the  net  intensity  increases, 
while  maintaining  approximately  the  same  spatial  profile 
height  dependence. 

If  one  were  to  observe  the  emission  intensity  from 
Cal  and  Call  with  an  instrument  with  a  fixed  observation 
height  and  relatively  small  vertical  entrance  slit  (such 
that  a  narrow  spatial  region  of  the  plasma  was  observed) 
we  would  observe  the  data  presented  in  Figure  13.  This 
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Figure  13.  Emission  intensity  of  Cal  and  Call  as  a  function 
of  RF  Power  at  an  observation  height  of  14  mm. 
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figure  provides  a  plot  of  emission  intensity  at  a  height 
of  14  mm  above  the  load  coil  for  Cal  and  Call  obtained 
from  Figures  11  and  12  for  various  RF  powers.  It  can  be 
seen  that  the  intensity  of  Cal  decreases  with  applied  power, 
and  the  intensity  of  Call  increases  with  applied  power 
plateauing  at  higher  levels  of  power.  This  is  precisely 
the  type  of  differing  behaviour  described  by  Boumans  and 
DeBoer  for  the  "soft"  and  "hard"  lines  respectively.  In 
all  likelihood  their  observations  were  probably  the  result 
of  the  differing  spatial  behaviour  of  the  different  species. 
Thus  the  division  of  line  behaviour  into  the  classes  of 
hard  and  soft  can  be  based  on  the  spatial  profile  behaviour 
in  response  to  changes  in  the  applied  power.  Soft  lines 
show  the  type  of  spatial  behaviour  exhibited  by  Cal, 
whereas  hard  lines  are  typified  by  the  type  of  spatial 
behaviour  of  Call.  It  should  be  noted  that  the  relative 
intensities  of  the  two  lines  are  not  the  same  and  vary  with 
observation  height.  This  information  is  not  apparent  from 
the  plots  provided.  For  example,  at  a  height  of  6  mm  the 
ratio  of  Call  emission  to  Cal  emission  is  1.54;  at  20  mm 
this  ratio  is  15.5,  at  0.85  1pm  aerosol  flow  rate  and  1.5 
kw  RF  power. 

On  the  basis  of  the  previous  discussion  we  can  refer 
back  to  the  plots  of  Boumans  and  DeBoer [42]  and  predict  that 
the  670.8  nm  Li  atom  line  will  show  marked  spatial  variation 
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with  power,  whereas  Cdl  (228.8  nm) ,  Znl  (213.9  nm) ,  Mnll 
(403.1  nm)  and  Cdll  (226.5  nm)  will  not.  The  accuracy  of 
this  prediction  is  reflected  in  the  spatial  profiles  pro¬ 
vided  in  Figures  14,  15,  16,  17,  and  18  for  the  above  lines. 
Lil  plotted  in  Figure  14  shows  spatial  variation  in  the 
region  4  to  18  mm  above  the  load  coil  as  power  is  changed 
from  1.25  kw  to  2.0  in  0.25  kw  increments  corresponding  to 
plots  2  to  4.  These  were  taken  at  0.85  1pm  aerosol  flow 
rate.  On  the  other  hand,  Cdl,  Cdll,  Mnll  and  Znl  do  not 
show  spatial  shifts  with  changes  in  RF  power.  For  the  sake 
of  completeness  the  spatial  variation  of  emission  intensity 
with  power  and  flow  rate  have  been  included  for  Cdl,  Cdll, 
Mnll,  and  Znl.  Also,  in  the  plot  for  Znll  provided  in 
Figure  19,  Znll  is  seen  to  behave  like  a  hard  line.  For 
the  plots  given  in  Figures  15  to  19,  the  frames  (a),  (b) , 
and  (c)  represent  aerosol  flow  rates  of  0.9,  1.0,  and  1.1 
lpm  respectively.  The  numbers  on  each  frame  represent  the 
power  applied  in  the  following  manner:  (1)  1.25  kw,  (2) 

1.50  kw,  (3)  1.75  kw,  (4)  2.0  kw,  and  (5)  2.25  kw. 

In  all  cases,  the  "hard"  lines  are  seen  to  peak  at 
about  18-20  mm  above  the  load  coil,  and  to  decrease  in  net 
intensity  with  a  decrease  in  applied  RF  power,  but  not 
shift  in  spatial  position.  The  effect  of  an  increase  in 
aerosol  flow  rate  is  to  shift  the  emission  intensity  peak 
only  slightly  higher  in  spatial  position,  while  simul- 
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Figure  14.  Vertical  spatial  profiles  for  Lil  670.8  nm 
power;  (1)  1.25  kw,  (2)  1.50  kw,  (3)  1.75 
and  (4)  2.0  kw. 
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Figure  17.  Vertical  spatial  profiles  for  Mnll  403.1  nm.  Aerosol  flow  rates;  (a)  0.9  1pm 

(b)  1.0  lpm ,  and  (c)  1.1  1pm.  RF  power;  (1)  1.25  kw,  (2)  1.50  kw,  (3)  1.75 
kw,  (4)  2.0  kw,  and  (5)  2.25  kw. 
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Figure  18.  Vertical  special  profiles  for  CdH  226.5  nm.  Aerosol  flow  rates;  (a)  0.9  1pm 

(b)  1.0  1pm,  and  (c)  1.1  1pm.  RF  power;  (1)  1.25  kw,  (2)  1.50  kw,  (3)  1.75 
kw,  (4)  2.0  kw,  and  (5)  2.25  kw. 
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taniously  reducing  the  overall  emission  intensity. 

It  may  be  noted  that  while  the  emission  intensity  for 
hard  type  lines  is  expected  to  increase  with  an  increase  in 
applied  power  regardless  of  the  spatial  observation  zone, 
the  behaviour  for  soft  lines  is  not  quite  so  predictable. 
For  example,  the  spatial  profiles  of  Nal  588.9  nm, 
a  soft  line,  are  provided  in  Figure  20.  The  aerosol  flow 
rate  is  0.8  1pm  for  curve  (a),  and  0.9  1pm  for  curve  (b) . 
The  format  of  the  plots  is  analogous  to  those  previously 
discussed.  Curves  1,  2,  3,  and  4  represent  1.25  kw,  1.5 
kw,  1.75  kw,  and  2.0  kw  respectively.  A  plot  of  the  net 
emission  intensity  versus  power  is  provided  in  Figure  21 
for  observation  heights  of  13  mm,  19  mm,  and  25  mm  above 
the  load  coil.  The  figure  demonstrates  three  different 
types  of  behaviour  for  the  Na  atom  emission  intensity  in 
response  to  applied  RF  power  depending  on  the  height  of 
observation.  At  25  mm  the  intensity  decreases  with  an 
increase  in  power,  at  19  mm  the  intensity  has  a  maximum 
at  1.5  kw,  and  at  13  mm  the  intensity  increases  with  power. 
Clearly  the  spatial  profile  information  is  important  in 
characterizing  the  nature  of  the  response  of  emission 
intensity  to  changes  in  ICP  operating  parameters. 

The  vertical  spatial  profiles  of  the  418.2  nm  Ar  atom 
line  are  provided  in  Figure  22  in  the  region  10  to  28  mm 
above  the  load  coil.  The  two  frames  (a)  and  (b)  are  for 
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Figure  21.  Emission  intensity  as  a  function  of  power  for  Nal 
588.9  nm,  at  observation  heights  of  13  mm,  19  mm, 
and  25  mm  above  the  load  coil. 
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Figure  22.  Vertical  spatial  profiles  for  Arl  418.2  nm.  Aerosol  flow  rate;  (a)  0.9  1pm 

(b)  1.0  1pm.  RF  power;  (1)  1.0  kw,  (2)  1.25  kw,  (3)  1.50  kw,  (4)  1.75  kw, 

and  (5)  2.0  kw. 
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argon  aerosol  flow  rates  of  0.9  1pm  and  1.0  1pm  respectively. 
The  plots  from  1  to  5  represent  applied  powers  of  1.0  to 
2.0  kw  in  increments  of  0.25  kw.  The  argon  line  intensity 
is  seen  to  decrease  with  a  decrease  in  RF  power  as  would 
be  expected,  and  shows  little  spatial  response  to  changes 
in  aerosol  flow  rate.  The  emission  intensity  decreases  by 
a  factor  of  about  0.75  at  its  maximum  as  the  aerosol  flow 
rate  increases.  This  is  the  result  of  cooling  of  the  aerosol 
channel  leading  to  a  decrease  in  argon  emission  from  this 
region . 

In  summary,  the  following  features  characterize  spatial 
behaviour  of  analyte  species  in  the  analyte  channel: 

(1)  Emission  lines  may  be  subdivided  into  two  categories 
based  on  their  spatial  profile  response  to  changes  in 
applied  power.  These  two  categories  have  been  given 
the  names  "soft"  and  "hard"  following  the  terminology 
of  Boumans  [ 94 ] . 

(2)  The  vertical  spatial  profiles  of  soft  lines  are  much 
more  sensitive  to  changes  in  power  than  the  spatial 
profiles  of  hard  lines. 

(3)  For  soft  lines  an  increase  in  RF  power  causes  the 

emission  intensity  peak  to  shift  lower  in  the  plasma. 

For  hard  lines  an  increase  in  RF  power  causes  very 

little  shift  in  the  position  of  maximum  emission. 
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(4)  In  all  cases  ion  lines  have  demonstrated  hard  line 
behaviour.  Atom  lines  may  display  hard  or  soft  line 
behaviour,  depending  upon  the  species  being  observed. 

(5)  Hard  lines  are  observed  to  peak  at  positions  higher 

above  the  load  coil  than  soft  lines.  This  is  in  agree¬ 
ment  with  the  observations  of  Edmonds  and  Horlick  [69] , 

Kawaguchi  et  al.  [71],  and  Savage  and  Hieftje  [96]. 

The  question  which  arises  is:  "What  is  the  cause  for 
the  two  different  types  of  spatial  behaviour  for  analyte 
species  in  the  ICP?"  Edmonds  and  Horlick  have  roughly 
correlated  the  spatial  position  of  maximum  intensity  with 

the  excitation  potential  of  the  line  being  observed.  This, 

however,  cannot  be  the  only  criterion  for  determining  spatial 
behaviour  since  Cal  (422.7  nm)  and  Call  (393.3  nm)  have 
similar  excitation  potentials  (2.9  ev  and  3.1  ev  respectively) . 
The  process  of  ionization  of  atom  species  must  play  an 
important  role  in  determining  spatial  behaviour  since  it 
will  define  the  relative  population  of  atom  and  ion  species 
in  the  ICP.  If  the  excitation  and  ionization  characteristics 
of  the  line  being  observed  are  important  in  determining  the 
spatial  position  of  the  atom  and  ion  lines,  how  can  we  use 
these  characteristics  to  predict  spatial  behaviour?  These 
questions  are  addressed  in  the  next  chapter  of  this  thesis. 

In  a  more  practical  vein;  what  use  can  be  made  of  the 
spatial  information  presented  in  this  chapter?  The  use  of 


internal  standards  in  emission  spectroscopy  is  well  known. 
Internal  standards  are  chosen  such  that  they  will,  as 
closely  as  possible,  mimic  analyte  emission  behaviour  in 
response  to  changes  in  various  parameters  which  will  affect 
the  emission  intensity  of  the  analyte.  The  theoretical 
aspects  of  choosing  internal  standards  have  been  discussed 
by  Barnett  et  al.  [95] .  In  the  ICP  the  use  of  the  internal 
standard  principle  can  improve  precision  [97,98].  The 
spatial  behaviour  of  various  analyte  species  in  response 
to  applied  power  has  been  discussed  in  this  chapter,  and  it 
has  been  shown  that  different  species  can  behave  quite 
f srently .  This  spatial  behaviour  should  certainly  be 
taken  into  account  in  choosing  an  appropriate  internal 
standard.  For  example,  a  soft  line  should  not  be  chosen 
as  a  ratio  partner  with  a  hard  line.  Of  course,  the  choice 
of  an  internal  standard  with  similar  excitation  character¬ 
istics  as  the  analyte  will  probably  assure  similar  spatial 
behaviour . 

In  a  simultaneous  multielement  analysis  situation, 
information  such  as  that  presented  in  this  chapter  should 
be  considered  carefully.  The  procedure  of  parameter  op¬ 
timization  could  be  considerably  simplified  with  the  aid  of 
this  spatial  information. 
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CHAPTER  IV 


SPATIAL  BEHAVIOUR  AND  ANALYTE  EXCITATION  CHARACTERISTICS 
Introduction 

In  the  previous  chapter  it  was  pointed  out  that  the 
spatial  behaviour  of  analyte  emission  may  be  strongly  spe¬ 
cies  dependent.  A  trend  emerged  in  which  the  response  of 
the  spatial  profile  to  changes  in  applied  RF  power  divided 
into  two  classes;  hard  or  soft,  depending  on  the  analyte 
species  being  observed.  The  purpose  of  this  chapter  is 
to  investigate  the  manner  in  which  analyte  ionization  and 
excitation  characteristics  determine  spatial  behaviour  of 
analyte  emission  in  the  ICP. 

Soft  Lines 

Edmonds  and  Horlick  [69]  demonstrated  a  rough  cor¬ 
relation  between  the  position  of  the  intensity  maximum 
in  the  vertical  spatial  profile  and  analyte  excitation 
characteristics  for  atom  lines  (soft  lines)  in  the  region 
12  mm  to  24  mm  above  the  load  coil.  Atom  lines  with  low 
excitation  potential  tended  to  peak  lower  in  the  analyte 
channel  and  those  with  higher  excitation  potential  tended 
to  peak  higher  up.  In  an  effort  to  more  clearly  define  the 
nature  of  this  spatial  structure  a  similar  study  was  carried 
out  on  a  set  of  atom  lines  with  excitation  potentials 
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ranging  from  1.6  ev  to  4.4  ev.  The  spatial  emission 
profiles  for  the  various  lines  in  the  region  4  to  18 
mm  above  the  load  coil  are  provided  in  Figure  23,  which 
plots  "emission  intensity"  versus  "height  above  the  load 
coil"  for  these  lines.  Some  of  the  profiles  have  been 
truncated  and  offset  in  the  vertical  direction  for  clarity. 

The  conditions  chosen  for  the  study  were  1.4  kw  RF  power, 

0.85  1pm  aerosol  flow  rate,  15  1pm  plasma  gas  flow  rate,  and 
less  than  1  1pm  of  auxiliary  gas.  These  conditions  are  close 
to  those  normally  used  for  quantitative  multielement  analysis. 
The  analyte  lines  studied  are  listed  in  Table  III  along  with 
their  respective  excitation  potentials  and  observed  position 
of  peak  maximum.  The  latter  was  measured  by  locating  the 
position  where  the  emission  intensity  is  greatest.  The 
position  of  maximum  intensity  is  plotted  against  excitation 
potential  in  Figure  24.  The  plot  indicates  a  crude  cor¬ 
relation  between  excitation  potential  and  position  of  peak 
maximum,  however  there  are  two  features  of  the  plot  which 
indicate  that  excitation  potential  is  not  the  sole  criterion 
for  determining  the  peak  position.  The  first  of  these  is 
the  fact  that  Na,  K,  and  Li  (670.8  nm)  form  one  straight 
line,  and  Ca,  Cr,  and  Mg  form  another  of  differing  slope. 
Secondly,  the  Li  line  at  610.4  nm  does  not  fall  on  either 
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Height  Above  Load  Coil  (mm) 


Figure  23.  Vertical  spatial  emission  profiles  for;  (a)  Nal 
588.9  nm,  (b)  Lil  610.4  nm,  (c)  Lil  670.8  nm, 

(d)  CrI  357.8  nm,  (e)  KI  766.5  nm,  (f)  Cal 
422.7  nm,  (5)  Mg I  285.3  nm. 


TABLE  III 


EXCITATION  POTENTIALS  AND  PEAK  MAXIMUM  POSITION 
FOR  ANALYTE  LINES  IN  FIGURE  23 


Excitation 

Element  Line  (nm)  Potential  (ev)  Peak  Position 


KI 

766.5 

1.62 

8.9 

Lil 

570.8 

1.85 

10.5 

Lil 

610.4 

3.87 

11.3 

Nal 

588.9 

2 . 11 

11.8 

Cal 

422.7 

2.93 

11.8 

CrI 

357.8 

3.46 

12.7 

Mg  I 

285.2 

4 . 34 

14.4 

* 

mm 

above  load  coil 

Position  of  peak  maximum  for  lines  listed  in 
Table  III,  as  a  function  of  excitation  potential 
of  the  line. 


Figure  24. 
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line  but,  rather,  appears  as  an  isolated  data  point. 

In  order  to  more  accurately  determine  which  analyte 
and  plasma  characteristics  are  important  in  determining 
peak  position  the  problem  will  be  approached  from  a  theor¬ 
etical  standpoint.  This  theoretical  approach  will  include 
analyte  parameters  such  as  analyte  excitation  potential 
and  ionization  potential,  as  well  as  ambient  channel 
temperature  ,  a  plasma  parameter.  The  basis  of  the  ap¬ 
proach  is  the  assumption  that  temperature  increases  as  one 
moves  from  a  position  just  above  the  load  coil  up  through 
the  analyte  channel  away  from  the  load  coil.  This  assumption 
has  not  been  based  on  measured  values  of  spatial  temperature 
distributions  which  have  appeared  in  the  literature  but 
rather  on  an  intuitive  feeling  that  the  analyte  spatial 
structure  suggests  this.  Atomic  species  which  are  easily 
excited  consistently  peak  lower  in  the  analyte  channel  than 
those  which  are  more  difficult  to  excite.  In  addition  all 
ion  lines  observed  thus  far  have  peaked  high  in  the  analyte 
channel;  this  in  spite  of  the  fact  that  ionization-exci¬ 
tation  for  these  species  is  a  very  energetic  process.  It  is 
acknowledged  that  some  of  these  observations  may  be  the 
result  of  "non-thermal"  processes  (such  as  Penning  ion¬ 
ization)  ,  however,  for  the  sake  of  simplicity  this  has  not 
been  considered  in  the  following  treatment.  This  point  will 
be  returned  to  in  the  discussion  at  the  end  of  this  chapter. 
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Theoretical  Treatment  of  Analyte  Emission  Intensity 

The  processes  of  desolvation,  vaporization,  and  dis¬ 
sociation  of  the  analyte  in  the  region  of  the  load  coil  result 
in  ground  state  atomic  species  traveling  up  the  analyte 
channel.  The  temperature  increase  experienced  by  the  ground 
state  atom  will  influence  the  emission  intensity  observed 
from  the  atomic  species  in  two  opposing  ways.  As  the  atom 
experiences  a  steadily  increasing  temperature  the  emission 
intensity  will  increase  as  the  excited  state  population 
increases  according  to  the  Boltzmann  distribution.  Opposing 
this  is  a  decrease  in  population  of  atom  species  due  to 
increased  ionization  as  described  by  the  Saha  equation. 

This  will  tend  to  decrease  emission  intensity  of  atom 
species.  At  some  point  in  the  analyte ' s  progress ,  the  effects 
of  ionization  will  overtake  the  effects  of  increased  ex¬ 
cited  state  population  and  the  emission  intensity  will 
reach  a  peak  and  begin  to  decrease,  forming  an  emission 
intensity  maximum  at  some  position  above  the  load  coil.  The 
temperature  at  which  this  maximum  occurs  has  been  called 
the  "normal  temperature"  [94]. 

This  process  may  be  quantitatively  described  in  the 
following  way.  The  observed  emission  intensity  of  an  atomic 
species,  a,  in  excited  state  1  is  given  by: 


(24) 


, 

. 


4-i 


where  h  is  Planck's  constant/  v  is  the  frequency  of  the 

al 

transition,  A  is  the  transition  probability  and  N  is 

1  al 

the  number  of  atoms  in  the  excited  state.  The  number 

of  atoms  in  the  excited  state  can  be  described  by  the 

Boltzmann  distribution: 


N 


Na  9ai 

w 


exp  (-E  /kT  ) 
a^  a 


(25) 


where  N  is  the  total  number  of  atomic  species,  g  is 

al 

the  statistical  weight  of  the  excited  state,  Q  (T  )  is 

a  a 

the  electronic  partition  function,  E  is  the  excitation  po 

al 

tential  of  the  excited  state,  k  is  the  Boltzmann  constant, 
an|3  Ta  atom;*-c  excitation  temperature.  The  total  num 

ber  of  atoms,  N  ,  is  given  by: 


N  =  N  ( 1-X ) 

a.  L. 


(26) 


where  N  is  the  total  number  of  analyte  species  present  in 
all  forms,  ionic  and  atomic,  and  X  is  the  degree  of  ioni¬ 
zation.  Substituting  the  expression  given  in  Equations 
25  and  26  into  Equation  24  yields; 


g  A 

Si  cl 

Ia,  "  hv  Nt(l-X)  Q-^jry  exp(-E  /kTa>  (27) 

■!  i.  a  a  l 


This  expression  gives  the  emission  intensity  in  terms 
of  the  degree  of  ionization  of  the  analyte,  its  excitation 


‘ 
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potential,  and  the  temperature  experienced  by  the  atom  such 
that : 


I 


(1-X) 

Qa(Ta> 


exp ( -E  /kT  ) 
a^  a 


(28) 


if  the  total  number  of  atoms  N  is  constant. 

The  ionization  can  be  expressed  in  terms  of  electron 

pressure,  pg ,  and  the  Saha  equilibrium  constant,  S  ,  by  [56]: 

P 


(1-X) 


p  +-S 
e  p 


(29) 


where : 


-7  5/2  ^-j  ) 

sp  =  6.58  X  10  Ti  7  Q-jijry  exp  (-V ,/kT..  )  (30) 

a  a 

and  T±  is  the  ion  temperature,  Q± (T± )  is  the  ion  partition 
function,  and  is  the  ionization  potential  of  the  atom. 

The  partial  pressure  of  electrons  is  related  to  electron 
density  (ng)  by; 

pe  =  1.36  x  10'22  neT  (31) 

For  the  sake  of  simplicity,  it  is  assumed  that  LTE 
exists,  that  is,  the  temperature  defining  the  populations  of 
excited  states  and  the  temperature  defining  the  relative 
distribution  of  atom  and  ion  species  is  the  same.  Under 
this  condition  the  relative  emission  intensity  of  atom  a 


30 


from  excited  state  1  can  be  expressed  as; 


exp  (-E  /kT) 
al 


(32) 


where  T  is  the  LTE  temperature. 

Equation  32  may  be  solved  to  yield  relative  emission 


intensity  for  an  atom  line  at  various  values  of  temperature. 


The  logarithm  of  the  partition  functions  Q  (T)  and  Q.(T)  are 

3  1 

assumed  to  vary  linearly  with  temperature  [56];  according  to 
the  formula; 


log  Q  =  log  Q'  -  (5|±°)  £ 


(33) 


where  Q'  and  £  are  constants  for  a  particular  atom  or  ion  and 
Q  is  the  partition  function  at  temperature  T.  The  values 
of  Q'  and  £  for  the  set  of  atom  lines  studied  are  listed  in 
Table  IV  for  atom  (subscript  a)  and  ion  (subscript  i) . 

These  values  have  been  extracted  from  Boumans  [56] . 

The  procedure  employed  was  to  solve  Equations  30,  31, 
and  33  for  a  series  of  temperature  values  and  inserting  the 
solutions  into  Equation  32  to  arrive  at  relative  emission  in¬ 
tensities  for  these  temperatures.  A  plot  of  emission 
intensity  as  a  function  of  temperature  will  yield  a  result 
such  as  that  provided  in  Figure  25  for  KI  (766.5  nm)  and 
Mg I  (285.3  nm)  for  an  electron  density  of  10^*5  cm~^.  The 


overall  shape  and  relative  position  of  these  curves  are 
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TABLE  IV 


Element 

ATOM  AND  ION 

Q'a 

PARTITION  FUNCTION 

Q'i 

CONSTANTS 

£ 

£  . 

K 

4.374 

1.0 

0.30 

1 

0.0 

Li 

2.636 

1.0 

0.10 

0.0 

Ba 

22.946 

9.463 

0.95 

0.35 

Na 

2.586 

1.0 

0.10 

0.0 

Sr 

4.951 

3.0595 

0.6 

0.15 

Ca 

3.325 

3.516 

0.45 

0 .20 

Cr 

32.873 

17.966 

0.50 

0.40 

Mn 

10.181 

12.249 

0.20 

0.20 

A1 

5.9 

1.0 

0.0 

0.30 

Mg 

1.146 

2 . 0 

0.05 

0.0 

Cd 

1.0 

2.0 

0 . 0 

0 . 0 

Co 

67.211 

59.387 

0.30 

0 . 30 

. 
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Figure  25.  Theoretically  generated  profiles  of  emission  intensity  as  a  function  of  height 
above  the  load  coil  for  Mgl  285.3  nm  and  KI  766.5  nm. 
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Mg  (I) 
285.3  nm 


K  (I) 

766.5  nm 


Figure  26. 


Vertical  spatial  profiles  for  Mg I  285.3  nm  and 
KI  766.5  nm. 


8  4 


markably  similar  to  those  which  are  obtained  experimentally 
from  the  ICP  as  can  be  seen  by  comparison  of  Figures  25  and 
26.  The  curves  shown  in  Figure  26  are  plots  of  relative 
emission  intensity  as  a  function  of  height  above  the  load 
coil  for  these  same  K  and  Mg  lines,  for  an  RF  power  of  1.4  kw, 
aerosol  flow  rate  of  0.8  1pm,  and  plasma  gas  flow  rate  of 
15  1pm.  This  indicates  a  strong  correlation  between  the 
temperature  axis  in  Figure  25  and  the  height  above  the  load 
coil  axis  in  Figure  26.  Both  the  position  of  the  peak 
maximum  and  the  profile  width  are  in  theory  dependent  on 
temperature,  and  in  practice  dependent  on  the  vertical 
position  above  the  load  coil.  As  was  mentioned  previously, 
the  temperature  at  which  the  maximum  in  emission  intensity 
occurs  in  Figure  25  is  called  the  normal  temperature.  The 
comparison  of  Figures  25  and  26  indicates  that  the  spatial 
position  of  the  peak  maximum  observed  in  the  ICP  is  probably 
related  to  the  normal  temperature  of  the  element  being  ob¬ 
served.  This  proposal  is  further  substantiated  by  the 
plot  provided  in  Figure  27.  In  this  figure  the  normal 
temperature  for  the  lines  listed  in  Table  III  are  plotted  as 
a  function  of  the  excitation  potential  of  the  line.  The 
electron  density  assumed  is  10  cm  .  This  plot  is 

analogous  to  that  provided  in  Figure  24.  It  is  clear  from 
these  two  plots  that  normal  temperature  and  peak  position 
are  strongly  related  to  each  other. 
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Figure  27 


Plot  of  normal  temperature  versus  excitation 
potential  for  lines  listed  in  Table  III. 
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In  order  to  verify  this  correlation  of  peak  position 
and  calculated  normal  temperature,  the  normal  temperature 
was  calculated  for  a  set  of  16  atom  lines  of  varying  ex¬ 
citation  and  ionization  potentials.  These  are  listed  in 
Table  V.  The  electron  density  was  taken  as  1013'5  cm"3  for 
this  calculation.  Kalnicky  et  al.  [36]  have  recently  mea¬ 
sured  the  electron  density  using  the  Saha  equation  to  be 
between  10  cm  and  10  cm  in  the  region  15  mm  to  25  mm 
above  the  load  coil.  The  channel  electron  density  demon¬ 
strated  a  height  dependence  decreasing  with  increasing  ob¬ 
servation  height.  For  simplicity  a  uniform  value  of  1013*5 
-3 

cm  was  adopted  for  the  values  of  normal  temperature  cal¬ 
culated.  The  effect  of  an  increase  in  electron  density  is 
to  increase  the  normal  temperature.  A  plot  of  normal  temp¬ 
erature  as  a  function  of  electron  density  is  provided  in 
Figure  28  for  Mgl ,  CrI,  Bal ,  and  KI.  It  can  be  seen  that 
regardless  of  the  electron  density,  the  normal  temperature 
for  the  lines  will  be  proportionally  related  to  each  other. 

The  plot  provided  in  Figure  29  demonstrates  that  the 
emission  peaks  observed  in  the  ICP  are  strongly  related  to 
the  normal  temperature  of  the  line  being  observed.  In  this 
figure  the  peak  position,  a  value  obtained  from  experiment 
is  plotted  as  a  function  of  the  calculated  normal  temperature 
for  the  lines  tabulated  in  Table  V.  The  spatial  position  is 
plotted  with  a  bar  to  represent  the  finite  width  of  the 
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Figure  28.  Normal  temperature  for  Mgl  285.3  nm,  CrI  425.4  nm,  Bal  553. 
nm,  as  a  function  of  electron  density. 
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TABLE  V 

NORMAL  TEMPERATURE  FOR  SEVERAL  ANALYTE  LINES 

Normal 


Species 

Line 

Tempe 

(nm) 

(K) 

K 

766 . 5 

2870 

Li 

670.1 

3578 

Ba 

553 . 4 

3100 

Li 

610 . 3 

3878 

Na 

588 . 9 

3470 

Sr 

460.7 

3510 

Ca 

422 . 7 

3801 

Cr 

425.4 

4372 

Mn 

403.0 

4671 

Al 

396.2 

4465 

Cr 

357 . 8 

4448 

Co 

345.8 

5078 

Mg 

383.8 

5048 

Mg 

285.3 

4837 

Mn 

279.4 

4838 

Cd 

228.8 

5655 

' 


. 

Mn  279.4 
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spatial  profile  at  its  peak.  This  plot  confirms  the  notion 
that  the  peaks  in  the  spatial  profiles  of  these  lines  are 
the  result  of  the  excitation-ionization  competition  which 
results  in  a  normal  temperature  maximum.  The  plot  is  a 
rather  indirect  way  of  arriving  at  an  LTE  temperature  pro¬ 
file  of  the  aerosol  channel  in  the  region  7  mm  to  16  mm 
above  the  load  coil.  The  temperatures  in  this  region  are 
seen  to  range  from  3000  K  to  6000  K.  These  are  consistent 
with  the  values  of  temperature  which  have  been  measured  in 
our  own  lab  and  by  other  workers  [36,38].  Decker  [100]  has 
used  a  similar  treatment  to  rationalize  spatial  behaviour  in 
a  D.C.  plasma,  with  equally  successful  results. 

The  usefulness  of  the  formulae  which  have  been  presented 
in  this  chapter  is  that  they  may  be  used  to  predict  much 
of  the  observed  spatial  behaviour  of  many  atom  lines  in  the 
ICP .  It  is  not  universally  useful,  however,  since  it  has 
been  observed  that  above  16-18  mm  the  spatial  profiles  for 
many  lines  peak  regardless  of  the  value  for  the  normal 
temperature.  In  fact,  the  position  occupied  by  Cd  on  the 
plot  in  Figure  29  represents  the  point  at  which  spatial 
position  no  longer  follows  normal  temperature.  This  sug¬ 
gests  that  above  this  region  the  spatial  position  of  the  peak 
is  independent  of  analyte  excitation  and  ionization  charac¬ 
teristics.  In  addition,  it  would  be  useful  to  know  how 
these  analyte  parameters  affect  the  spatial  profiles  of  ion 


. 
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lines  which  are  not  "covered"  by  the  normal  temperature 
calculation.  These  points  are  addressed  in  the  next  section. 

Hard  Lines 

The  spatial  profiles  of  Znl  (213.8  nm)  and  Znll  (202.5 
nm)  provided  in  Figure  30  typify  the  type  of  behaviour 
observed  for  hard  type  lines  in  response  to  changes  in 
applied  power.  These  plots  are  similar  to  those  which  have 
been  discussed  in  Chapter  III.  The  numbers  refer  to  the 
applied  RF  power:  (1)  2.0  kw,  (2)  1.75  kw,  (3)  1.5  kw,  and 
(4)  1.25  kw.  The  position  of  maximum  intensity  for  these 

two  lines  is  seen  to  appear  in  approximately  the  same  place 
in  spite  of  the  fact  that  the  energy  requirements  for  the 
two  are  quite  different.  This  type  of  behaviour  is  also 
exhibited  for  Cdl  (228.8  nm)  and  Cdll  (214.4  nm)  as  is  shown 
in  Figure  31.  The  numbers  on  these  plots  are  analogous  to 
those  in  Figure  30.  Cdl  and  Cdll  clearly  behave  very  sim¬ 
ilarly  in  response  to  changes  in  RF  power.  In  fact  the 
maximum  in  emission  intensity  is  the  same  for  both  lines 
under  the  same  conditions  as  can  be  seen  in  Figure  32.  The 
spatial  profile  peaks  higher  in  this  figure  than  that  in 
Figure  29  because  the  aerosol  flow  rate  is  higher,  and  the 
power  in  this  case  is  2.0  kw.  This  close  coincidence  is 
puzzling  since  the  two  lines  have  widely  differing  energy 
requirements  to  reach  the  excited  state.  These  are  sum- 
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Emission  Intensify 


Cdl  228.8  nm 


Height  Above 


Cd II  226.5  nm 


Vertical  spatial  profiles  for  Cdl  228.8  nm  and 
Cdl I  226.5  nm.  RF  power;  (1)  2.0  kw,  (2)  1.75 
kw,  (3)  1.5  kw,  (4)  1.25  kw. 


Figure  31. 


Emission  Intensity 


Height  Above  Load  Coil  (mm) 

Figure  32.  Superimposed  vertical  spatial  profiles  of  Cdl 
and  Cdl I,  at  2.0  kw  RF  power. 


95 


marized  in  Table  VI.  Cdl  requires  5.41  ev  and  Cdll,  14.77 
ev.  This  information  suggests  that  there  may  be  a  wide  range 
of  energies  available  at  this  position  above  the  load  coil. 

This  position  also  represents  the  region  where  we  cannot 
use  the  "thermal"  normal  temperature  calculation  to  predict 
spatial  behaviour. 

The  relative  insensitivity  of  the  spatial  profile  to 
excitation  characteristics  is  further  demonstrated  in  Figure 
33.  The  spatial  profiles  of  five  Ba  ion  lines  with  ionization- 
excitation  sums  between  7.72  ev  and  11.21  ev  are  plotted. 

Table  VII  summarizes  the  energies  of  the  various  lines 
involved.  The  experimental  conditions  for  this  set  of 
profiles  are  1.5  kw  RF  power,  0.9  1pm  aerosol  flow  rate,  15 
lpm  plasma  gas  flow,  and  less  than  1  1pm  auxiliary  gas  flow. 
These  plots  confirm  that  relative  independence  of  spatial 
position  on  energy  requirements  of  the  different  lines, 
since  the  493.4  Ball  line  peaks  at  the  same  position  above 
the  load  coil  as  the  233.5  Ball  line.  Ca  ion  lines  are  seen 
to  behave  in  a  similar  manner  as  can  be  seen  in  Figure  34. 

In  this  case  the  ionization-excitation  sums  are  (a)  9.26  ev, 

(b)  12.6  ev ,  and  (c)  13.2  ev .  The  three  lines  peak  at  the 

same  position  above  the  load  coil. 

Discussion 


The  examples  discussed  in  the  previous  sections  lead  one 
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BaH 


Figure  33.  Vertical  spatial  profiles  for  Ball  lines;  (a) 

455.4  nm,  (b)  493.4  nn,  (c)  413.1  nn ,  (d)  489.2 

nm,  and  (e)  233.5  nm. 
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Height  Above  Load  Coil  (mm) 


Figure  3^.  Vertical  spatial  profiles  for  Call  lines;  (a) 
39.3  nm,  (b)  373.6  nm,  and  (c)  315.9  nm. 
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TABLE  VI 

WAVELENGTHS,  EXCITATION  POTENTIALS,  AND  SUM  OF  EX¬ 
CITATION  POTENTIAL  AND  IONIZATION  POTENTIAL  FOR  Cd 


Species  Line  (ran) 
Cdl  228.8 

Cdll  214.4 


EP  (ev) EP  +  IP  (ev) 

5.41  - 

5.78  14.77 


. 


TABLE  VII 


WAVELENGTHS,  EXCITATION  POTENTIALS  AND 

EXCITATION  POTENTIALS  AND 

IONIZATION  POTENTIALS  FOR  Ball 

SUM  OF 

EP  +  IP 

Line  (nm) 

EP  (ev) 

(ev) 

493.4 

2.51 

7.72 

455.4 

2.72 

7.93 

489.2 

5.69 

10.9 

413.1 

5.69 

10.9 

233.5 

6.00 

11.21 

100 


to  believe  that  the  aerosol  channel  may  be  divided  into  two 
regions.  The  first  is  a  "thermal"  region  where  soft  line 
behaviour  is  observed,  and  spatial  structure  may  be  predicted 
on  the  basis  of  a  normal  temperature  calculation.  This  is 
predominantly  an  atom  line  emission  region  and  extends  from 
just  above  the  load  coil  to  about  15  mm  above  the  load  coil 
at  1.5  kw.  In  this  region  it  is  probable  that  the  soft  line 
species  are  excited  predominantly  by  collisions  of  the  first 
kind  between  analyte  atoms  and  electrons.  The  effect  of  an 
increase  in  power  is  to  shift  this  zone  toward  the  load  coil 
and  to  increase  the  energy  of  the  electrons  slightly  leading 
to  a  downward  spatial  profile  shift  and  an  increase  in  the 
overall  emission  intensity.  The  temperature  in  this  region 
appears  to  increase  as  one  moves  up  the  central  channel. 

This  is  in  all  likelihood  the  result  of  the  diffusion  of  hot 
argon  gas  from  the  annular  region  toward  and  into  the  channel 
region.  "Cold"  argon  gas  injected  into  the  bottom  of  the 
plasma  moves  up  the  aerosol  channel  becoming  hotter  as  it 
mixes  with  "hot"  annular  argon  gas.  An  increase  in  aerosol 
flow  rate  delays  this  process  leading  to  the  observed  upward 
shift  in  the  spatial  emission  profile  of  the  analyte. 

A  "non-thermal "  region  where  hard  lines  peak,  and  spatial 
behaviour  is  insensitive  to  analyte  excitation  parameters , 
characterizes  the  upper  region  of  the  analyte  channel.  By 
the  time  the  analyte  reaches  this  region,  which  extends  from 
about  15  mm  to  25  mm  above  the  load  coil,  annular  and  aerosol 


. 
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gas  should  be  thoroughly  mixed.  The  species  responsible  for 
excitation  in  this  region  appears  capable  of  providing  the 
high  energy  required  to  ionize  and  excite  species  such  as 
Cd,  Zn,  Ba,  and  Ca.  Since  this  region  is  well  away  from  the 
load  coils  it  is  unlikely  that  electrons  would  possess  suf¬ 
ficient  energy  to  excite  these  species.  Ion  lines  emitting 
in  this  region  are  often  the  most  sensitive  lines  for 
analytical  determinations  in  the  ICP .  This  ion  line  advantag 
has  been  attributed  to  an  excitation  process  involving  argon 
metastable  species  as  mentioned  in  Chapter  I.  It  has  been 
suggested  that  these  metastable  species  may  have  lifetimes 
sufficiently  long  to  exist  in  this  region  [101] . 

In  the  case  of  emission  from  Cd  ion  the  value  of  14.77 
ev  for  ionization  plus  excitation  precludes  the  direct 
Penning  reaction, 

met a  +*  — 

Ar  tcl  +  Cd  +  Ar  +  Cd  +  e  (34) 

since  the  metastable  species  have  energies  of  11.55  ev  and 

11.71  ev.  However  the  Penning  process, 

Armeta  +  cd  Ar  +  Cd+  +  e~  (35) 

followed  by  a  collision  of  the  first  kind  with  an  electron 
of  the  proper  energy, 

Cd+  +  e  ->  Cd+  +  e"  (36) 

is  possible.  The  close  correspondence  of  the  spatial  pro¬ 
files  for  Cdl  and  Cdll  also  suggests  that  the  mode  of  ex- 
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citation  of  these  two  species  is  related.  Cdl  may  be  sig¬ 
nificantly  populated  by  the  ion-electron  recombination 
reaction  [5] , 

Cd+  +  e”  +  Cd*  +  hv  (37) 

such  that  Penning  produced  Cd+  may  act  as  a  reservoir  for 
both  Cdl  and  Cdll  emission. 

In  the  case  of  Ba  the  energies  involved  make  the  ex¬ 
citation  reaction 

met a  +*  — 

Ba  +  Ar  +  Ba  +  e  +  Ar  (38) 

possible  for  all  of  the  energy  levels  which  were  observed. 
This  could  be  used  to  explain  the  close  correspondence  in 
spatial  behaviour  for  these  lines. 

The  division  between  these  two  regions  is  not  a  distinct 
one.  The  regions  probably  overlap  significantly.  There  is 
also  the  possibility  that  ionization  in  the  lower  channel 
region  is  partially  the  result  of  a  Penning  ionization  pro¬ 
cess  such  as  that  outlined  in  Equation  35.  If  this  is  the 
case  it  is  difficult  to  assess  the  degree  to  which  it  will 
affect  the  spatial  structure  in  the  lower  channel.  The 
density  of  argon  metastable  species  is  undoubtably  closely 
related  to  the  temperature.  Thus,  the  correlation  of  the 
spatial  structure  of  soft  lines  with  a  thermal  ionization- 
excitation  treatment  only/  might  be  misleading,  even  though 
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it  is  quite  successful.  Further  experiments  are  required 
to  unravel  the  problem.  A  detailed  spatial  map  of  argon 
metastable  densities  would  be  very  useful.  The  effects  of 
a  metastable  quencher  should  also  be  extensively  investigated. 


CHAPTER  V 


IMPLEMENTING  ABEL  INVERSIONS  USING  A 
PHOTODIODE  ARRAY  MEASUREMENT  SYSTEM 


Introduction 


Up  to  this  point  the  discussion  of  analyte  spatial 
behaviour  has  been  limited  to  observations  on  a  vertical 
spatial  zone  up  the  centre  of  the  plasma.  Using  2:1  imaging 
this  means  that  there  is  only  a  50  ym  wide  spatial  obser¬ 
vation  window.  The  analyte  channel  is  much  wider  than  this, 
however,  and  it  would  be  useful  to  measure  emission  intensity 
across  the  entire  width  of  the  channel.  In  addition, 
since  the  ICP  is  viewed  from  the  side,  the  emission  in¬ 
tensity  observed  is  the  result  of  contributions  from  all 
depths  within  the  source,  and  is  not  solely  limited  to  axial 
information.  From  a  fundamental  point  of  view  it  would  be 
useful  to  be  able  to  measure  the  emission  from  the  ICP  with 
radial  resolution,  that  is,  such  that  one  can  separately 
observe  emission  from  varying  depths  within  the  source 
without  contributions  from  other  depths  within  the  optical 
path. 

The  most  common  way  of  solving  this  problem  has  been 
to  acquire  a  number  of  data  points  at  right  angles  to  the 
source  axis,  or  "laterally".  These  points  may  subsequently 
be  operated  on  using  a  mathematical  process  known  as  an 
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"Abel  Inversion"  [78],  This  process  converts  measured, 
spatially  integrated  lateral  information,  and  converts  it 
to  spatially  resolved  radial  information.  The  measurement 
of  radial  emission  intensities  from  the  ICP  using  this 
Procedure  is  not  new.  Kornblum  and  deGalan  designed  a 
measurement  system  to  measure  radial  spatial  distributions 
of  emission  intensity  in  the  ICP  [3] ,  and  used  it  to  cal¬ 
culate  temperature,  and  electron  densities  [35], as  well  as 
to  study  matrix  effects  [76].  Kalnicky,  et  al.  [36],  Jarosz, 
et  al.[37],  and  Alder,  et  al.  [38]  have  all  used  this 
technique  to  measure  emission  intensities  for  the  purpose 
of  calculating  radially  resolved  temperatures  and  electron 
densities . 

The  generation  of  a  useful  radial  intensity  distri¬ 
bution  from  a  lateral  intensity  distribution  imposes 
several  requirements  on  the  quality  of  the  measured  lateral 
data.  These  requirements  have  been  discussed  previously 
[3] .  One  of  these  requirements  is  that  the  lateral  data  be 
as  detailed  as  possible;  the  greater  number  of  data  points, 
the  more  accurate  will  be  the  generated  radial  distribution. 
The  acquisition  of  detailed  spatial  emission  data  from  a 
spectrochemical  source  such  as  the  ICP  requires  a  very  large 
number  of  measurements  if  complete  information  is  to  be 
obtained.  Often  the  tedium  of  such  measurements  limits 
the  range  of  conditions  and  number  of  analyte  species  which 
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can  be  studied  in  a  reasonable  amount  of  time.  This  problem 
has  been  circumvented  in  this  study  by  using  the  photodiode 
array-based  measurement  system.  The  purposes  of  this  chapter 
are  to  detail  modifications  to  the  previously  described 
profiling  system  such  that  it  can  be  used  to  acquire  lateral 
emission  intensities ,  and  to  outline  how  this  data  is  con— 
verted  into  radial  emission  intensity  information.  Portions 
of  this  chapter  have  been  previously  published  [83] . 

Experimental 

The  monochromator/diode  array  measurement  system  has 
been  previously  described  in  Chapter  II.  In  order  to  obtain 
complete  detailed  maps  of  analyte  emission  throughout  the 
ICP  this  system  has  been  modified  so  that  a  series  of  hor¬ 
izontal  emission  profiles  can  be  acquired.  In  order  to  obtain 
horizontal  emission  profiles  the  monochromator  was  simply 
rotated  90°  so  that  the  entrance  slit  was  horizontal  rather 
than  vertical.  This  was  implemented  by  constructing  a 
monochromator  mount  from  aluminium  schematically  outlined 
in  Figure  35.  This  mount  holds  the  monochromator  securely 
in  place  and  mates  with  the  optical  rail  bed.  In  order  to 
measure  horizontal  slices  of  the  plasma  at  different 
vertical  heights  the  monochromator  and  lens  were  mounted  on 
a  3  point  screw  adjustable  kinematic  mount  which  could  be 
raised  or  lowered  with  a  precision  of  0.01  mm  measured  using 
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ALIGNMENT 

LASER  MONOCHROMETER 


Figure  35.  Horizontal  profiling  system. 
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a  precision  dial  guage  (Mitutoyo  No.  3052,  Japan).  The 
plasma  is  imaged  onto  the  entrance  slit  of  the  monochromator 
using  a  single  spherical  lens  with  10  cm  focal  length  and  a 
magnification  of  1/2.  Since  the  individual  diodes  are  on 
25.4  mm  centres  and  25.4  mm  wide  the  array  samples  the 
emission  from  the  plasma  in  256  separate  50.8  ym  steps. 

This  measurement  capability  of  the  diode  array-based  system 
results  in  high  quality  input  data  for  the  Abel  inversion. 
The  acquisition  of  256  separate  lateral  intensities  at  such 
fine  spacing  by  horizontal  translation  of  a  conventional 
measurement  system  would  indeed  be  tedious.  Using  this 
system  and  the  Abel  inversion  technique  it  is  possible  to 
obtain  horizontal  profiles  from  0  to  30  mm  above  the  load 
coil,  and  thus  build  a  complete  3-dimensional  picture  of 
emission  in  the  plasma. 

Abel  Inversion  of  Lateral  Data 

The  process  of  calculating  a  set  of  radial  data  points 
from  a  corresponding  set  of  lateral  data  points  is  known  as 
the  Abel  inversion.  The  application  of  this  technique  to 
spectroscopic  data  is  well  known  [78]  and  good  discussions 
of  the  technique  as  applied  to  spark  sources  [1,2,81],  and 
ICP  [7,35-37]  are  available  and  will  not  be  repeated  in 
this  thesis.  A  short  discussion,  however,  is  warranted. 

The  relationship  between  the  measured  lateral  in¬ 
tensities  ln(x)  as  measured  with  the  diode  array  and  the 
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desired  radial  intensities  i  (r)  calculated  by  the  Abel 
inversion  is  shown  schematically  in  Figure  36.  The  diode 
number  index  is  n  and  the  radial  index  is  k.  We  use  a  256 
element  diode  array  centered  on  the  plasma;  thus  the  actual 
maximum  value  of  N  is  128. 

Following  the  method  of  Nestor  and  Olsen  [79]  the 
radial  intensity  i^(r)  may  be  expressed  in  terms  of  the 
measured  lateral  intensity  In(x)  by  the  summation  [78]: 


(39) 


where  k  and  n  are  integer  position  indicies  for  the  radial 
and  lateral  intensities,  and  a  is  the  width  of  the  slice 


sampled  by  a  single  photodiode.  The  coefficient  B 


k,n 


can 


be  evaluated  by  using  the  expressions 


B 


k,n 


-A. 


k ,  k 


for  n  =  k 


(40) 


and 


(41) 


where 


A. 


k,  n 


=  [n2- (k-1 ) 2 ] 1/2  -  [ (n-1 ) 2  -  (k-l)2]1/2 

2n-l 


(42) 


To  test  the  accuracy  of  the  inversion  calculation  a 


set  of  radial  test  data  was  generated  using  the  polynomial 
expression  shown  in  Equations  (43)  and  (44)  [78]. 
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i  (r )  =  0.75  +  12r2  -  32r3  O^rlO.25 
i (r )  =  16/27 (l+6r-15r2+8r3)  0.25<r-1.0 


(43) 

(44) 


The  corresponding  set  of  lateral  intensities  may  be  obtained 
by  direct  integration  of  the  expression 


I  (x)  =2  / 


1  i  (r )  rdr 

*  (r2-xV/2 


(45) 


The  set  of  lateral  data  points  thus  generated  were 
used  to  calculate  a  radial  array  using  the  Abel  inversion 
program.  This  program  runs  on  the  PDP-8  and  is  listed  in 
Appendix  C.  The  procedure  utilized  is  to  read  lateral  data 
stored  on  magnetic  tape  into  the  computer,  perform  the  Abel 
inversion,  and  store  the  resultant  radial  data  back  on 
magnetic  tape.  The  result  of  the  test  is  shown  in  Figure 
37.  Curve  (a)  is  the  radial  test  data  generated  using 
Equations  (43)  and  (44).  Curve  (b)  is  the  lateral  test 
data  which  is  calculated  from  the  integral  in  Equation  (45), 
and  curve  (c)  is  the  radial  profile  which  is  calculated 
using  our  Abel  inversion  program.  Curve  (c)  shows  a  small 
amount  of  scatter  near  the  centre  of  the  trace  (r=0)  which 
is  due  to  the  recursive  nature  of  the  calculation  but 
generally  it  compares  very  well  with  the  theoretical  profile 
of  curve  (a) . 

An  example  calculation  using  spatial  data  acquired  from 
an  ICP  is  shown  in  Figure  38  in  order  to  illustrate  the 
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sequence  of  processing  from  raw  horizontal  data  to  cal¬ 
culated  radial  data.  The  spectral  signal  monitored  was  the 
Mg I I  279.5  nm  line  measured  at  a  height  of  12  mm  above  the 
load  coil.  The  horizontal  spatial  emission  profile  as 
measured  with  the  diode  array  measurement  system  is  shown 
in  Figure  38a.  It  takes  only  about  30  sec  to  acquire  this 
signal.  An  average  one-sided  profile  is  then  calculated 
by  summing  the  intensities  of  diodes  equidistant  from  the 
center  and  dividing  by  two.  We  found  this  step  to  be 
necessary  because  the  Abel  inversion  is  very  sensitive  to 
slight  asymmetries  and  a  discontinuity  results  at  the 
center  of  the  radial  data  for  slightly  asymmetric  lateral 
data.  The  averaged  one-sided  lateral  profile  is  shown  in 
Figure  38b.  The  radial  profile  calculated  from  this  lateral 
profile  is  shown  in  Figure  38c.  Note  the  difference  in  shape 
of  the  lateral  and  radial  emission  profiles.  The  lateral 
profile  of  Figure  38b  implies  that  emission  peaks  at  the 
center  of  the  analyte  channel  whereas,  the  radial  profile 
in  Figure  3_8c  indicates  that  the  emission  in  fact  dips 
slightly  at  the  center  of  the  channel  peaking  just  off  axis. 

Results  and  Discussion 

This  system  has  been  used  to  study  analyte  emission 
from  the  ICP,  by  constructing  3-dimensional  "maps"  of  radial 
emission  intensity.  This  is  carried  out  by  acquiring 
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Figure  38. 
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lateral  profiles  at  different  heights  above  the  load  coil 
and  storing  them  on  magnetic  tape.  These  profiles  may  be 
subsequently  inverted  and  plotted  in  the  manner  shown  in 
Figures  39,  40,  41  and  42.  In  these  figures  emission 
intensity  from  the  observed  analyte  is  plotted  on  the  ver¬ 
tical  axis  as  a  function  of  radial  coordinate.  The  radial 
coordinate  is  the  distance  in  mm  out  from  the  centre  of  the 
aerosol  channel.  A  series  of  these  are  plotted  at  dif¬ 
ferent  heights  above  the  load  coil  resulting  in  the  3-di¬ 
mensional  plots  provided. 

The  complete  radial  emission  profiles  for  Cal  (422.7 
nm) ,  Call  (393.3  nm) ,  Mgl  (285.3  nm) ,  and  Mgll  (279.5  nm) 

are  provided  in  Figures  39,  40,  41,  and  42  respectively. 

The  ICP  operating  conditions  are  RF  power  1.5  kw,  plasma 
gas  flow  rate  15  1pm,  coolant  gas  flow  rate,  less  than  1 
lpm,  and  aerosol  flow  rate,  0.8  1pm  for  the  Mg  profiles,  and 
0.9  lpm  for  the  Ca  profiles.  These  profiles  show  several 
interesting  spatial  features  which  seem  characteristic  of 
analyte  emission  from  the  analyte  channel.  Low  in  the 
plasma  all  four  lines  have  similar  spatial  structure. 
Emission  peaks  up  off  axis  with  the  peak  emission  cor¬ 
responding  roughly  to  the  edge  of  the  analyte  channel. 

This  is  good  evidence  that  the  excitation  lower  in  the 
channel  is  related  to  contact  with  hot  argon  gas  orig¬ 
inating  from  the  annular  region  of  the  plasma.  It  is  also 
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Figure  39.  Complete  radial  emission  map  of  the  ICP  discharge 
for  Cal  422.7  nn. 
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Figure  40.  Complete  radial  emission  map  of  the  ICP  discharge 
for  Call  393.3  nm. 
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Figure  41.  Complete  radial  emission  map  of  the  ICP  discharge 
for  Mg I  285.3  nm. 
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Figure  42.  Complete  radial  emission  map  of  the  ICP  discharge 
for  Mg I I  279.5  nm. 
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verification  of  an  assumption  made  by  Human  and  Scott  [43] 
in  explaining  the  shapes  of  spectral  lines  in  the  ICP. 

They  noted: 

Low  in  the  plasma  (at  10  mm  height)  self —reversal 
of  spectral  lines  does  not  occur  for  high  concen¬ 
trations  of  analyte,  but  the  lines  are  considerably 
broadened.  The  unusual  shape  of  the  lines  can  be 
explained  by  assuming  that  the  core  of  the  plasma 
is  a  predominantly  absorbing  region,  while  the 
emission  originates  from  the  annular  region  around 
the  core . " 

It  can  be  seen  that  this  is  indeed  the  case. 

With  an  increase  in  observation  height  the  position  of 
the  off  axis  peak  moves  closer  to  the  centre  of  the  channel 
until  at  a  height  of  about  12-14  mm  above  the  load  coil  the 
emission  peaks  in  the  center  of  the  channel.  This  can  be 
interpreted  as  an  increase  in  the  amount  of  mixing  of  cold 
channel  gas  and  hot  annular  gas  with  an  increase  in  height, 
such  that  complete  mixing  occurs  at  this  height.  Emission 
continues  to  peak  at  higher  observation  heights.  The 
emission  from  ion  lines  is  characterized  by  a  somewhat 
broader  analyte  emission  plume  at  heights  between  14  and 
22  mm  above  the  load  coil,  when  compared  to  emission  from 
atom  lines.  This  might  be  the  result  of  ambipolar  diffusion 
of  ion  species.  This  phenomenon  will  be  discussed  in 
detail  in  Chapter  VI  of  this  thesis.  Also,  in  agreement 
with  the  data  presented  in  previous  chapters  the  neutral 
atom  lines  are  observed  to  peak  lower  in  the  analyte  channel. 
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In  general  the  analyte  channel  is  characterized  by  a 
structure  which  starts  out  low  in  the  plasma  somewhat  like 
a  hollow  bullet  and,  then  higher  in  the  plasma,  in  the 
region  normally  observed  for  analytical  purposes,  becomes 
a  solid  emitting  plume.  This  structure  suggests  that  the 
temperature  of  the  channel  is  hotter  near  the  edges  low 
in  the  channel  and  hotter  in  the  middle  at  higher  observation 
heights.  To  test  this,  radially  resolved  excitation  temp¬ 
eratures  have  been  measured. 

Spatial  Distribution  of  Excitation  Temperatures 

Excitation  temperature  has  been  measured  using  the 
two  line  ratio  method  [55,56,102-105] .  The  method  involves 
the  measurement  of  emission  intensity  for  two  lines  of  the 
same  species,  and  calculating  the  temperature  from  the  ratio 
of  these  two  lines,  according  to  the  formula: 

g  A  g  A 

ln  1rD'/1n  =  ln  I  ]  +  (E  -E  )/kT  (46) 

p  g  xp  Ag  g  P 

where 

=  emission  intensity  from  excited  state  p 
!g  =  emission  intensity  from  excited  state  g 
9  A  =  transition  probability  of  p 

ir  ir 

g  A  =  transition  probability  of  g 
y  y 

A  =  wavelength  of  p  transition 
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X  =  wavelength  of  g  transition 
y 

E  =  excitation  energy  of  state  g 
y 

Ep  =  excitation  energy  of  state  p 
k  =  Boltzmann's  constant 
T  =  excitation  temperature 

The  two  lines  used  for  this  measurement  were  Pbl 
280.1  nm  and  Pbl  283.3  nm.  Transition  probabilities  were 
taken  from  Moore  [106],  The  wavelength,  transition  pro¬ 
babilities  and  excitation  potentials  for  these  two  lines 
are  listed  in  Table  VIII.  Radial  profiles  for  these  two 
lines  were  acquired  at  various  heights  above  the  load  coil 
and  the  ratio  calculated  to  provide  the  temperature  profiles 
given  in  Figure  43.  The  conditions  used  for  this  study  were: 
1.5  kw  RF  power,  15  1pm  plasma  gas  flow  rate,  less  than  1 
lpm  auxiliary  flow  rate,  and  0.8  1pm  aerosol  flow  rate.  In 
Figure  43  temperature  is  plotted  as  a  function  of  radial 
distance  from  the  centre  of  the  plasma  from  0  to  3  mm,  at 
the  different  heights  noted  on  the  figure.  Temperatures 
ranged  from  3500  K  to  4900  K  depending  on  spatial  position. 
Errors  in  these  values  have  not  been  included,  so  they  should 
be  taken  as  approximate  only;  however,  the  relative  mag¬ 
nitudes  may  be  compared.  At  low  observation  heights  between 
4  mm  and  8  mm  there  may  be  inaccuracies  in  the  Abel  inversion 
as  the  result  of  off  axis  peaks,  so  these  values  near  the 
analyte  channel  may  not  be  representative. 
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TABLE  VIII 


WAVELENGTHS,  TRANSITION  PROBABILITIES,  AND  EX¬ 
CITATION  ENERGIES  FOR  Pbl  USED 
TEMPERATURE  DETERMINATION 


A  (A° ) 

gA(sec  1) 

E  (ev) 

Pbl 

P 

2801.99 

43  x  io8 

5.7439 

Pbl 

g 

2833.06 

00 

o 

r— 1 

X 

00 

• 

1 — 1 

4 . 37492 

5000 


124 


(»)  ejr^DJaduuaj^ 


0) 

tn 

P 

fd 

A 

U 

to 

•H 

Td 

P4 

u 

H 

CL) 

£ 

-P 

P 

O 

M-l 

CO 

0) 

I — I 
•H 
4-f 
O 
P 
& 

CD 

P 

P 

-P 

fd 

p 

d) 

cu 

e 

<u 

■p 


fd 

•H 

Td 

(d 

& 


n 

<D 

P 

P 

tn 

■P 


125 


The  results  confirm,  however,  that  low  in  the  analyte 
channel,  temperature  is  higher  near  the  edges 
as  can  be  seen  in  the  plots  for  4  to  12  mm  above  the  load 
c°il-  At  higher  observation  heights  this  trend  is  reversed 
and  temperatures  peak  near  the  centre  of  the  analyte  channel 
dropping  with  increasing  radial  distance,  as  can  be  seen  in 
the  plots  for  18  to  24  mm  above  the  load  coil. 

Conclusion 


Clearly  the  ability  to  make  the  kind  of  measurements 
described  in  this  chapter  has  added  greatly  in  the  spatial 
characterization  of  analyte  emission  in  the  ICP .  The  ease 
and  speed  of  acquisition  of  spatial  information  from  the 
plasma  can  be  tremendously  increased  using  the  photodiode 
array .  High  resolution  lateral  data  are  simultaneously 
acquired  in  an  ideal  format  for  Abel  inversion.  Thus  the 
sometimes  tedious  task  of  collecting  a  large  number  of 
lateral  points  to  ensure  accuracy  of  the  inversion  is 
eliminated.  The  ease  of  acquisition  means  that  complete 
3-dimensional  pictures  of  emission  may  be  quickly  built  up 
which  can  provide  detailed,  quantitative  information  about 
analyte  emission.  These  pictures  have  proved  invaluable 
in  unraveling  the  complex  problem  of  the  effect  of  easily 
ionizable  elements  on  analyte  emission  intensity.  This  is 
discussed  in  the  next  chapter. 
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CHAPTER  VI 


INTERFERENCE  FROM  EASILY  IONIZABLE  ELEMENT  MATRICES 
Introduction 

A  literature  survey  of  the  effects  of  an  excess  of  an 
easily  ionizable  element  (EIE)  such  as  Na,  K,  or  Cs  as  a 
concomitant  on  analyte  emission  intensity  in  ICP 1 s  is  a  study 
in  confusion.  Virtually  every  ICP  spectroscopist  has  studied 
the  effect  of  this  type  of  matrix,  and  the  literature  abounds 
with  tables  and  graphs  on  the  magnitude  and  characteristics 
of  the  effect.  The  driving  force  for  these  studies  is  the 
fact  that  EIE's,  especially  Na  and  K,  turn  up  in  a  variety  of 
natural  matrices  for  which  emission  spectroscopy  is  used  as 
an  analytical  probe.  These  include  animal  and  plant  materials, 
blood,  seawater,  and  geological  samples. 

The  sources  of  much  of  this  confusion  are  twofold. 

First,  much  of  the  data  which  has  appeared  has  not  been  con¬ 
sistent,  some  investigators  observing  analyte  emission  "en¬ 
hancements",  some  observing  "depressions",  and  some  "no 
effect".  Secondly,  there  is  still  no  satisfactory  explanation 
for  the  observed  effects  at  a  fundamental  level.  To  clarify 
some  of  these  observations  and  offer  rational  arguments 
on  the  nature  of  the  effect  a  spatial  study  of  the  effect 
has  been  undertaken  and  forms  the  basis  of  this  chapter.  Some 
preliminary  results  have  been  published  [77] . 
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Survey  of  the  EIE  Effect  in  Plasmas 

The  effect  of  EIE's  in  analytical  flame  spectroscopy  has 
been  well  characterized.  Enhancements  in  emission  intensity 
from  atom  lines  is  observed  when  signals  with  and  without  the 
presence  of  an  EIE  are  compared  [107,108],  These  observations 
have  been  rationalized  on  the  basis  of  a  shift  in  the  ionization 
equilibrium  for  analyte  atom  (x) ,  ion  (x  )  and  electrons  (e  ) . 

x  t  x+  +  e~ 

The  addition  of  the  EIE  has  the  effect  of  increasing  the 
electron  density  in  flames  such  that  the  equilibrium  shifts 
toward  the  neutral  atom  species  causing  enhancement  of  neutral 
atom  line  emission  [107,109],  The  magnitude  of  the  effect 
has  been  shown  to  be  related  to  the  ionization  potential  (IP) 
of  the  added  EIE,  being  larger  for  those  of  lower  IP  and 
reducing  in  magnitude  as  the  IP  of  the  added  element  increases 
[110] .  Thus,  in  terms  of  the  effect,  the  sequence, 

Cs  >  Rb  >  K  >  Na  >  Sr  >  Li  >  Mg 

is  observed.  Quite  a  different  situation  is  observed  in  an 
arc  discharge.  The  addition  of  an  EIE  to  the  arc  increases 
the  electrical  conductivity  of  the  arc  but  leaves  the  thermal 
conductivity  essentially  unchanged  [56].  The  arc  responds 
to  this  by  a  lowering  in  temperature  such  that  the  power 
supplied  to  the  arc  is  once  again  balanced  by  energy  loss  due 
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to  thermal  conduction.  Thus  analyte  emission  intensities  are 
observed  to  decrease  due  to  the  presence  of  EIE's. 

The  situation  is  not  quite  so  well  defined  in  the  ICP , 
partially  because  of  the  variability  in  plasma  operating 
conditions  between  investigators,  and  partially  because  the 
influence  of  EIE's  has  been  observed  to  be  complex  in  nature. 
Hoare  and  Mostyn  [28]  observed  that  1000  yg/ml  Na  did  not 
affect  the  emission  intensity  of  B  and  that  10,000  yg/ml  pro¬ 
duced  slight  enhancements  for  atomic  lines  and  slight  de¬ 
pressions  for  ionic  lines.  Similarly  Kirkbright  et  al. 

[39,40]  reported  that  the  emission  intensities  of  S,  P,  As, 

I,  Se,  and  Hg  were  not  affected  by  the  addition  of  K  or  Na 
in  50  fold  weight  excess  representing  0 . 05  to  0.5%  by  weight 
salt.  In  the  first  detailed  study,  Larson  et  al.  [48]  in¬ 
vestigated  the  effects  of  up  to  a  600  to  1  molar  ratio  excess 
of  Na  on  the  emission  intensity  of  Ca,  Cd,  and  Cr.  The  ion 
lines  were  all  observed  to  be  depressed  by  the  addition  of 
Na  at  heights  above  15  mm,  the  magnitude  of  the  depression 
increasing  with  increasing  Na  concentration.  The  magnitude  of 
the  depression  was  for  the  most  part  independent  of  the  ob¬ 
servation  height  between  15  to  25  mm  above  the  load  coil,  at 
low  aerosol  flow  rates  (1.0  1pm) .  An  increase  in  RF  power 
from  1025  watts  to  1250  watts  likewise  had  little  effect  on 
the  magnitude  of  the  depression.  At  1025  watts  and  1.3  1pm 
flow  rate,  however,  ion  emission  was  enhanced  between  15  and 
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17  mm  above  the  load  coil  and  depressed  thereafter.  This 
was  an  indication  that  there  may  be  complex  spatial  features 
associated  with  the  effect.  The  emission  behaviour  of  the 
atom  lines  was  more  complex.  Ca  atom  emission  was  in  general 
enhanced,  the  magnitude  increasing  with  increasing  molar 
ratios  of  Na  to  Ca.  The  enhancement  increased  with  increasing 
height  above  the  load  coil.  An  increase  in  power  increased 
the  magnitude  of  the  effect  as  did  an  increase  in  the  aerosol 
flow  rate.  Cr  atom  emission  on  the  other  hand  was  depressed 
at  lower  observation  heights  and  enhanced  at  higher  observation 
heights.  An  increase  in  both  power  and  aerosol  flow  rate 
increased  the  effect.  Cd  atom  emission  was  depressed  at  all 
heights,  insensitive  to  power  changes,  but  changed  to  en¬ 
hancements  at  higher  aerosol  flow  rates.  They  concluded  that 
there  was  some  evidence  to  suggest  an  ionization  supression 
mechanism,  but  they  also  commented  that  other  factors  must 
play  a  role  since  some  of  the  data  conflicted  with  this 
interpretation . 

Abdallah  et  al.  [45]  carried  out  a  similar  study  in 
which  the  emission  intensities  of  PI,  Cal,  Call,  Cdl,  Cdll, 

CrI,  and  CrII  were  observed  with  additions  of  Na  up  to  0.1%, 
at  observation  heights  of  5,  15,  and  25  mm  above  the  load 
coil.  Cal  was  enhanced  and  Call  depressed  only  slightly  in 
agreement  with  the  observations  of  Larson  et  al.  Cdl  and  Cdll 
were  both  observed  to  be  depressed  slightly  at  all  heights. 
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CrI  was  relatively  insensitive  to  Na  addition,  whereas  CrII 
showed  slight  enhancement  and  then  a  larger  depression  as 
Na  concentration  was  increased.  The  effect  became  less 
pronounced  at  higher  observation  heights.  Emission  from 
atomic  phosphorus  was  enhanced  significantly  at  5  mm  above 
the  load  coil,  and  slightly  at  25  mm  above  the  load  coil.  The 
variation  of  aerosol  flow  rate  was  found  to  cause  complex 
changes  in  the  nature  of  the  effect.  An  increase  in  applied 
power  was  found  to  reduce  the  effect  on  Cal  and  Call  until 
at  2  Kw  no  effect  could  be  observed.  By  measuring  the  ratio 
of  the  Mgl  (285.2  nm)  and  Mgll  (279.6  nm)  emission  lines,  and 
the  ratio  of  two  Till  lines  (322.2  and  322.4  nm)  the  effects 
of  added  Na  on  electron  density  and  temperature  were  found  to 
be  negligible.  It  was  also  observed  that  Na  exerted  the  same 
influence  on  emission  and  absorption  signals  which  led  them 
to  believe  that  the  effect  of  Na  was  largely  due  to  a 
volatalization  interference.  Veillon  and  Margoshes  [113]  also 
performed  an  emission/absorption  measurement  on  Cal  and  Call 
and  found  no  change  in  the  absorbance  signal  for  Cal  in  the 
presence  of  Na  even  though  the  emission  signal  increased  by 
a  factor  of  2.7.  The  plasma  studied,  however,  was  a  tear  drop 
shape,  rather  than  an  annular  plasma.  Boumans  and  DeBoer 
[46]  investigated  the  effects  of  K  on  analyte  emission  in¬ 
tensity  as  part  of  the  optimization  procedure  for  simultan¬ 
eous  multielement  analysis  using  ICP .  For  Lil  a  fivefold  en- 
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enhancement  dropped  to  a  20%  depression  as  the  aerosol  flow 
rate  was  decreased  from  1.7  1pm  to  1.1  1pm.  Similar  results 
were  observed  for  Mnl  and  Ball .  The  effect  decreased  sharply 
with  an  increase  in  power  from  0.4  Kw  to  1.0  Kw,  such  that 
at  1.0  Kw  the  effect  was  very  small.  Under  the  optimum  con¬ 
ditions  for  multielement  analysis,  0.7  Kw  power  input,  1.3 
lpm  aerosol  flow  rate  and  15  mm  observation  height, the  ef¬ 
fects  were  found  to  be  minimal,  +20%  for  Lil,  +  15%  for  Sri, 

-3%  for  All,  -2%  for  Mnl,  0%  for  Cdl ,  and  -6%  for  SrII  and 
Lall,  for  0.1%  K  as  interferent.  Similar  magnitudes  were 
found  for  Cs  as  an  interferent. 

In  a  subsequent  publication  [47]  using  ultrasonic  neb- 
ulization,  Boumans  and  DeBoer  separated  out  nebulizer/des¬ 
olvation  effects  and  plasma  effects.  In  general  enhancements 
of  3%  to  36%  were  observed  for  atom  line  emission,  the  excep¬ 
tion  being  Mn  where  a  depression  of  14%  was  observed.  Ion 
lines  were  depressed  from  2%  to  22%  for  all  the  elements 
studied.  Nebulizer/desolution  effects  were  found  to  be 
significant  in  only  3  of  the  9  elements  studied,  the  effect 
being  in  the  range  4%  to  7%  compared  with  21%  to  29%  for  the 
analogous  plasma  effect.  The  magnitude  of  the  effect  correl¬ 
ated  with  the  ionization  potentials  of  the  elements  Ba,  Sr, 

V,  Mn,  Fe,  decreasing  with  increasing  ionization  potential. 

The  results  for  Pd,  Zn,  and  Cd  did  not  follow  this  trend  how¬ 
ever  .  The  effect  of  K  on  Bal  was  found  to  be  related  to  molar  con- 
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centration  of  K  rather  than  the  molar  ratio  of  K  to  Ba .  They 
suggested  that  ionization  effects  are  partially  the  cause 
of  the  effect  of  K,  but  pointed  to  a  possible  volatilization 
effect  as  the  result  of  water  "included"  in  desolvated 
particles;  or  retarded  volatilization  due  to  the  matrix 
causing  "lateral  diffusion  effects". 

Watson  and  Steele  [111]  studied  the  effects  of  up  to 
3%  added  Na  on  the  emission  intensities  from  several  elements 
in  an  ICP  operating  at  5.5  kw,  2.0  1pm  aerosol  flow  rate  with 
nitrogen  as  the  coolant  gas.  Fell,  Call,  BI,  Sal,  All,  Cdll, 
and  CrII  all  showed  enhancements  at  10  mm  above  the  load  coil, 
and  depressions  at  15  mm  above  the  load  coil.  The  magnitude 
of  the  interference  was  found  to  increase  with  increasing  Na 
concentration,  except  for  BI  which  showed  a  sharp  initial 
enhancement  of  8%  which  decreased  to  about  2%  at  a  Na  con¬ 
centration  of  30  g/1.  The  magnitude  of  the  depressions  at 
15  mm  above  the  load  coil  were  much  larger  than  the  analogous 
enhancements  at  10  mm  above  the  load  coil.  For  the  most  part 
an  increase  in  power  was  found  to  decrease  the  matrix  effect. 
No  explanations  were  offered  for  the  observations. 

Clearly  the  exact  nature  of  the  effect  of  easily  ion- 
izable  elements  on  analyte  emission  is  complex.  The  picture 
is  further  clouded  by  the  variety  of  operating  conditions 
chosen  by  different  investigators.  Nevertheless,  we  would 
not  expect  the  physical  interactions  among  EIE ,  analyte  and 
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plasma  gas  to  be  substantially  different  from  system  to 
system  and  some  postulate  which  can  unify  the  observations 
is  required.  Boumans  and  DeBoer  [47]  noted  that,  "Fund¬ 
amental  studies  using  spatial  resolution  will  have  to  clear 
up  the  picture  of  plasma  interferences...." 

Kornblum  and  deGalan  [76]  undertook  the  first  de¬ 
tailed  spatial  study  of  the  effects  of  EIE's  on  analyte 
behaviour  in  the  ICP .  They  used  spatially  resolved  emis¬ 
sion  and  absorption  measurements  to  investigate  the  effect 
of  Cs  on  ground  state  population  and  emission  from  excited 
states  of  Ca  and  Mg  atom  and  ion  species.  They  made  the 
that  since  high  aerosol  flow  rates  have  been  obser¬ 
ved  to  magnify  the  effect,  they  would  use  this  condition 
such  that  the  effect  could  be  more  easily  studied.  The 
observations  could  then  be  extrapolated  to  analytical  con¬ 
ditions,  hopefully  with  no  changes  in  the  basic  underlying 
mechanism.  A  horizontal  profiling  system  was  used  which 
provided  radial  distributions,  after  Abel  inversion,  from 
the  channel  centre  to  5  to  6  mm  off  axis.  These  radial 
profiles  were  collected  at  three  different  observation 
h^icfhts^  7.5  mm,  15  mm,  and  30  mm  above  the  load  coil. 

For  Cal  (422.7  nm),  emission  enhancements  were  obser¬ 
ved  at  7.5  and  15  mm,  but  changed  to  a  depression  at  30  mm. 
The  analogous  absorption  measurements  showed  that  ground 
state  atom  number  densities  decreased  at  all  observation 
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heights.  A  similar  depression  was  observed  for  ground  state 
ion  (Call  393.3  nm)  populations  at  all  observation  heights. 
Ion  emission  was  enhanced  at  15  mm  but  depressed  at  both 
7.5  mm  and  30  mm  above  the  load  coil.  Kornblum  and  deGalan 
observed  that  this  data  indicated  that  there  were  probably 
three  major  contributions  to  the  effect  of  EIE's.  Radial 
data  showed  that  the  effect  of  Cs  was  to  shift  the  off  axis 
absorption  maximum  outward  from  the  central  channel.  They 
interpreted  this  as  a  volatilization  effect,  the  presence  of 
excess  matrix  resulting  in  a  slower  release  of  analyte  atoms. 
It  is  unclear  whether  this  was  the  result  of  larger  par¬ 
ticle  size  or  the  formation  of  more  stable  molecular 
species.  The  other  postulates  were  an  increase  in  the 
excitation  temperature  due  to  the  presence  of  Cs  resulting 
in  increased  excitation,  and  a  shift  in  ionization  equ¬ 
ilibrium,  enhancing  Cal  emission  and  depressing  Call  em¬ 
ission.  The  complex  interaction  of  these  three  effects  lead 
to  the  observed  behaviour.  The  fact  that  high  aerosol  flow 
rates  were  used  for  these  studies  may  make  some  of  the 
conclusions  suspect,  for  extrapolation  to  "normal"  operating 
conditions.  This  must  be  considered  especially  in  light  of 
observations  by  Kalnicky  et  al.  [136]  which  indicated  that; 
"The  surprisingly  small  changes  in  n^  and  temperature  pro¬ 
files  at  15  and  20  mm  upon  addition  of  Na  as  an  EIE  suggest 
that  changes  in  the  total  composition  of  the  sample  do  not 
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^■^tsct  the  radial  excitation  temperature  nor  degree  of 
ionization  of  analyte  species  in  a  dominant  manner". 

Clearly  more  spatial  studies  under  normal  operating  con¬ 
ditions  are  required. 

There  is  considerable  variability  in  the  observations 
of  workers  who  have  studied  the  effect  of  EIE's  on  analyte 
emission.  In  part  these  are  probably  due  to  a  considerable 
variability  in  operating  conditions  under  which  these  studie 
have  been  made.  The  observations  by  Boumans  and  DeBoer 
[46,47],  and  Larson  et  al.  [48],  indicate  that  the  mag¬ 
nitude  of  the  effect  is  a  function  of  power.  In  addition 
Boumans  and  DeBoer  demonstrated  that  a  change  in  aerosol 
flow  rate  could  change  the  effect  from  enhancement  to 
depression.  Couple  these  observations  with  those  of 
Kornblum  and  deGalan  [76],  who  observed  a  vertical  spatial 
dependence  whereby  enhancements  changed  to  depressions 
with  increasing  height,  and  at  least  part  of  the  problem 
is  clear.  It  seems  that  by  the  appropriate  choice  of 
operating  and  observation  conditions  any  matrix  effect 
desired  can  be  generated.  Of  course  this  would  not  be  done 
on  purpose  but  it  points  out  that  care  should  be  exercised 
interpreting  the  cause  of  the  effect  using  a  single  set 
of  conditions  for  the  investigation.  Also,  the  observation 
that  under  optimum  analytical  conditions  the  matrix  effect 
is  negligible  is  not  necessarily  true.  It  may  not  be 
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"observed",  but  instruments  which  measure  spatially  in¬ 
tegrated  emission  intensities  can  cloud  the  reality  of  the 
situation  [77] . 

Possible  Mechanisms  of  EIE  Interference 

Several  possible  mechanisms  have  been  postulated  for 
the  interference  effect  of  EIE's  on  analyte  emission  which 
has  been  surveyed  in  the  previous  section.  These  mechanisms 
manifest  themselves  at  different  states  in  the  analytes* 
progress  from  nebulization  to  emission  of  radiation.  A 
brief  discussion  of  each  of  these  is  warranted  as  a  basis 
for  a  discussion  of  the  data  presented  in  this  chapter. 

(1)  Shifts  in  ionization  equilibrium.  The  nature  of  this 
interference  has  already  been  discussed  in  the  previous  sec¬ 
tion.  If  ionization  equilibrium  shifts  were  a  dominant 
mechanism  in  the  ICP  we  would  expect  to  observe  the  same  type 
of  behaviour  as  that  exhibited  in  flames.  An  increase  in 
the  electron  density  of  the  plasma  would  shift  equilibrium 
toward  neutral  atom  species  resulting  in  enhancement  of 
atomic  emission  and  depression  in  ionic  emission.  The  de¬ 
gree  to  which  this  mechanism  will  influence  emission  in¬ 
tensities  in  the  plasma  will  depend  upon  how  much  the 
electron  number  density  in  the  plasma  is  changed  by  in¬ 
troduction  of  the  EIE.  It  has  been  calculated  [112]  and 
measured  [36]  that  the  electron  density  does  not  change 
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significantly  when  Na  is  introduced.  However,  the  measurement 
of  electron  density  extended  only  from  15  mm  to  25  mm  above 
the  load  coil.  An  extrapolation  of  the  plot  of  electron 
density  as  a  function  of  observation  height  to  lower  re¬ 
gions  indicates  that  this  may  not  be  the  case  between  0  mm 
and  15  mm  above  the  load  coil.  More  measurements  in  this 
region  are  warranted  before  generalizations  can  be  made. 

(2)  Enhanced  Collisional  Excitation.  Vellion  and  Margoshes 
[113]  suggested  enhancements  in  emission  intensity  may  be 
c^ue  to  increased  collision  rate  due  to  an  increase  in 
electron  density .  They  evoked  this  mechanism  to  explain 
increases  in  emission  intensity  due  to  added  Na,  which  did 
not  cause  corresponding  increases  in  the  ground  state  pop¬ 
ulation  of  the  analyte  species.  Franklin  et  al .  [85]  sug¬ 

gested  that  this  could  be  used  to  explain  simultaneous 
enhancements  in  the  emission  intensity  from  both  Ca  atom 
and  ion  lines  in  an  ICP .  If  this  is  the  case,  then  the 
effects  of  increased  collision  rate  due  to  increased  elec¬ 
tron  density  must  outweigh  the  effects  of  an  ionization  eq¬ 
uilibrium  shift.  Presumably  this  would  depend  upon  the  relative 
velocity  of  the  electrons.  The  greater  the  kinetic  energy 
of  the  electron,  the  more  likely  an  ion-electron  collision 
will  result  in  excitation  rather  than  recombination.  The 
outcome  of  an  increase  in  the  channel  electron  density,  then, 
will  depend  on  relative  kinetic  energies  of  the  electrons. 
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Kitagawa  and  Takeuchi  [122,144]  suggested  that  the  degree 
of  ionization  of  analyte  may  also  be  enhanced  by  increased 
collision  rate  due  to  added  EIE  in  a  microwave  plasma 
torch.  In  the  same  vein,  they  stated  that  the  degree  of  ion¬ 
ization  could  also  be  depressed  by  charge  transfer  re¬ 
actions  between  concomitant  atom  (c)  and  analyte  ions  (X+) 
in  the  following  manner, 

X+  +  c  Z  c+  +  X  (41) 

It  should  be  noted  that  energy  balance  requires  that  the 
ionization  potentials  of  the  two  species  closely  match 
each  other  for  this  process  to  take  place. 

(3)  Voltailization  effects.  Boumans  and  DeBoer  [47] ,  and 
Kornblum  and  deGalan  [76]  discussed  the  possibility  of 
volatilization  effects  in  the  ICP  as  a  mechanism  to  explain 
the  effects  observed  due  to  EIE's.  In  this  mechanism, 
analyte  atoms  are  "included"  in  a  desolvated  particle  matrix 
which  may  take  longer  to  volatilize  than  analyte  particles 
alone.  This  may  either  be  the  result  of  larger  particles, 
or  particles  with  atomic  bonds  which  are  more  difficult  to 
break.  The  net  effect  is  a  delay  in  the  release  of  the 
analyte  atoms,  when  compared  to  the  situation  without 
matrix.  This  delayed  volatilization  can  result  in  shifts 
in  the  spatial  behaviour  of  analyte  emission.  "Lateral 
diffusion  effects"  are  based  on  this  mechanism,  and  have  been 
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observed  in  flame  spectroscopy  [114,115],  and  discussed  by 
L  vov  et  al .  [116] .  In  the  ICP  a  lateral  diffusion  effect 

would  most  likely  result  in  a  depression  of  emission  in  the 
central  zone  of  the  plasma,  and  an  enhancement  nearer  the 
edges  of  a  radial  analyte  distribution,  in  lower  regions 
of  the  analyte  channel.  One  should  also  observe  vertical 
spatial  effects  as  a  result  of  volatilization  effects. 

These  should  manifest  themselves  as  a  depression  of  emission 
in  the  lower  analyte  channel.  Shifts  in  radial  emission 
profiles  [48],  and  absorption  profiles  [76]  in  the  ICP  have 
been  explained  as  due  to  volatilization  effects. 

(4)  Ambipolar  diffusion.  Boumans  and  DeBoer  [47]  have 
mentioned  ambipolar  diffusion  as  a  possible  source  of 
emission  intensity  modification  in  the  ICP.  Ambipolar 
diffusion  is  interpreted  as  the  diffusion  of  charged  par¬ 
ticles  under  the  influence  of  an  electric  field  [117] .  In 
a  gas  which  contains  an  appreciable  quantity  of  ions  and 
electrons,  electrons  tend  to  diffuse  faster  than  ions  which 
is  related  to  their  smaller  mass.  The  resultant  separation 
of  charge  produces  an  electric  field  which  exerts  its  force 
in  such  a  way  to  increase  the  drift  velocity  of  the  ion 
species  and  retard  the  drift  velocity  of  electrons.  Neutral 
particles  are  not  affected  to  any  significant  degree  by  this 
process.  The  diffusion  is  in  the  direction  from  higher  to 
lower  electron  densities.  Murayama  and  his  colleagues  [72,73] 
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have  interpreted  the  enhancement  of  emission  from  rare 
earth  elements  in  a  microwave  discharge  as  the  result  of 
ambipolar  diffusion. 

(5)  Nebulizer  effects.  The  process  of  nebulization  although 
simple  in  concept,  can  provide  several  different  mechanisms 
through  which  an  EIE  might  ultimately  affect  the  emission 
from  analyte  species.  The  most  straightforward  of  these 
are  solution  viscosity  effects.  An  increase  in  the  molar 
concentration  of  a  dissolved  salt  can  cause  increases  in 
the  viscosity  of  the  solution,  which  in  turn  will  decrease 
the  rate  of  transport  of  analyte  solution  to  the  nebulizing 
tip.  This  would  have  the  overall  result  of  decreasing  the 
amount  of  analyte  reaching  the  plasma  and  would  tend  to 
increasingly  depress  emission  signals  as  the  quantity  of 
salt  increased.  In  the  present  study  solution  uptake  was 
depressed  by  3%  at  NaCl  concentrations  of  2%,  and  the  changes 
from  this  source  are  not  considered  significant. 

In  a  recent  study,  Borowiec  et  al .  [118]  discussed  the 

phenomenon  of  "Aerosol  Ionic  Redistribution  (AIR)",  as  a 
mechanism  for  interference  due  to  concomitents  in  atomic 
spectroscopy.  They  have  shown  that  in  systems  involving 
cation  pairs,  the  nebulized  droplets  demonstrate  an  en¬ 
richment  in  the  minor  component  as  a  result  of  the  presence 
of  high  concentrations  of  the  major  component.  In  a  system 
involving  1%  Mg  and  Na  in  a  200:1  weight  ratio,  enrichments 
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as  high  as  220%  were  found  for  Na  in  the  aerosol  droplets. 
The  enrichment  is  generated  in  the  droplet  formation  stage 
and  the  degree  of  enrichment  is  related  to  the  droplet  size. 
The  enrichment  is  related  to  both  the  major  cation  con¬ 
centration  and  the  major  to  minor  cation  ratio  becoming 
negligible  at  ratios  less  than  50:1.  It  is  difficult  to 
assess  to  what  degree  this  effect  will  influence  obser¬ 
vations  of  the  effect  of  EIE's  in  the  plasma;  since  the 
magnitude  of  the  AIR  effect  is  related  to  droplet  size, 
it  will  depend  on  the  nebulizer  type,  and  operating  con¬ 
ditions.  The  effect  is  most  pronounced  in  very  small 
droplet  sizes  which  usually  form  a  low  percentage  of  the 
overall  droplet  distribution.  Dual  nebulizer  experiments 
by  Boumans  and  DeBoer  [47]  provide  some  evidence  that  AIR 
may  be  operative;  however  on  a  relative  basis  the  contri¬ 
bution  to  the  EIE  effect  from  the  nebulizer  was  small  com¬ 
pared  to  that  from  the  plasma.  Johnson  et  al.  [75]  found 
that  the  effect  of  added  EIE  in  a  DC  plasma  was  the  same 
when  analyte  and  concomitant  were  nebulized  from  separate 
nebulizers  or  from  the  same  nebulizer,  indicating  an  absence 
of  AIR. 

Skogerboe  and  Olsen  [119]  investigated  aerosol  trans¬ 
port  effects  due  to  the  presence  of  high  concentrations  of 
Na  in  a  microwave  induced  plasma  system  with  and  without  a 
desolvation  chamber  attached.  They  found  that  less  Pb  was 
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delivered  to  the  plasma  when  Pb  was  nebulized  in  the  pre¬ 
sence  of  Na  than  when  Pb  was  nebulized  alone  when  the  de¬ 
solvation  chamber  was  used.  They  attributed  the  decrease 
to  the  aerodynamic  properties  of  aerosol  particles,  the 
transport  efficiency  being  affected  mainly  by  gravitational 
settling  and  inertial  impingement  losses.  The  transfer 
efficiencies  for  wet  aerosol  systems  were  found  to  be  in¬ 
dependent  of  Na  concentration.  This  observation  was  sup¬ 
ported  in  a  subsequent  study  of  the  interference  of  high  Na 
concentrations  in  a  DC  plasma  system  [75] .  These  studies 
indicate  that  care  should  be  taken  in  comparing  results  from 
systems  employing  desolvation  and  those  that  do  not.  The  present 
study  did  not  involve  the  use  of  a  desolvation  chamber,  so 
it  is  expected  that  transport  effects  do  not  play  a  sig¬ 
nificant  role  in  causing  the  effects  observed. 

Although  there  may  be  other  mechanisms  which  may  con¬ 
tribute  to  the  EIE  effect,  those  which  have  been  mentioned 
are  those  which  have  most  frequently  appeared  in  the  lit¬ 
erature,  and  as  such  are  considered  the  prime  candidates. 

It  is  entirely  possible  that  no  single  mechanism  is  dominant, 
but  rather  the  observations  of  many  researchers  are  the 
result  of  a  complex  interaction  of  several  of  them.  This 
could  be  particularly  confusing  since  the  effects  of  one 
mechanism  may  tend  to  oppose  the  effects  of  another. 

Different  mechanisms  are  most  likely  to  manifest  them- 
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selves  in  different  spatial  regions  of  the  plasma.  It  has 
already  been  shown  that  different  spatial  regions  in  the 
plasma  display  different  types  of  emission  behaviour.  This 
spatial  inhomogeneity  could  result  in  different  spatial 
regions  having  matrix  effects  due  to  different  mechanisms. 

The  results  of  Kornblum  and  deGalan  suggest  this  may  be  the 
case  [76].  Thus  a  detailed  spatial  study  of  different 
analytes  and  interferents  under  various  operating  conditions 
has  been  undertaken.  The  purpose  of  the  study  presented 
in  the  following  pages  is  to  clarify  many  of  the  observations 
on  the  effect  of  EIE's  and  to  present  data  which  will, 
hopefully,  lead  to  a  "universal"  interpretation  of  this 
effect. 

Experimental 

The  elements  chosen  for  study  were  the  effects  of  Li, 

Na,  and  K  as  interferents  on  Ca,  Mg,  Cd,  and  Cr.  In  all 
cases  the  molar  concentrations  of  Ca,  Mg,  Cd,  and  Cr  were 
1.25  mmolar  corresponding  to  50,  30,  140,  and  65  ppm 
respectively.  Solutions  containing  these  elements  and  the 
interferent  in  the  form  of  a  chloride  salt  in  molar  ratios 
of  10,  30,  100,  and  300  to  1  were  made  up  using  distilled 
water.  These  represent  salt  concentrations  (by  weight)  in 
the  range  of  0.053  to  1.59%  for  Li,  0.073  to  2.19%  for  Na, 
and  0.077  to  2.8%  for  K.  Corresponding  blank  solutions 
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were  made  up  which  contained  the  interfering  salt  only, 
at  the  same  concentrations,  for  the  purpose  of  background 
subtraction . 

Effect  of  EIE's  on  Vertical  Profiles 

Vertical  spatial  profiles  have  been  acquired  in  the 
same  manner  as  those  presented  in  Chapter  III  in  the  region 
7  to  24  mm  above  the  load  coil  for  Ca,  Mg,  and  Cd  with  Na  as 
the  concomitent.  These  are  provided  in  Figures  44,  45,  and 
46.  The  numbers  on  each  set  of  figures  refer  to  the  molar 
ratio  of  Na  to  analyte;  these  are:  (1)  0:1,  (2)  10:1,  (3) 

30:1,  (4)  100:1,  and  (5)  300:1.  Curve  sets  a  and  b  refer  to 

emission  from  the  atom  and  ion  line  of  the  analyte  species 
respectively.  The  operating  conditions  chosen  for  the  study 
were,  1.5  kw  RF  power,  0.85  1pm  aerosol  flow  rate,  15  1pm 
plasma  gas  flow,  and  less  than  1  1pm  auxiliary  gas  flow.  These 
conditions  are  close  to  those  normally  used  for  analytical 
determinations  using  the  ICP. 

Figures  44a  and  44b  provide  the  vertical  spatial  pro¬ 
files  for  Cal  (422.7  nm)  and  Call  (393.3  nm)  respectively. 

It  can  be  seen  that  the  effect  of  Na  is  to  generally  enhance 
the  emission  of  the  atom  and  ion  lines  in  the  lower  regions 
of  the  aerosol  channel,  and  to  depress  emission  in  regions 
higher  above  the  load  coil.  The  magnitude  of  this  enhance¬ 
ment  or  depression  is  dependent  on  either  the  molar  ratio 
of  Na  to  Ca,  or  the  total  concentration  of  Na  present.  A 
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plot  of  the  magnitude  of  the  matrix  effect  as  a  function  of 
height  above  the  load  coil  is  provided  in  Figure  47  for  a 
molar  ratio  of  100:1  Na  to  Ca.  The  matrix  effect  is  de¬ 
fined  as  the  ratio  of  the  emission  intensity  with  Na  added 
to  the  emission  intensity  without  Na  added.  It  can  be  seen 
that  the  magnitude  of  the  effect  is  greatest  at  a  height  of 
7  mm  above  the  load  coil  and  decreases  with  an  increase  in 
height,  until  at  about  16  mm  the  signal  begins  to  show  a 
depression  which  continues  up  to  24  mm. 

Even  though  Cal  and  Call  behave  quite  similarly  at  a 
particular  height  above  the  load  coil,  the  response  at  the 
optimum  height  for  emission  for  the  two  is  quite  different. 
For  Cal,  enhancements  are  observed  at  the  emission  intensity 
peak.  Overall  the  profile  shifts  slightly  in  the  direction 
of  the  load  coil.  On  the  other  hand.  Call  is  depressed  at 
the  emission  intensity  peak,  the  profile  appearing  to  pivot 
about  a  crossover  point  where  there  is  no  effect  due  to 
excess  Na.  The  position  of  this  crossover  point  is  in¬ 
dependent  of  the  concentration  of  Na.  This  behaviour  appears 
to  be  the  result  of  the  two  lines  peaking  in  different  spa¬ 
tial  regions  in  the  ICP. 

Using  Figure  47  it  can  be  seen  why  variable  results 
might  be  reported  in  observing  the  effect  of  Na  on  Ca  using 
an  instrument  with  a  fixed  observation  height,  particularly 
one  with  a  relatively  small  slit  height.  In  this  case 
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Figure  47.  Matrix  effect  of  100:1  molar  ratio  Na  to  Ca  as 
a  function  of  height  above  the  load  coil. 
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Height  above  load  coil  (mm) 

Figure  48.  Matrix  effect  of  100:1  molar  ratio  Na  to  Mg  as 
a  function  of  height  above  the  load  coil. 
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Figure  49.  Matrix  effect  of  100:1  molar  ratio  Na  to  Cd  as 
a  function  of  height  above  the  load  coil. 
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measurements  made  low  in  the  plasma  (below  17  mm)  would 
result  in  the  observation  of  an  enhancement.  Measurements 
centered  about  17  mm  would  likely  result  in  an  observation 
of  "no  effect",  due  to  the  enhancement  effect  being  offset 
by  the  depression  effect,  and  measurements  above  17  mm  would 
indicate  a  depression  in  emission  intensity.  The  magnitude 
of  the  effect  would  depend  on  the  relative  size  and  pos¬ 
itioning  of  the  spatial  observation  window.  It  should  be 
noted,  also,  that  a  broad  spatial  window  (10-15  mm)  centered 
around  18-20  mm  would  likely  result  in  the  observation  that 
there  is  no  effect  of  easily  ionizable  elements,  which  might 
be  satisfying  analytically,  but  is  certainly  not  "true"  in 
this  case. 

Profiles  analogous  to  those  for  Ca  are  provided  in 
Figure  45  for  Mgl  (285.3  nm)  and  Mgll  (279.5  nm) ,  and  Figure 
46  for  Cdl  (228.6  nm)  and  Cdll  (214.4  nm) .  We  can  observe 
an  overall  trend  in  these  profile  sets  which  is  similar  to 
that  observed  for  Cal  and  Call.  In  general,  enhancements  are 
observed  at  lower  observation  heights,  changing  to  depres¬ 
sions  at  higher  observation  heights.  Atom  lines  and  ion 
lines  show  the  same  general  trends,  as  can  be  seen  in 
Figures  48  for  Mg  and  49  for  Cd.  These  are  plots  of  matrix 
effect  as  a  function  of  height  above  the  load  coil  analogous 
to  Figure  47,  for  a  molar  ratio  of  100:1  of  Na  to  analyte. 
These  plots  have  features  similar  to  those  for  Ca.  A 
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depression  is  not  observed  for  Cdl  below  24  mm  above  the 
load  coil,  but  the  trends  of  the  profiles  provided  in 
Figure  46  indicate  that  they  are  likely  to  be  encountered  at 
higher  observation  heights.  The  case  for  Mg  is  similar  to 
that  for  Ca.  The  atom  line  is  enhanced  at  its  peak  of 
emission  intensity  and  the  peak  position  shifts  slightly 
toward  the  load  coil.  The  ion  line  is  depressed  at  its 
peak  in  emission  intensity  and  also  exhibits  a  downward 
profile  shift.  The  crossover  point  for  Mgll  occurs  at  a 
higher  height  than  that  for  Call.  This  trend  continues  for 
Cdl  and  Cdll. 

There  are  several  features  of  the  data  presented  in 
Figures  44  to  49  which  warrant  consideration  in  an  inter¬ 
pretation  of  the  matrix  effect.  A  major  point  is  the  fact 
that  emission  from  both  atom  and  ion  lines  is  enhanced  at 
lower  observation  heights.  Even  though  the  overall  effect 
on  the  spatial  profile  is  different  between  these  two,  at  a 
fixed  height  above  the  load  coil  the  magnitude  of  the  effect 
is  very  similar  for  atom  and  ion  lines.  At  very  low  heights, 
ion  lines  tend  to  be  enhanced  more  than  atom  lines  of  the 
same  species.  Also  at  some  height,  the  matrix  effect 
changes  from  being  an  enhancement  to  a  depression  of  the 
analyte  emission  intensity.  The  plot  provided  in  Figure  50 
demonstrates  that  the  height  at  which  this  happens  for  ion 
lines  can  be  correlated  with  the  ionization  potential  of  the 
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Figure  50.  Spatial  position  where  enhancements  change  to 

depressions  as  a  function  of  ionization  potential 
for  Ca,  Mg,  and  Cd. 
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Figure  51. 


Ratio  of  emission  intensities  (I)  for  Call  393.3  nm 
and  Cal  422.7  nm.  Molar  ratio  of  Na  to  Ca;  (a)  0:1, 
(b)  100:1. 
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Figure  52.  Ratio  of  emission  intensities  (I)  for  Mgll  279.5  nm 
and  Mgl  285.3  nm.  Molar  ratio  of  Na  to  Mg;  (a) 

0:1,  (b)  100:1. 


Height  above  load  coil  (mm) 


Figure  53.  Ratio  of  emission  intensities  (I)  for  Cdll  214.4  nm 
and  Cdl  228.8  nm.  Molar  ratio  of  Na  to  Cd;  (a) 

0:1,  (b)  100:1. 
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element.  This  is  a  plot  of  the  spatial  position  of  the 
crossover  point  of  ion  line  enhancement  to  depression  as 
a  function  of  the  ionization  potential  of  the  element. 

Atom  lines  closely  follow  this  behaviour.  More  data  are 
required  to  confirm  the  correlation  but  the  trend  is  cer¬ 
tainly  present.  The  independence  of  spatial  position  of 
the  crossover  point  as  a  function  of  Na  concentration  is  a 
curious  feature  (Figures  44  to  46) . 

The  effect  of  Na  on  the  ratio  of  ion  emission  to  atom 
emission  is  plotted  in  Figures  51,  52,  and  53  for  Ca,  Mg, 
and  Cd  respectively.  The  data  for  these  plots  have  been 
extracted  from  Figures  44  to  46.  Curve  (a)  corresponds  to 
a  0:1  molar  ratio  of  Na  to  analyte,  and  curve  (b) ,  to  a 
100:1  molar  ratio  of  Na  to  analyte.  The  increase  in  the 
value  of  this  ratio  with  height  above  the  load  coil 
signals  an  increase  in  the  degree  of  ionization.  The 
magnitude  of  the  slope  in  each  plot  is  proportional  to  the 
ionization  potential  of  the  respective  analyte.  The  slope 
is  greatest  for  Ca  with  an  ionization  potential  of  6.11  ev, 
and  least  for  Cd  with  an  ionization  potential  of  8.99  ev. 

For  the  case  of  Ca,  the  addition  of  Na  depresses  the  ratio 
at  heights  less  than  17  mm,  and  enhances  it  at  higher 
heights.  This  depression  in  the  lower  region  of  the  channel 
might  be  the  result  of  some  small  amount  of  ionization 
suppression.  It  has  already  been  noted,  however,  that  overall, 
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enhancement  in  ion  emission  at  these  heights  is  observed. 
Therefore ,  even  if  ionization  equilibrium  shifts  are  taking 
place,  they  are  not  a  dominating  mechanism. 

In  contrast  to  the  slight  depression  observed  for  the 
ratio  of  intensities  for  Ca  up  to  17  mm,  the  ratio  is  en¬ 
hanced  for  Cd  and  Mg  at  heights  below  17  mm  due  to  excess 
Na.  It  seems  ionization  supression  does  not  even  play  a 
minor  role  for  these  two  elements.  The  shift  indicates 
either  an  increase  in  the  degree  of  ionization,  or  in¬ 
creased  excitation  of  the  ion  line  over  the  atom  line.  Above 
the  level  of  17  mm,  the  ratio  is  not  appreciably  affected 
by  the  addition  of  excess  Na. 

The  data  presented  in  this  section  certainly  do  not 
support  an  interpretation  of  the  matrix  effect  in  terms  of 
ionization  suppression  mechanism.  The  observations  for 
Cs  indicate  that  some  small  degree  of  this  must  be  taking 
place  lower  in  the  channel,  but  its  effect  is  dominated  by 
some  other  process  which  seems  to  affect  Ca,  Mg,  and  Cd  in 
a  similar  manner  for  both  atom  and  ion  emission.  The  fact 
that  enhancements  in  emission  are  observed  lower  in  the 
channel,  and  depressions  are  observed  in  the  upper  region, 
points  to  at  least  two  different  mechanisms,  or  one  mechanism 
which  redistributes  the  spatial  structure  of  analyte  emis¬ 
sion.  The  point  of  crossover  where  signal  enhancements  con¬ 
vert  to  depressions  is  more  dependent  on  analyte  ionization 
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parameters  than  the  concentration  of  Na  present,  even 
though  the  electron  density  at  this  point  should  increase 
with  increasing  Na  concentration.  This  might  indicate  that 
it  is  the  velocity  distribution  of  electrons  which  is  more 
important  than  the  actual  number  density,  velocity 
being  a  function  of  distance  from  the  load  coil.  Elements 
with  different  ionization  potentials  might  be  affected  in  a 
way  which  is  proportional  to  this  velocity  distribution. 

The  data  presented  thus  far  have  all  been  for  the 
interference  effect  of  Na  on  analyte  emission  intensity. 
Figures  54,  55,  and  56  demonstrate  that  similar  results  are 
obtained  for  the  matrix  effect  of  Li  on  Ca  and  Cr  emission. 
These  figures  are  analogous  to  Figures  44,  45,  and  46  for 
the  spatial  region  4  to  18  mm  above  the  load  coil.  Figure 
54  provides  the  effect  of  Li  on  Cal  (422.7  nm)  ,  Figure  55, 
the  effect  of  Li  on  Call  (393.3  nm)  ,  and  Figure  56,  the 
effect  of  Li  on  CrI  (357.8  nm)  and  CrII  (283.5  nm)  .  These 
figures  show  the  same  general  features  as  those  already 
presented.  To  investigate  any  subtle  changes  which  might 
ensue  when  the  EIE  is  changed, a  study  of  the  effect  of  Li, 
K,  and  Na  on  Ca,  Mg  and  Cd  has  been  carried  out. 

Effect  of  Varying  Interferent 

To  assess  the  effect  of  different  EIE's  on  the  mag¬ 
nitude  and  spatial  features  of  the  matrix  effect,  the  same 
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Figure  54.  Vertical  spatial  profiles  for  Cal  422.7  nm 
Molar  ratio  of  Li  to  Ca;  (a)  0:1,  (b)  10:1 

(c)  30:1,  (d)  100:1,  (e)  300:1. 
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Figure  55.  Vertical  spatial  profiles  for  Call  393.3  nm. 

ratio  af  Li  to  Ca;  (a)  0:1,  (b)  10:1,  (c)  30: 

(d)  100:1,  (e)  300:1. 
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Figure  56.  Vertical  spatial  profiles  for  Crl  357.8  nm  and  CrII  283.5  nm.  Molar  ratio 
of  Li  to  Cr;  (a)  0:1,  (b)  10:1,  (c)  30:1,  (d)  100:1,  (e)  300:1. 


Ca.,  Mg ,  and  Cd  atom  and  ion  lines  discussed  in  the  previous 
section  were  studied.  Spatial  profiles  in  the  region  7  to 
24  mm  above  the  load  coil  were  acquired  using  a  100:1 
molar  ratio  of  Li,  Na  and  K  to  analyte.  These  are  presented 
in  Figures  57,  58,  and  59  for  Ca,  Mg,  and  Cd  respectively. 
Curve  sets  1  and  2  represent  the  atom  and  ion  of  each  analyt 
species  respectively.  The  curves  are  labelled  in  the  fol¬ 
lowing  manner:  (a)  no  interferent,  (b)  Li  added,  (c)  Na 

added,  and  (d)  K  added.  The  conditions  are  1.5  Kw  RF 
power,  0.85  1pm  aerosol  flow  rate,  15  1pm  coolant  gas  flow 
rate,  and  less  than  1  1pm  auxiliary  gas  flow  rate. 

The  same  general  features  of  the  matrix  effect  are 
apparent  in  going  from  Li  to  Na  to  K  as  the  interfering 
element.  Analyte  emission  intensity  tends  to  be  enhanced 
lower  in  the  aerosol  channel  and  depressed  at  higher  heights 
The  magnitude  of  the  enhancement  increases  with  decreasing 
ionization  potential  of  the  EIE.  The  ionization  potentials 
of  Li,  Na,  and  K  are  5.39,  5.14,  and  4.34  ev  respectively. 
This  evidence  tends  to  support  the  notion  that  it  is  some 
species  produced  as  the  result  of  ionization  of  the  EIE 
which  produces  the  enhancements  observed,  either  electrons 
or  the  EIE  ion.  However,  the  effect  of  these  is  more  sig¬ 
nificant  in  the  lower  regions  of  the  plasma,  even  though  the 
eiE  should  be  more  fully  ionized  as  one  moves  up  the  analyte 
channel.  Either  the  relative  contribution  of  electrons  from 
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Figure  57.  Vertical  spatial  profiles  for;  (1)  Cal  422.7  nm,  and  (2)  Call  393.3  nm 
Concomitant  present  in  100:1  molar  ratio;  (a)  none,  (b)  Li,  (c)  Na , 
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the  EIE  to  the  overall  electron  density  decreases  with 
height ,  or  the  strength  of  the  interaction  of  electron  and 
RF  field  decreases  with  increasing  height.  Both  of  these 
are  probably  important. 

Effect  of  Power  on  Matrix  Effect 

It  was  pointed  out  in  the  introduction  that  several 
researchers  have  found  that  the  effect  of  increasing  the 
power  was  to  decrease  the  matrix  effect.  In  order  to 
clarify  this  observation  spatial  profiles  of  Cal  (422.7  nm) 
and  Call  (393.3  nm)  were  acquired  with  the  presence  of  a 
100:1  molar  ratio  Na  to  Ca,  and  without  the  presence  of  Na. 
The  spatial  observation  region  extends  from  7  to  24  mm 
above  the  load  coil.  The  conditions  used  for  this  study 
were  0.85  1pm  aerosol,  15  1pm  coolant  flow  and  less  than 
1  1pm  auxiliary  flow.  Figures  60 (i)  and  61 (ii)  provide  the 
spatial  profiles  for  Cal,  at  RF  powers  of  (a)  2.0  kw,  (b) 
1.75  kw,  (c)  1.50  kw,  (d)  1.25  kw,  and  (e)  1.0  kw.  The 
profiles  1  and  2  represent  the  emission  intensity  of  Ca 
without  and  with  the  presence  of  Na  respectively.  All 
profiles  are  on  the  same  scale.  Figures  61  (i)  and  61  (ii) 
provide  the  spatial  profiles  for  Call  in  the  same  format  as 
those  for  Cal.  The  plots  in  61  (ii),  frame  e,  are  on  a 
scale  which  is  4  times  as  sensitive  as  those  in  frames  a 
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The  basic  nature  of  the  interference  effect  does  not 
appear  to  change  as  the  power  is  changed  from  1.0  kw  to 
2.0  kw.  Figures  62  and  63  plot  the  magnitude  of  the  inter¬ 
ference  effect  as  a  function  of  applied  power  for  Cal  and 
Call  respectively.  The  data  for  these  figures  have  been 
taken  from  Figures  60  and  61.  Curves  a,  b,  c,  and  d  refer 
to  observation  heights  of  12,  14.5,  16.5,  and  21  mm  above 
the  load  coil  respectively.  It  can  be  seen  that,  indeed, 
the  effect  of  an  increase  in  power  is  to  decrease  the  mag¬ 
nitude  of  the  matrix  effect.  The  effect,  however,  does  not 
go  away.  It  is  still  apparent  at  all  power  levels.  A 
significant  observation  with  respect  to  an  interpretation 
of  the  decrease  in  matrix  effect  with  increasing  power  is 
the  variation  in  the  spatial  position  of  the  crossover  point 
where  enhancements  change  to  depressions.  As  power  is  in¬ 
creased  from  1.0  kw  to  2.0  kw  the  crossover  point  demon¬ 
strates  a  downward  shift  toward  the  load  coil.  If  one  were 
to  use  an  instrument  with  a  fixed  observation  height  and 
relatively  small  entrance  slit  height  to  perform  an  analogous 
experiment  the  results  could  be  misinterpreted.  For  example, 
if  the  observation  window  were  centred  at  15  mm,  the  en¬ 
hancement  would  decrease  as  RF  power  was  increased  in  a 
manner  similar  to  the  curves  provided  in  Figure  62.  At  2.0 
kw,  it  could  be  observed  that  the  matrix  effect  has  "dis¬ 
appeared".  In  fact  this  is  not  the  case.  At  this  height  the 


* 

' 


169 


Ajisuajuj  uoissjoig 


to  Ca 


170 


Ajisu94U|  uojssjujg 


171 


E 

E 


o 

U 


<D 

> 

O 

_jQ 

< 


g> 

'a> 

X 


^  03 
2 
o 

•  4-1 
CN  0 

—  0 
(C  "H 
-P 

03 

•V  U 

5-1 

03  P 
5  03 

0  i — i 

a  o 
E 
Eh 

PC  rH 

•  • 

o 

•  o 

E  rH 

c 

ro  cn 

• 

ro 

- 

m  re 
2 
H 

H  0 
03  2 
U 

0 

4-1 

U3  • 
Q)  £ 

rH  ^ 
•H 

4-1  O 

0  in 

U  • 

O-i  * — l 


03  U 
•H  *• — 

-P 

03  - 

a  2 

in 

iH  LD 

o3  r- 

u  •  • 

•Hi — I  03 
4->  U 

u  — 

03  JQ  0 
>  ~  -P 


•H 


4D 

03 

U 

a 

•H 

Eh 


X4isua^u|  uoissiujg 


172 


Xjisuejui  uoissiiug 


igure  61 (11) .  Vertical  spatial  profiles  for  Call  393.3  nm.  RF  power;  (d) 

(e)  1.0  kw.  (1)  No  Na,  (2)  100:1  molar  ratio  of  Na  to  Ca. 
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the  enhancements  lower  in  the  channel  coil  would  be  offset 
by  depressions  higher  up  to  give  the  impression  that  there 
is  no  matrix  effect  at  this  power.  For  Cal  these  ob¬ 
servations  could  further  be  compounded  by  the  overall  shift 
in  spatial  profile,  which  has  been  discussed  in  Chapter  III 
It  is  this  shift  which  produces  the  increase  in  effect 
between  1.25  and  1.5  kw  for  an  observation  height  of  12  mm. 
Clearly  the  spatial  information  is  required  for  a  complete 
interpretation  of  analyte  behaviour  in  the  ICP. 

Horizontal  Spatial  Profiles  of  the  Matrix  Effect 

To  assess  any  possible  radial  spatial  effects  of  the 
presence  of  excess  EIE  on  analyte  emission,  horizontal 
profiles  of  Ca  atom  and  ion  emission  have  been  acquired 
with  and  without  the  presence  of  a  100:1  molar  ratio  of  Na 
These  were  acquired  using  the  profiling  system  described 
in  Chapter  V  from  0  to  22  mm  above  the  load  coil,  in  2  mm 
increments.  The  corresponding  radial  emission  profiles  are 
Presented  in  Figures  64,  65,  66,  67,  68,  and  69  for  Cal 
(422.7  nm)  and  Call  (393.3  nm) .  In  each  figure  curve  (a) 
is  the  spatial  emission  profile  of  Ca  without  the  presence 
of  Na,  and  curve  (b)  is  the  profile  with  addition  of  Na. 

In  each  figure,  curve  (a)  is  given  an  amplitude  of  1.0  and 
curve  (b)  scaled  relative  to  this. 

The  spatial  profiles  for  Ca  at  0  mm,  2  mm,  4  mm,  and  6 
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Figure  64.  Radial  spatial  profiles  of  Ca  emission 

(i)  0  mm,  and  (ii)  2  mm  above  the  load 
(a)  no  Na  added,  (b)  100:1  molar  rati 
to  Ca . 
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Figure  65. 


Ca  I  422.7  nm  Ca  II  393.3  nm 


Radial  Coordinate  (mm) 


Radial  spatial  profiles  of  Ca  emission  at,  (i) 

4  mm  and  (ii)  6  mm  above  the  load  coil.  (a)  no 
Na  added,  (b)  100:1  molar  ratio  of  Na  to  Ca . 
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Figure  66. 


Ca  I  422.7  nm  Call  393.3  nm 


Radial  Coordinate  (mm) 


Radial  spatial  profiles  of  Ca  emission  at  (i) 
8  mm  and  (ii)  10  mm  above  the  load  coil. 

(a)  No  Na  added,  (b)  100:1  molar  ratio  of  Na 
to  Ca . 
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Figure  67.  Radial  spatial  profiles  of  Ca  emission  at,  (i) 
12  mm  and  (ii)  14  mm  above  the  load  coil.  (a) 
No  Na  added,  (b)  100:1  molar  ratio  of  Na  to 
Ca. 
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Figure  68. 


Radial  spatial  profiles  of  Ca  emission  at,  (i)  16 
mm  and  (ii)  18  mm  above  the  load  coil.  (a)  No  Na 
added,  (b)  100:1  molar  ratio  of  Na  to  Ca. 
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Cai 422.7  nm 


Ca  1393.3  nm 


Radial  Coordinate  (mm) 


Figure  GO. 


Radial  spatial  profiles  of  Ca  emission  at,  (i)  20 
mm  and  (ii)  22  mm  above  the  load  coil.  (a)  No  Na 
added,  (b)  100:1  molar  ratio  of  Na  to  Ca. 
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above  the  load  coil  are  provided  in  Figures  64  (i) ,  64  (ii)  , 
and  65  (i)  ,  and  65  (ii)  respectively.  The  behaviour  of  Ca 
atom  and  Ca  ion  in  this  region  is  very  similar.  An  off  axis 
emission  intensity  peak  is  exhibited  for  both  Cal  and  Call 
at  about  1  to  1.5  mm  from  the  channel  axis.  At  these 
heights,  the  presence  of  excess  Na  enhances  emission  intensity, 
but  does  not  substantially  alter  the  spatial  position  of 
the  peak  maximum.  This  is  an  indication  that  the  lateral 
diffusion  effect  from  delayed  volatilization  is  probably  not 
important  in  this  region. 

The  similar  behaviour  of  Cal  and  Call  in  this  region 

suggests  that  the  enhancement  is  probably  due  to  either  an 

increase  in  the  amount  of  analyte  transferred  to  the  plasma, 

or,  increased  excitation  due  to  the  presence  of  Na.  A  dual 

nebulizer  experiment  in  which  analyte  and  concomitant  were 

nebulized  separately,  then  mixed,  and  introduced  into  the 

plasma  demonstrated  that  similar  enhancements  are  observed. 

This  leads  one  to  believe  that  increased  collisional  excitation 

is  the  root  cause  of  the  signal  enhancements  observed  in 

this  region  of  the  plasma.  The  spatial  profiles  for  Ca  at 

8  mm  and  10  mm  above  the  load  coil  are  provided  in  Figures 

66 (i)  and  66 (ii)  respectively.  Analogous  enhancements  to 
those  observed  in  the  lower  region  are  observed  in  this 

region  for  both  Cal  and  Call.  Similar  off  axis  behaviour  is 

observed  for  both  Cal  and  Call  and  there  is  no  evidence  to 
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suggest  that  the  mechanism  causing  the  enhancement  in  this 
region  is  any  different  from  that  in  the  lower  region. 

In  the  region  12  mm  to  22  mm  above  the  load  coil  some 
interesting  differences  begin  to  appear  between  the  behaviour 
of  Cal  and  Call  emission,  when  Na  is  present.  At  a  height 
of  12  mm  above  the  load  coil  (Figure  67  (i) )  ,  Cal  no  longer 
displays  the  off  axis  peak  whether  or  not  Na  is  present. 

Call  however  retains  the  peak  at  this  height,  and  the  pre¬ 
sence  of  Na  shifts  the  position  of  maximum  intensity  more 
toward  the  edge  of  the  analyte  channel.  Whereas  at  this 
height  an  enhancement  is  observed  for  Cal  across  the  entire 
profile,  the  emission  for  Call  is  depressed  near  the  centre 
of  the  analyte  channel,  and  enhanced  at  distances  greater 
than  about  1  mm  from  the  centre.  Similar  behaviour  is  ob¬ 
served  at  14  mm  (Figure  6  7  ( ii ) )  and  16  mm  (Figure  6  8  ( i ) ) 
above  the  load  coil,  although  the  off  axis  peak  becomes  less 
pronounced  for  Call  at  these  heights.  This  seems  to  in¬ 
dicate  that  at  a  height  of  12  mm  above  the  load  coil  an 
additional  interference  mechanism  begins  to  make  itself 
evident  which  selectively  affects  Call  and  not  Cal.  The 
enhancements  for  Call  at  the  channel  edge  become  less 
pronounced  as  the  height  of  observation  increases  from  12 
to  16  mm  above  the  load  coil.  The  difference  in  spatial 
behaviour  between  Cal  and  Call  continues  at  heights  of  18 
mm  (Figure  68 (ii))  and  20  mm  (Figure  69  (i))  above  the  load 
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coil.  At  these  heights  the  emission  from  Cal  is  depressed 
at  all  radial  distances,  but  Call  still  demonstrates  de¬ 
pressions  at  the  channel  centre  and  enhancements  near  the 
edges.  The  spatial  profile  for  Cal  at  22  mm  above  the  load 
coil  is  not  available  but  clearly  the  trend  is  established. 

These  observations  suggest  that  at  observation  heights 
greater  than  12  mm,  ambipolar  diffusion  may  be  important  as 
a  matrix  effect  mechanism.  At  this  height  the  electron 
density  in  the  channel  from  ionization  of  Na  may  be  sig¬ 
nificant  compared  to  the  electron  density  at  the  channel 
edges  due  to  argon  atom  ionization.  Under  this  condition 
electrons  will  diffuse  toward  the  channel  edge  setting  up 
an  electric  field  which  will  "drag"  ions  along  with  them. 
This  reduces  the  overall  ion  population  at  the  channel 
center  and  increases  it  at  the  channel  edges,  resulting  in 
a  reduction  in  emission  intensity  in  the  centre  and  an 
enhancement  at  the  edge.  The  field  will  not  affect  neutral 
atom  populations  to  any  significant  degree  so  the  profile 
should  remain  relatively  unchanged.  The  data  presented 
in  Figures  67,  68,  and  69  support  this  interpretation. 

Significant  enhancements  of  emission  intensity  at  12 
mm,  14  mm,  and  16  mm  above  the  load  coil  indicate  that  the 
effects  of  increased  excitation  may  still  be  important,  but 
decrease  progressively  with  increasing  height.  At  18  mm, 

20  mm,  and  22  mm  above  the  load  coil,  ambipolar  diffusion 
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appears  to  be  the  dominant  mechanism.  The  depression  in 
Cal  emission  intensity  at  these  heights  is  probably  related 
more  to  an  excitation  effect  than  a  direct  diffusion  effect. 
Certainly  Cal  does  not  demonstrate  the  spatial  redistribution 
of  emission  intensity  that  is  observed  for  Call.  The  de¬ 
pression  may  be  the  result  of  an  enhancement  in  ionization. 

As  Call  ion  diffuses  outward  the  ionization  equilibrium 
will  shift  to  produce  more  ions  in  order  to  re-establish 
equilibrium.  This  will  reduce  the  overall  number  density 
of  Ca  atom.  The  relative  ratio  of  ion  to  neutral  atom 
populations  in  this  spatial  region  is  such  that  even  a  1% 
increase  in  the  degree  of  ionization  could  radically  reduce 
the  Ca  population,  without  "significantly"  adding  to  the 
ion  population. 

Nebulizer  Effects 

To  investigate  the  possibility  of  aerosol  ionic  re¬ 
distribution  as  a  mechanism  for  producing  the  enhancements 
observed,  an  experiment  was  performed  in  which  Mg  analyte 
and  Na  concomitant  were  independently  nebulized.  Mg  was 
nebulized  using  an  ultrasonic  nebulizer  and  Na  was  nebulized 
through  the  crossflow  nebulizer  used  in  the  previous  studies. 
The  aerosol  streams  were  mixed  and  introduced  into  the  plasma 
concurrently,  such  that  the  molar  ratio  of  Na  to  Mg  entering 
the  plasma  was  222  to  1 .  At  the  same  overall  flow  rate 
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the  mixed  solution  was  introduced  through  the  same  setup 
using  only  the  crossflow  nebulizer  at  molar  ratios  of 
100:1  and  300:1,  and  the  emission  intensity  for  Mg I  was 
interpolated  to  a  molar  ratio  of  222:1.  The  results  of  this 
experiment  showed  an  enhancement  of  1.87  for  the  dual 
nebulizer  system,  and  2.3  for  the  single  nebulizer  system, 
at  a  height  of  10  mm  above  the  load  coil.  These  results 
indicate  that  under  these  conditions  23%  of  the  observed 
enhancement  could  be  due  to  an  AIR  effect.  This  is  sig¬ 
nificant  but  is  still  not  sufficient  to  completely  explain 
the  observed  enhancements,  and  is  certainly  not  the  dominant 
interference  mechanism. 


Ionization  Equilibrium  Shifts 

There  is  very  little  indication  that  any  of  the  effects 
observed  result  from  a  shift  in  the  ionization  equilibrium 
of  the  analyte  species.  This  has  been  confirmed  by  perform 
ing  an  atomic  absorption  experiment  on  Mgl  under  plasma 
operating  conditions  and  in  a  spatial  region  where  large 
signal  enhancements  are  observed  due  to  the  presence  of  Na. 
The  measurement  was  made  by  monitoring  the  285.3  nm  Mgl 
line  from  a  hollow  cathode  lamp  focussed  on  the  analyte 
channel.  Measurements  were  performed  using  a  photomultiplier 
tube  coupled  to  the  data  acquisition  system  described  in 
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Appendix  A.  Various  molar  ratios  of  Na  to  Mg  were  in¬ 
troduced  and  the  emission  and  absorption  signals  measured. 
The  results  are  summarized  in  Table  IX.  It  can  be  seen  that 
while  the  emission  signal  increases  with  increasing  Na 
concentration,  the  absorption  signal  shows  only  a  slight 
change  indicating  an  increase  in  excited  state  population 
but  no  corresponding  increase  in  ground  state  population 
of  the  Mg  atom.  This  data  supports  an  interpretation  of 
the  effect  in  terms  of  increased  collisional  excitation,  and 
discounts  any  mechanism  which  increases  overall  atom 
populations . 

Summary  and  Conclusions 

Clearly  the  spatial  study  of  the  effect  of  EIE's  on 
analyte  emission  intensity  has  clarified  the  observations 
of  many  workers  in  the  field,  and  has  provided  data  to 
support  a  rational  interpretation  of  the  nature  of  the 
effect.  This  type  of  spatial  information  is  indespensable 
in  a  description  of  the  matrix  effect.  It  has  been  demon¬ 
strated  that  the  magnitude  and  net  result  of  the  presence 
of  excess  EIE  is  spatially  dependent.  This  spatial  de¬ 
pendence  has  been  used  to  satisfactorily  explain  the  vari¬ 
ability  in  observations  by  several  investigators  on  the 
same  analyte-concomitent  system.  The  matrix  effect  does 
not  disappear  under  an  appropriate  choice  of  operating  and 
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TABLE  IX 

EMISSION  INTENSITY  AND  ABSORBANCE  OF 
Mg I  285.3  nr:  FOR  VARIOUS  MOLAR  RATIOS  OF  Na  to  Mg 


Solution 

Emission 

Intensity 

Absorbance 

0  : 1 

Na :  Mg 

27.14 

0.2805 

10:1 

Na :  Mg 

38.63 

0.3129 

30:1 

Na :  Mg 

49.95 

0.3132 

100:1 

Na  :  Mg 

69.52 

0.3170 

300  : 1 

Na  :Mg 

95.51 

0.2929 
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observation  conditions.  The  spatial  redistribution  of 
analyte  emission  can  make  this  appear  to  be  the  case;  how¬ 
ever,  it  is  simply  not  an  "exact"  description  of  the  situation. 

In  general,  the  addition  of  excess  EIE  enhances  em¬ 
ission  in  the  lower  regions  of  the  analyte  channel  for  both 
atom  and  ion  lines.  In  the  upper  regions  of  the  analyte 
channel,  emission  intensity  is  depressed  for  both  atom  and 
ion  species.  The  spatial  position  at  which  enhancements 
switch  over  to  depressions  is  not  a  function  of  Na  con¬ 
centration  but  rather  seems  correlated  to  the  ionization 
potential  of  the  analyte  species.  Thus  it  is  felt  that 
the  addition  of  Na  defines  the  excitation  environment,  and 
different  analyte  species  respond  in  slightly  different 
ways  to  this  environment.  The  cause  of  the  crossover 
seems  to  be  the  result  of  a  switch  in  the  mechanism  of 
interference.  In  the  lower  regions  of  the  channel,  the 
mechanism  causing  an  enhancement  appears  to  be  the  result 
of  increased  collisional  excitation,  as  a  result  of  an 
increase  in  the  number  of  electrons  with  the  appropriate 
energy  to  cause  excitation  of  atom  or  ion  levels.  The 
energy  of  these  electrons  will  be  a  function  of  distance 
from  the  load  coil.  Thus,  the  decrease  in  overall  en¬ 
hancement  with  increasing  height  above  the  load  coil.  This 
interpretation  also  correlates  well  with  the  observation 
that  the  enhancement  increases  slightly  as  the  ionization 
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potential  of  the  concomitant  decreases  .  The  effect  of  increasing 
coilisional  excitation  appear  to  outweigh  any  effects  of  a 
shift  in  ionization  equilibrium  which  could  take  place  with 
an  increase  in  electron  density.  The  evidence  from  the  AA 
experiment  on  Mgl,  confirms  the  relative  absence  of  any  such 
ionization  shift  mechanism  at  heights  where  emission  en¬ 
hancements  occur.  As  the  effects  of  increased  excitation 
decrease  with  height,  the  influence  of  ambipolar  diffusion 
begin  to  exert  an  influence  on  the  spatial  distribution  of 
analyte  species.  This  eventually  leads  to  the  observed 
depressions  in  the  analyte  emission.  There  is  not  enough 
data  to  interpret  the  plot  provided  in  Figure  50,  which 
indicates  the  position  of  the  switchover  as  related  to 
ionization  potential  of  the  analyte.  This  may  be  related 
to  the  relative  width  of  the  radial  profiles  of  Ca,  Mg,  and 
Cd  at  the  various  heights.  Certainly  a  detailed  radial 
spatial  study  of  a  variety  of  analyte  species  and  concom¬ 
itants  is  warranted  to  further  clarify  this  point. 

The  important  message  of  this  chapter  is  that  workers 
should  take  care  in  investigating,  reporting,  and  interpre¬ 
ting  this  complex  system.  Much  confusion  could  have  been 
avoided  if  early  investigators  had  more  carefully  studied 
the  spatial  dependence  of  the  matrix  effect.  Different 
observations  by  different  workers  gave  the  impression  that 
the  nature  of  the  effect  varied  from  plasma  to  plasma.  As 
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has  been  stated  previously,  it  is  the  author's  opinion  that 
the  interaction  of  argon  plasma  gas,  analyte,  and  con¬ 
comitant  is  basically  the  same  in  all  argon  ICP's, 
regardless  of  RF  power,  aerosol  flow  rate,  etc. 


CHAPTER  VII 


CONCLUSIONS 

The  ease  and  speed  with  which  spatial  information  may  be 
acquired  using  a  photodiode  array  detector  has  facilitated  a 
detailed  study  of  ICP  spatial  structure.  This  structure  is 
a  complex  function  of  instrumental  operating  parameters,  and 
analyte  excitation  and  ionization  characteristics.  The  emis¬ 
sion  lines  studied  divide  into  two  categories,  based  on  re¬ 
sponse  of  the  spatial  profile  to  changes  in  applied  power. 

The  first  of  these  is  one  in  which  the  position  of  maximum 
emission  intensity  is  highly  dependent  on  RF  power.  These 
are  predominantly  atom  lines  whose  emission  peaks  in  the 
spatial  region  below  15-17  mm  above  the  load  coil  (at  1.5  kw 
power  and  0.8  to  0.9  1pm  aerosol  flow  rate).  The  vertical 
position  at  which  these  lines  peak  in  the  analyte  channel  can 
be  correlated  with  their  normal  temperature.  The  normal 
temperature  can  be  easily  calculated  using  a  set  of  equations 
if  LTE  is  assumed.  The  ability  to  predict  spatial  structure 
based  on  an  LTE  calculation,  suggests  that  in  this  spatial 
region  LTE  may  exist.  This  region  is  not  normally  observed  in 
analytical  applications.  Thus  it  has  not  been  studied  in  as 
much  detail  as  the  region  between  15  and  25  mm  by  workers 
using  temperature  and  electron  density  measurements  as  probes 
to  test  for  LTE.  This  lower  region  should  be  further  studied 
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since  it  is  fundamentally  very  interesting.  Such  studies 
would  confirm  or  deny  the  existence  of  LTE .  Lines  whose 
emission  peaks  between  17  and  25  mm  above  the  load  coil  form 
the  other  category.  The  spatial  structure  of  these  lines  is 
not  as  highly  dependent  on  ICP  operating  parameters  or 
analyte  excitation  characteristics.  It  has  been  stated  that 
LTE  does  not  exist  in  this  region  and  the  spatial  structure 
seems  to  support  this. 

The  upper  and  lower  regions  are  characterized  by  two 
different  types  of  structure  for  radial  emission  profiles  as 
well.  In  the  lower  region,  between  0  and  12  to  15  mm  above 
the  load  coil,  analyte  lines  show  off  axis  peaks.  These  peaks 
roughly  mark  the  edge  of  the  aerosol  channel.  Excitation  in 
this  region  seems  mainly  the  result  of  contact  with  hot  argon 
gas  originating  from  near  the  channel  edge.  Higher  up,  this 
structure  collapses  and  emission  peaks  in  the  centre  of  the 
channel.  It  is  likely  that  at  these  heights  there  is  complete 
mixing  between  analyte  and  hot  argon  gas  originating  from 
annular  regions  of  the  plasma.  The  width  of  the  radial 
profiles  for  ion  lines  is  greater  than  those  for  atom  lines. 
This  may  be  the  result  of  ambipolar  diffusion  of  ion  species. 
The  relationship  between  profile  width  and  analyte  species 
(ion  or  atom)  and  excitation  characteristics  should  be  in¬ 
vestigated.  It  may  be  that  the  degree  of  ionization  increases 
as  radial  distance  increases.  This  would  also  broaden  the  ion 
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line  relative  to  the  neutral  atom  line.  A  detailed  study  of 
ion  and  atom  ground  state  populations  would  certainly  aid  in 
an  overall  understanding  -  of  this  spatial  structure.  In  fact, 
three  dimensional  spatial  maps  of  ground  state  populations 
analogous  to  those  for  emission  presented  in  Chapter  V  would 
be  invaluable  in  characterizing  ICP  phenomena.  Included  in 
these  should  be  a  map  of  argon  metastable  densities.  To 
this  point  spatial  studies  of  the  aerosol  channel  have  sug¬ 
gested  that  these  species  reach  maximum  density  at  heights 
between  18  and  22  mm  above  the  load  coil.  A  spatially  re¬ 
solved  absorption  experiment  would  confirm  or  deny  this  sug¬ 
gestion  . 

The  dependence  of  emission  intensity  on  excess  EIE  has 
certainly  been  clarified  by  the  spatial  study.  It  has  been 
demonstrated  that  many  of  the  conflicting  statements  which 
have  appeared  on  the  nature  of  the  effect  can  be  rationalized 
on  the  basis  of  changes  in  the  spatial  structure  of  analyte 
emission.  Several  suggestions  have  been  offered  on  the 
reasons  for  these  changes  in  spatial  structure.  The  enhance¬ 
ments  in  emission  observed  in  the  lower  regions  of  the  analyte 
channel  appear  to  be  the  result  of  an  increase  in  the  channel 
electron  temperature.  Measurement  of  electron  temperature 
from  the  absolute  continuum  emission  intensity,  with  and  without 
the  presence  of  EIE,  would  clarify  this  suggestion.  De¬ 
pressions  in  emission  intensity  appear  to  be  the  result  of 
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ambipolar  diffusion.  The  dependence  of  the  position  where 
enhancements  convert  to  depression  on  analyte  ionization 
characteristics  is  significant  and  should  be  further  in¬ 
vestigated.  This  may  be  related  to  radial  profile  width,  atom 
and  ion  mobilities,  or  complex  electric  field  effects. 

There  are  several  areas  which  warrant  spatial  in¬ 
vestigation  which  have  not  been  undertaken  as  part  of  this 
thesis.  The  spatial  distribution  of  species  introduced  as 
part  of  the  sample  aerosol  is  one  of  these.  The  dissociation 
products  of  water  such  as  OH,  H,  and  oxygen  radical  may  be 
important  in  excitation  and/or  de-excitation  processes  and 
should  be  studied  in  detail.  This  should  be  combined  with 
a  study  of  any  spatial  redistribution  of  analyte  emission 
which  may  result  from  the  use  of  a  desolvation  chamber.  In 
addition,  it  is  curious  that  as  ion  emission  drops  off  in  the 
area  between  20  and  35  mm  above  the  load  coil,  atom  emission 
does  not  correspondingly  increase.  This  may  be  the  result 
of  oxide  formation.  A  spatial  study  of  analyte  oxides  would 
be  helpful  in  deciding  if  these  species  are  important  in  the 
removal  of  atom  and  ion  species. 

Finally,  as  has  been  illustrated  a  number  of  times  through¬ 
out  this  thesis,  a  detailed  knowledge  of  analyte  spatial 
behaviour  has  clarified  observations  which  have  been  made  by 
researchers  using  instruments  which  measure  spatially  integrated 
emission  intensities.  It  has  been  pointed  out  that  care  should 
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be  exercised  in  interpreting  plasma  phenomena  based  on 
measurements  of  this  type  alone.  The  confusion  on  the  nature 
of  the  matrix  effect  is  a  good  example  of  this.  Spatially 
resolved  measurements  are  the  key  element  in  the  investigation 
and  understanding  of  a  complex  source  like  the  ICP .  The  use 
of  a  linear  photodiode  array  is  a  rapid,  detailed,  and 
quantitative  method  for  making  these  measurements. 
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APPENDIX  A 


The  photomultiplier  tube  is  the  most  widely  used  de¬ 
tector  for  spec trochemical  measurements  in  the  ultraviolet  and 
visible  spectral  regions.  Although  many  diverse  experimental 
situations  exist, a  rather  large  number  of  applications  have 
very  similar  basic  measurement  requirements.  In  our  own  lab¬ 
oratory  we  had  a  need  for  a  general  purpose  wide  dynamic  range 
integrating  readout  system.  One  of  the  main  applications  for 
this  system  was  the  acquisition  of  basic  signal  to  noise  ratio 
data  on  spec trochemical  sources  such  as  the  inductively 
coupled  plasma.  For  this  application  background  subtraction 
capability  was  required  in  order  that  a  blank  (i.e.  plasma 
spectral  background)  and/or  photomultiplier  tube  dark  current 
signal  could  be  automatically  subtracted  from  an  analyte 
emission  reading.  In  addition,  to  facilitate  the  signal  to 
noise  ratio  studies,  automatic  acquisition  of  replicate  in¬ 
tegrated  measurements  was  required  along  with  subsequent 
calculation  of  standard  deviations  and  signal  to  noise 
ratios.  Finally  it  was  desirable  to  have  all  data  logged  in 
a  systematic  fashion  to  a  hardcopy  readout  device. 

These  basic  requirements  essentially  dictated  the  de¬ 
velopment  of  a  microprocessor  based  readout  system.  In  fact 
we  developed  two  types  of  wide  dynamic  range  integrating 
readout  systems  both  of  which  operate  under  microprocessor 
control.  In  both  systems  the  output  current  signal  from  the 
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PM  tube  is  first  converted  to  a  voltage  with  an  operational 
amplifier.  Then,  in  the  first  measurement  system  (Figure  70) 
this  voltage  is  converted  to  a  frequency  with  a  wide  dynamic 
range  (6  orders  of  magnitude)  voltage-to-f requency  converter. 
The  voltage-to-f requency  converter  was  interfaced  to  an  8080 
based  single  board  microcomputer,  the  Intel  SDK-80.  The  key 
subsystem  in  this  interface  is  the  8253  programmable  interval 
timer  chip  of  Intel.  Utilizing  two  of  those  chips  in  con¬ 
junction  with  the  SDK-80  microcomputer  one  has,  in  a  sense,  a 
very  flexible  software  programmable  frequency  meter  for  the 
measurement  of  v-to-f  converter  output.  In  the  second  system 
(Figure  70)  a  monolithic  10  bit  successive  approximation  ADC 
was  interfaced  to  the  SDK-80.  Wide  dynamic  range  integrating 
measurement  capability  was  achieved  by  summing  a  large  number 
of  conversions. 

In  both  systems  the  microprocessor  controls  the  ac¬ 
quisition  of  replicate  integrated  measurements  and  software 
is  included  for  background  subtraction,  the  calculation  of 
standard  deviations  of  the  replicate  values  and  signal  to 
noise  ratios.  Development  of  this  software  capability  in  a 
small  system  was  greatly  facilitated  by  use  of  the  AM9511 
Arithmetic  Processing  Unit  of  Advanced  Micro  Devices.  This 
chip  (24  pin  package)  is  capable  of  performing  fixed  point 
single  and  double  precision  (16  or  32  bit)  arithmetic  and 
floating  point  single  precision  (24  bit  mantissa,  7  bit  ex- 
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Figure  70.  Block  diagram  of  the  voltage-to-f requency  con¬ 
verter  and  analog-to-digital  converter  based 
measurement  systems. 
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ponent)  arithmetic.  Operations  performed  include  add, 
subtract,  multiply,  divide, trigonometric  and  inverse  trig¬ 
onometric  functions,  square  roots,  logarithms  and  exponen¬ 
tiation.  The  use  of  this  chip  greatly  expedited  overall 
software  development  and  execution  time  compared  to  soft¬ 
ware  floating  point  packages.  Both  systems  will  now  be 
described  in  detail. 

Experimental 

Components  common  to  both  measurement  systems  are  the 
operational  amplifier  and  the  SDK-80  single  board  micro¬ 
computer.  The  OA  is  Model  52K  (Analog  Devices),  a  low  noise, 
low  drift,  precision  FET  amplifier.  It  was  operated  with  a 
10  Mfi  feedback  resistor  in  parallel  with  a  450  pf  capacitor. 

-7 

Thus  the  output  voltage  was  1  V  for  a  10  amp  input  cur¬ 
rent  from  the  PM  tube. 

The  heart  of  both  measurement  systems  is  the  Intel 
SDK-80  single  board  microcomputer.  The  basic  hardware 
features  of  this  board  include  an  8080  A  central  processing 
unit,  1  k  of  read/write  memory  (8111's),  4  k  of  2708 
eraseable  programmable  read  only  memory  (EPROM) ,  an  8251 
programmable  peripheral  interface  chip  providing  48  digital 
I/O  lines.  In  addition  a  small  but  very  useful  breadboard¬ 
ing  area  is  provided  on  the  board.  All  additional  circuitry 
for  both  applications  was  placed  in  this  breadboard  area, 
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except  for  the  voltage-to-f requency  converter  and  op¬ 
erational  amplifier,  which  were  located  near  the  PMT . 

Voltage-to-Freauency  Converter  Based  System 

The  voltage-to-f requency  converter  used  was  Model 
460  k  (Analog  Devices) ,  a  high  accuracy  unit  with  a  six 
decade  (1  Hz  to  1  MHz)  dynamic  range. 

The  key  subsystem  of  this  measurement  system  is  the 
Intel  8253  programmable  interval  timer.  A  block  diagram  of 
this  24  pin  integrated  circuit  is  shown  in  Figure  71.  It 
consists  of  three  independent  16  bit  pre-settable  down 
counters,  read/write  logic,  a  control  word  register  and  a 
data  bus  buffer.  The  read/write  logic  controls  the  overall 
device  operation  such  as  loading  and  reading  of  the  various 
counters  as  well  as  loading  of  the  control  word  register. 

The  counters  can  be  software  programmed  to  act  as  either 
binary  or  BCD  counters  and  each  counter  can  be  independently 
programmed  to  operate  in  any  of  six  different  modes.  The 
basic  operation  of  each  mode  is  illustrated  schematically 
in  Figure  72. 

Mode  0  is  a  software  triggered  interupt  in  which  an 
interupt  signal,  a  rising  edge  at  the  output  terminal,  is 
generated  N  clock  pulses  after  software  loading  the  counter 
with  N.  In  Mode  1  (programmable  one  shot  or  monostable)  an 
external  trigger  (rising  edge)  applied  to  the  gate  input 
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Block  diagram  of  the  8253  programmable  interval 
timer . 
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Schematic  representations  of  the  8253  modes 
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generates  a  zero  level  pulse  at  the  output.  The  duration  of 
the  pulse  is  N  clock  pulses  where  N  is  again  a  software 
loaded  value.  Mode  2  and  Mode  3  both  divide  by  N 
counter  modes.  In  Mode  2  the  output  is  a  zero  level  pulse, 
one  clock  pulse  in  duration, every  Nth  clock  pulse  while  in 
Mode  3  the  output  is  a  symmetrical  divided  by  N  waveform. 

Mode  4  provides  a  software  triggered  strobe,  a  zero  level 
pulse  of  one  clock  period  duration  that  occurs  N  clock 
pulses  after  a  software  load  of  the  counter.  Mode  5 
(hardware  triggered  strobe)  provides  a  similar  operation 
but  requires  an  external  trigger,  a  rising  edge  at  the 
gate  input.  The  above  discussion  and  Figure  72  provide 
only  a  brief  and  basic  description  of  the  operation  and 
capabilities  of  the  8253.  For  complete  details  Intel's 
literature  should  be  consulted  [120] . 

The  complete  schematic  of  the  way  the  8253 's  are  connected 
to  the  SDK-80  is  shown  in  Figure  73.  One  of  the  8253' s 
functions  as  the  counter  for  the  measurement  of  the  voltage- 
to-frequency  converter  output  and  the  second  8253  generates 
the  counter  time  base  (integration  time)  for  the  first  8253. 
Simplified  diagrams  of  the  integration  timer  and  counter 
are  shown  in  Figure  74. 

The  integration  time  is  derived  from  the  2.048  MHz 
clock  of  the  microprocessor.  Counter  0  on  the  timer  8253 
is  operated  in  Mode  2,  binary  (divide  by  N  counter)  by 
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Circuit  diagram  of  the  complete  programmable 
frequency  meter  as  interfaced  to  the  .SDK-80. 
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loading  HEX  3C  into  the  control  register.  This  counter 
is  then  preset  to  HEX  5004  which  divides  the  2.048  MHz 
system  clock  down  to  100  Hz  which  is  available  at  the  output 
pin  of  counter  0.  This  100  Hz  signal  is  connected  to  the 
clock  input  of  the  second  counter.  Counter  1  (the  second 
counter)  on  the  same  IC  is  programmed  to  operate  in  Mode 
5  BCD  (hardware  triggered  strobe)  by  loading  the  control 
register  with  HEX  BB,  and  is  present  with  a  BCD  number  input 
from  the  teletype.  For  example,  presetting  this  counter 
with  0100  results  in  a  1  second  integration  time.  Thus 
integration  times  from  0.01  to  99.99  sec  in  0.01  sec 
intervals  can  be  set  by  TTY  input. 

The  other  8253  is  used  to  count  the  pulses  produced 
by  the  voltage  to  frequency  converter  for  the  integration 
time  generated  by  the  timer  8253.  In  order  to  obtain  the 
desired  resolution  two  counters  on  this  second  chip  are 
cascaded,  providing  32  bit  resolution.  Both  16  bit  cou¬ 
nters  are  set  to  operate  in  Mode  2,  binary  by  loading  the 
control  register  with  HEX  BC  and  7C.  The  counters  are 
preset  with  HEX  FF  for  full  count. 

An  8255  programmable  peripheral  interface  chip  that  is 
a  standard  feature  of  the  SDK-80  board  is  used  to  aid 
sequencing  of  the  measurement.  A  measurement  is  initiated 
by  outputting  a  logic  0  followed  by  a  logic  1  at  one  bit 
of  Port  B.  This  sets  the  flip-flop  (Q  =  1)  which  enables 
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the  gates  of  both  the  timer  and  counter  8253 's.  Note  that 
both  8253' s  have  NAND  gates  on  the  cascade  line  from  output 
to  the  succeeding  input  of  each  IC.  This  is  necessary 
because  we  found  that  the  8253 's  do  not  count  the  first 
pulse  which  is  used  to  set  up  the  counter.  For  this  reason 
a  charge  delay  monostable  is  used  to  set  up  the  counters 
at  the  start  of  the  measurement.  The  charge  delay  mono¬ 
stable  is  triggered  by  the  falling  edge  at  Q  of  the  flip- 
flop  . 

When  the  timer  8253  reaches  its  full  count  its  output 
undergoes  a  1  0  transition  which  toggles  the  flip-flop 

setting  Q  to  "0"  thus  disabling  all  gates  and  terminating 
the  count  measurement.  During  counting  the  status  of  the 
flip-flop  is  continuously  monitored  by  inputting  Q  of  the 
flip-flop  to  bit  0  of  Port  A  on  the  8255.  When  counting 
is  completed  (Q  =  "1")  the  microcomputer  enters  a  routine 
that  reads  and  stores  the  binary  numbers  in  the  counter  8253 
in  memory. 

Successive  Approximation  ADC  Based  System 

The  complete  schematic  of  the  successive  approximation 
ADC  based  system  is  shown  in  Figure  75.  An  Analog  Devices 
AD571  was  chosen  as  the  ADC.  It  is  a  10  bit  unit  with  a 
specified  conversion  time  of  25  psec.  The  ADC  was  front 
ended  with  an  Analog  Devices  AD582  sample-and-hold  amplifier. 


, 


AD  571  8255 


218 


SH8  SS3baav  08-  xas 


sng  lOb'iNOD  G8-xa$ 


sn8  vivaoe-xas 


u 

o 

4-> 

c 

•H 

CO 

(3 

C 

o 

4-) 

to 

to 

T, 

CD 

CO 

rd 

X! 

U 

a 

4-> 

u 

Q) 

> 

C 

o 

u 


03 

4-> 

-H 

cn 

-H 

'C 

I 

0 

4J 

I 

tT' 

0 

rH 

rd 

c 

rd 

CD 

XT 

4-J 

u_ 

0 

£ 

fd 

c 

c 

•H 

t: 

4-> 


O 

CO 

I 

in 

<D 

x: 

4-) 

o 

4J 


D  T3 

a  a) 

^  u 

■H  fd 

U  4-1 


LT) 


O 

^4 

P 

Co 

*r4 

En 


219 


The  sample-and-hold  amplif ier-ADC  subsystem  was  interfaced 
to  "the  SDK-80  via  the  8255  programmable  peripheral  interface 
chip.  Acquisition  of  a  data  point  was  initiated  by  a  pulse 
from  pin  C0  of  Port  C  and  upon  receipt  of  the  DATA  READY 
signal  the  converted  10  bit  value  is  read  into  the  micro¬ 
processor  using  Ports  A  and  B  of  the  8255.  Integration  is 
achieved  by  summing  a  large  number  of  converted  values.  For 
example  a  1  sec  integration  is  achieved  by  summing  10,923 
conversions.  Even  though  we  are  using  a  10  bit  converter 
the  integrated  value  has  a  resolution  of  about  16-17  bits 
as  a  result  of  the  summation  (signal  averaging)  step  [121]. 

Software  Development 

The  key,  of  course,  to  the  versatility  and  power  of  a 
microprocessor  based  system  is  software  control.  The  SDK-80 
single  board  microcomputer  has  room  for  1  k  of  read/write 
memory  (RAM)  and  4  k  of  ROM  (i.e.  4  2708  EPROM's) .  One 
EPROM  contains  the  standard  1  k  MONITOR  provided  with  the 
SDK-80  kit.  The  MONITOR  provides  software  to  communicate 
with  a  TTY,  enter  and  edit  user  programs,  examine,  list  or 
change  memory  locations,  and  initiate  execution  of  programs 
with  breakpoint  capability.  In  addition  to  providing 
general  system  control,  several  MONITOR  routines  are  called 
and  utilized  by  other  programs  developed  for  this  ap¬ 
plication.  Utility  programs  were  required  for  reading  and 
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writing  paper  tapes  on  the  TTY.  Unfortunately  Intel  did 
not  provide  these  programs  in  the  SDK-80  MONITOR.  Paper 
tape  input/output  programs  were  obtained  from  the  Insite 
[122]  user  library. 

In  order  to  carry  out  mathematical  operations  on  a 
microprocessor  beyond  the  simplest  binary  additions  and 
subtractions  a  set  of  so-called  floating  point  routines  is 
required.  We  initially  adapted  a  floating  point  package 
from  Findly  [123]  to  use  on  our  system.  However  it  soon 
became  both  too  time  and  memory  space  consuming  to  set  up 
a  complete  floating  point  package  on  this  small  a  micropro¬ 
cessor  system  beyond  routines  for  floating  point  normalizat¬ 
ion,  addition,  subtraction,  multiplication,  division  and 
TTY  output. 

Full  floating  point  capability  is  now  provided  on  our 
system  using  the  AM9511  Arithmetic  Processing  Unit  of 
Advanced  Micro  Devices.  The  chip  is  essentially  a  complete 
hardware  floating  point  package.  The  AM9511  contains  an 
internal  stack  onto  which  operands  are  loaded  from  the  8080 
CPU  over  the  system  data  bus.  The  appropriate  execution 
command  is  sent  to  the  AM9511  and  the  result  of  the 
operation  is  left  on  the  stack.  A  command  summary  for  the 
AM9511  is  shown  in  Table  X  and  schematic  for  the  hardware 
interface  to  the  SDK-80  is  shown  in  Figure  76.  It  is 
difficult  to  overemphasize  the  extent  to  which  the  use  of 
this  chip  expedited  overall  software  development  for  this 
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TABLE  X 


AM9511  COMAND  SUMMARY 


FIXED  POINT. ...  16  or  32  BIT 

ADD  MULTIPLY 

SUBTRACT  DIVIDE 


2.  FLOATING  POINT. . . .  25  BIT  MANTISSA, 


ADD 

SUBTRACT 

SQUARE  ROOT 

COS 

ASIN 

ATAN 

LN 

X  TO  POWER  Y 


MULTIPLY 

DIVIDE 

SIN 

TAN 

ACOS 

LOG 

LOG 

EXP 


3.  DATA  MANIPULATION 


FLOATING  POINT  TO  FIXED  POINT 
FIXED  POINT  TO  FLOATING  POINT 
CHANGE  SIGN 
READ/LOAD  STACK 


7  BIT  EXPONENT 


16  and 
32  BIT 
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76.  Interface  circuit  for  the  AM9511  to  the  SDK-80. 
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system. 

All  programs  were  placed  in  EPROM  (2708's).  The 
2708 's  were  programmed  using  the  PROMPT-80  system  of  Intel. 
Program  development  and  assembly  were  carried  out  using  a 
cross-assembler  running  on  a  PDP-11  minicomputer  in  our 
laboratory . 

Results  and  Discussion 

The  TTY  output  for  a  typical  run  (v-to-f  system)  for 

the  measurement  of  an  analyte  emission  signal  from  an  ICP 

source  is  illustrated  in  Figure  77.  The  operator  types 

in  the  desired  integration  time  (1.00  sec)  and  number  of 

replicate  measurements  (16  =  10,,).  Note  that  the  number 

1  b 

of  repeats  must  be  designated  in  HEX.  This  is  a  typical 
concession  that  one  often  has  to  make  in  dealing  with  a 
minimum  conf iguration  microprocessor  system  such  as  the 
single  board  SDK-80.  In  our  case  it  simply  was  not  worth 
the  programming  time  and  in  particular  the  programming 
space  to  include  a  conversion  routine.  The  integration 
time  could  be  input  in  decimal  format  because  the 
second  counter  on  the  timer  8253  was  being  operated  in  a 
BCD  mode.  Finally  the  operator  indicates  whether  this  is  a 
background  or  sample  (analyte)  run  and  types  a  sample  id¬ 
entification  code.  The  measurements  are  then  carried  out, 
in  this  case  16  consecutive  1  sec  integrations,  the  meas- 


I NT  TIM?  0100 
#  REPTS?  10 

back (B) /samp (s) ?  s 

ID?  1 

+0 . 2526000E+04 
+0. 2543000E+04 
+  0 . 2530000E+04 
+0. 2534000E+04 
+  0 .  25  36  0  00E+04 
+0 .2534000E+04 
+0 . 2539000E+04 
+0.2549000E+04 
+0 . 2534000E+04 
+0.2535000E+04 
+0 . 2539000E+C4 
+  0. 2538000 E +04 
+0 . 2545000E+04 
+  0 . 2534  000E+04 
+0 . 2537000E+04 
+  0 . 255.1000E+04 


SUM:  +0. 4060400 E+ 05 
AVE:  +0 . 2537750E+04 
BSB :  +0 . 24  61750E  +  04 
STD:  + 0 . 6 5 8 7 8 6  / E+ C 1 
S/N:  +0 . 373679 3E+0 3 


Typical  TTY  dialogue  and  output  for  a  run  with  the 
v-to-f  converter  based  Ficacurenent  system. 


Figure  77. 
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urements  are  summed,  averaged,  background  is  subtracted, 
standard  deviation  and  signal  to  noise  ratio  are  cal¬ 
culated,  and  all  data  are  logged  to  the  TTY.  The  background 
subtracted  is  that  measured  in  a  previously  acquired  and  so 
designated  run,  and  would  represent,  for  example,  the  output 
signal  while  running  a  blank  solution. 

The  TTY  dialogue  for  the  ADC  based  integrating  system 
is  essentially  identical  except  that  the  integration  time 
is  chosen  by  indicating  the  number  of  ADC  conversions  to  be 
summed.  A  1  sec  integration  is  equivalent  to  10,923  con¬ 
versions.  Based  on  signal  to  noise  ratios  (1/relative 
standard  deviation  when  n  =  32)  for  several  analyte  signals 
the  performance  of  the  v-to-f  and  ADC  based  systems  as  in¬ 
tegrating  measurement  systems  was  essentially  identical, 
i.e.  equivalent  signal  to  noise  ratios  were  obtained  for 
the  same  signal  and  integration  times  for  both  systems. 
However  practical  differences  do  exist  between  these  two 
systems  with  respect  to  factors  such  as  cost,  set-up  time 
and  expansion  to  multichannel  measurements. 

While  costs  are  admittedly  hard  to  pin  down  these 
days,  the  ADC  with  its  attendant  sample-and-hold  amplifier 
was  a  little  more  expensive  than  the  V-to-f  converter  and 
8253  counter/timer  chips  but  no  great  difference  occurs 
here.  In  terms  of  setup  time  (i.e.  interfacing  to  the 
SDK-80)  the  ADC  based  system  was  simpler,  as  it  was  carried 
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out  using  the  standard  on-board  chips  while  with  the  v-to-f 
system,  the  8253  based  programmable  frequency  meter  had  to 
be  developed  and  interfaced.  With  respect  to  expansion  to 
simultaneous  multichannel  measurements  (i.e.  for  a  direct 
reading  spectrometer  containing  several  photomultiplier 
tubes)  the  v-to-f  based  system  is,  perhaps,  more  effective. 

A  proposed  approach  is  illustrated  in  Figure  78.  Each 
channel  requires  an  operational  amplifier,  a  v-to-f  con¬ 
verter  and  a  12  bit  counter  (i.e.  2/3  of  an  8253).  A 
single  microprocessor  could  easily  control  the  complete 
system  as  once  the  counters  are  set  up  for  the  measurement 
they  need  not  be  read  until  the  integration  is  complete  and 
the  integrated  values,  i.e.  the  total  count,  can  be  held 
indefinitely,  it  being  a  digital  signal.  No  software  inter¬ 
vention  is  necessary  during  integration  and  in  general 
software  overhead  for  such  a  system  would  be  minimal.  Of 
course,  multichannel  ADC  systems  could  also  be  set  up  and 
in  fact  the  readout  systems  for  one  current  direct  reading 
spectrometers  [124]  is  based  on  analog  integration  at  each 
channel  with  subsequent  multiplexing  to  a  single  ADC.  In 
particular,  the  software  overhead  of  this  system  is  more 
complex  than  that  proposed  in  Figure  78,  in  part  because  of 
dynamic  range  limitations  associated  with  analog  inte¬ 
gration.  However,  our  goal  here  is  not  to  pick  a  "best" 
system,  but  simply  to  indicate  that  the  v-to-f  based  system 
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Figure  78.  Proposed  approach  for  a  multichannel  measurement 
system  based  on  the  v-to-f  converter-counter 
integration  system. 
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when  developed  with  software  programmable  counters  can 
potentially  provide  a  powerful  microprocessor  controlled 
multichannel  integrating  readout  system  for  PM  tubes. 

Conclusions 

Single  board  microprocessor  systems,  such  as  the  one 
used  in  this  application,  can  easily  provide  the  basis  for 
intelligent  data  acquisition  and  processing  systems.  The 
system  described  here  is  applicable  to  a  wide  variety  of 
spectrochemical  measurements  involving  photomultiplier 
tubes.  Versatility  is  built  in  as  a  result  of  the  cap¬ 
ability  of  programmable  chips  such  as  the  8253  counter/timer 
chip  and,  in  particular,  because  of  software  control  of  the 
system  via  EPROM  based  routines.  Indeed,  the  complete  single 
board  microprocessor  system  can  easily  be  dedicated  to  an 
entirely  new  task  simply  by  inserting  an  appropriate  EPROM. 

In  our  laboratory  this  system  has  been  used  for  ac¬ 
quisition  of  data  from  a  single  channel  inductively  coupled 
plasma  atomic  emission  spectrometer.  It  has  been  particularly 
useful  for  fundamental  signal  to  noise  ratio  measurements 
and  studies  where  it  has  relieved  the  tedium  of  collecting, 
calculating,  and  logging  the  large  amounts  of  statistical 
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APPENDIX  B 


This  program  is  used  to  collect  data  from  the  photo¬ 
diode  array.  The  program  runs  on  a  DECTAPE  based  PDP8/e 

minicomputer  with  an  OS/8  operating  system.  It  is  written 

in  FORTRAN  II  and  SABR.  Originally  written  by  E.G.  Codding, 

it  has  been  modified  by  E.D.  Salin  and  the  author. 
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611  FORMaTCFI  2.  4) 

DO  9  9  G  I  =  1 .»  J  ^ 

99  0  F  J  C I )  =  C  FL  J  AT  C  J  C I ) ) ) / 5  S 


GO  TJ  29 

9  63 

URI TEC Ij 502) 

502 

F 0  RM AT  C  ' BAG < GROU.JD  5 ’JBTRAC T I  0  Jj  Y  OR  ; J  ' ) 

READC  Ij  507) J3  CD 

,TRI  TEC  Ij  5  0  0  ) 

50  0 

FORi-laTC  'E.JTER  . 

JJ  OF  5Ca  J  S  A  J  D  P)I  JT3-.  214') 

A  o  i  i  D  C  1  j  j  u  1  )  J  j  .  ’> 

O  0  1 

FOPi'lATC  21  4) 

S 

r* 

o 

TAD  V) 

/EJTER  JO.  OF  SCAJ3 

s 

Cl  a 

/A  EG ATE  J 

r 

D  C  A  S  G  A  J 

/STJP.tL  IJ  SCAiJ 

s 

TaD  NOP 

/E.JTER  -JO.  OF  POINTS 

c 

CIA 

/.J  EG  ATE  JP 

S 

DC  A  POI.'J’r 

/STORE  IJ  PO  I M  T 

39 

CO  JTIO'TE 

s 

CLa  CEL 

r  —  n 

/REPEAT  CYCLE  STARTS  A"  39 

*  V  "■  J 

DO  1  I  =  1  j  .  J  -5 

1 

FJ  C  I  )  =  0 

') 


2 A  SI  4  ') 


& 


- 
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PEADC 1, 5 

C  7 )  J  G 

e  q  -  „  j  ~  . 

^  »  v>  *  -*-  •  — *  * 

T^J 

S3  A  J 

/Ghi"  JO 

.  OF  r-CAlc 

c 

DC  A 

1  2 

/ S FT  TA 

L  GCaJ  CO'IJTEP 

SG  J  , 

T  J 

•°OI.JT 

/J^~  TA 

2  J  J.  OF  aj  i  j  T  S 

q 

JCa 

i  •: 

/r'.3T  TA 

2  POI  IT  GJP.J^FP 

5 

TAD 

AG  177 

/GET  T;I 

L  c TAPT I  IC  .G)D'Arr 

s 

JG^ 

1  1 

/  rAT  TH 

£  storage  ADDRESS 

r- 

U  .  *  i  i 

r* 

o  I 

s 

oL  A 

S START  , 

SA 

0  :ip 

r  *r  \ 

S/  l  *  »  -  ‘  * 

/  Cil EGA 

FOR  STaTT  FLaG 

DOT  I 

/  P.EAD  I 

j  ~o\JT  3’»ftG~ 

DjjoI 

/CL AAA 

1. 1  PUT  L,tFF." 

s 

CF, 

_S 

J  M? 

f*  -t'  -x  »->  n* 

^  .  t-.  F\  1 

r* 

L*  A  4  —4 

r* 

_9 

SC  uC 

SCLOCA, 

SO  DF 

r» 

CLOG  A 

S 

SCGC 

PLCF 

r* 

CLL  °AP 

q 

lAn  I 

1  1 

/  3  TO  PE 

DtiTn  DJ  I  J  T 

r* 

J) 

I  S3 

1  0 

/ CALCA 

^0  I . I  T  GO  'I  J  "'E7-',  SA 1 

jap  clog  a 

C-L#  A  Ci-L 
I F C  A) 4 1 , 4  0 ,  4  1 


40 

DO  5 

0 

• — » 

•v 

■  ■  < 

II 

M 

5 

FJ  ( I  )  =  F 

loatc  j  c i )  >  +  :- j  ci ) 

r* 

G> 

JMP 

CAECA 

41 

D  J  4  2 

I  =  1 ,  J D 

42 

F  J  (  I  )  =  F 

J  C I ) -  FLOAT ( J  C I )  ) 

r*  j  -  ^  r 

•0  L/  •»  1  il«  Oi  'i  / 

I  S3 

12  /CAECA  SCaJ 

s 

J  A  n 

GO 

s 

v,EA  CLL 

IFC0  29 

30,  29 

30 

I  FCJ3-  1 

332) 29, 33, 29 

29 

P= FLOAT 

eg ) 

SI=  1 

DO  34  I  =  1,.JP 

34  FJ  (  I  )  =  FJ  (  I  )  /  P 

3  3  a.AaR=FJ(1) 

DO  17  1  =  2,  r3 

I  Ft  FJ  (  I  )  -it  i*vO  1  7,  1  7,  1  < 
.A-'a^  FJ  C  I  ) 

CJATI.J'JE 
AF1  J=FJ  C  1  ) 

DO  1  9 


I  F 


I  F 


13 
1  7 


I  =  2,  ;j  * 


. 


. 


I  FC  FJ  C I  )  -AMI  J )  20/  19^  19 
2  0  Ail  0  =  FJ  (19 

1 9  C  J JTI  JM£ 

?= 1 0  G  G . / C  AM  aX- AM I 0 ) 

I  FCX)  35/  36/  35 

35  ’’PI  TEC  1/  5  1  4)  S 

5  14  FORMATC  'SCALI03  FACT  J  P=  '/  F  1  2 .  4 ) 

DJ  231  1*1/. JP 

23  1  FJ  ( I  )  =FJ  C I  )  -AMID 

33  0  DO  23  2  I = 1 /  J P 

23  2  J  C I )  =FJ  C I  )  *  S 

CO  TO  1  5 

36  JO  23  1*1/ ,4? 

23  J  C  I  )  =  C  FJ  (  1  )  -  AM  I D  )  *  S 

16  M  R I  T  E  C  1/  503) 

5  0  3  FORMATC  '?  ') 

Hr^AJC  1  /  5  0  9  )00 
5  09  FJ  RMAT  (13) 

I  FC  JO)  16/  16/  37 

37  I  FC  3 -0  0)15/33/35 

33  COO  TI  DUE 

GO  TO  C  21/  22/  23/  24/  25/  26/  39/  7  0  0)  JO 
25  I  FC  J3-  1632)  3/  13/3 

13  READC 1/ 507) JG 

J  3=  0 
.<=  1 

S  JMP  SCAOS 

2  1  ’’FI  TEC  1/  701) 

701  FORMATC 'EJTER  FIRST  AO  D  LAST  DIODES  214.') 

READ  C  1/  7  02)  J  F/  JL 
70  2  FORMAT  C  21 4) 

MRI TEC  1/  5  0  5)  CJC  I  )/  I=OF/OL) 

505  FORMATC  1  0  C I  5/  IX)  ) 

MR  I  TEC  1/  7  0  3) 

7  0  3  FO RK  AT C  'FI  J I  SHED?  Y  OR  J •  ' ) 

^EAD  Cl/ 507)  JYO 
I  FC  JYO-  1  632)  21/  1  5/  21 
22  MR1TEC  1/5  10) 

510  FORMhTC  'ZD  TER  FILE  OaME/  A6') 

REaJC  1/  5  1  1  )  FILE 

511  FORMAT  CAS) 

Call  OOPE  J  C  'DTAl  '/ FILE) 

YRI  TEC4/  5C5)  C  J  C I  )  /  1=1/  J°) 

CALL  O CLOSE 
GO  TO  1 6 

2  6  mri  TEC  1/515) 

516  FORMaTC 'EJ TER  THE  THRESHOLD/ I  4  '  ) 

READ  C  1/ 50  1  )  JO 

MRI  TEC  1/517) 

517  FORMAT  C  'PO IDT  OO  .  AM PL I  TUDZ  SCALED  AM 
J  7  =  0  P-  3 


‘ 
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>  ,  :  i  3 


64 

65 

66 


67 


DO  63  1=6. 

I  f  C  J  ( I  )  -.JO)  6  3,  63,  64 
I  FC  J  C 1+  1 )  -J  C I  )  )  65,  6  5,  63 
I  F(  J  (  I  )  -J  ( I  -  1 )  )  63/  63,  66 
SJ  =  FJ  C I  ) 

I  3=1-7 
I  £= I  +  7 

sum=  o .  c 

DO  67  L=IS,  IE 

st;k=sum+fj  cl) 

vm  TEC  1  ,  5  1  6 )  I  ,  J  C  I  )  ,  SJ,  SUM 


513 

FO  PM  AT  C  1 X,  I  4,  9 

X,  I  4,  ax,  F6.  2,  4X,  FS.  0 

63 

COO  TI  HUE 

60  TO  16 

23 

VPI  TEC  1,  52  0) 
READ  C  1 ,  5  0  7 ) 0  S 

I  FC.JS-  1  532)  70, 

7  1,70 

7  1 

I  FCJb-  1  632)  75, 

76,  75 

76 

P.=  FLOAT  C.J ) 

DO  77  1=1,  OP 

77 

J  C  I  )  =  F  J  C  I  )  /  P 

60  TO  70 

75 

DO  7  2  I =  1 ,  J  P 

72 

J  C I  )  =  F  J  C  I  ) 

70 

O  C.JP+  1 )  =  50  0 

J  C.J  P+2)  =  10  23 

J  C  J  P+  3 )  =  5  0  0 

J3  =  OP+  3 

X=0 

S 

6 La  oLi_ 

s 

TAD  \.J  3 

5 

CIA 

S 

DCA  poi:j 

SSL  T, 

CL  A 

s 

TmD  POI J 

/GET  03.  OF  POIOTS 

S 

D  C  A  10 

/SE'r  POI  JT  COOOTLP. 

s 

tad  -<0177 

/GET  START I O  G  ADDS 

r* 

DCA  1  1 

/SET  STA-TIOG  ADD? 

S  PL  0  T  ^ 

CL  A  CM  A 

S 

Db  CO 

/CLEAR  O ’JTP’JT  LUFF 

c 

CL  A 

c 

TAD  I  1 1 

/GET  DATA  poiot 

s 

T  A  J  A  2  0  0  0 

/SET  OUTPUT  FLj-iG 

s 

DaSO 

/SET  0 UTP’.JT  JJFFEP 

r» 

D 

ASF 

/CHE CX  FOP  STOP  F 

j:;r!  cojt 

.iur 

\*/  —s  *  "i  WuL 


/C 


t.  r  -~'yS  \  ^  ’  r  ’ 

. .  .  i  —  .  k.'  u  A . .  ~  .  . 


■JO  D 


• 

' 
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I  F(  J  S- 

1  6  32)  7  3, 

74,  7  3 

n  /. 

t  ‘-X 

I FC  JB- 

1632) 330 

,  3  6,  3  3  0 

73 

CO  JTI  J 

UE 

Ji.P 

\  1  6 

/  2X1  T 

SCO J  T , 

I  n 

1  C 

/  CriEC/I  POIJT  CJUJTER, 

C 

-JAP 

PLOT 

/GET  J  EXT  POI.JT 

3 

JAP 

SET 

/START  OUEP 

24 

VPI  TEC 

1, 520) 

52  0 

FO.RMATC/'  SCAL 

E  PLOT?  y  OP  J  ') 

READ ( 1 

,  507)  J  S 

I  FC  4  S- 

1 632) 50, 

5  1,50 

5  1 

DO  5  2 

1  =  1,  JP 

52 

J  C  I  )  =  FJ  C  I  )  /2. 

50 

CU  JTI  J’JE 

S 

CL  A 

CLL 

r 

j 

TAD 

Y  J  p 

s 

CIA 

c 

DC  A 

PO I .  J  T 

r 

o 

CL  A 

CLL 

s 

TAD 

POI.JT 

/GET  .JO.  OF  POI  JTS 

s 

DCA 

1  0 

/SET  POIJT  CO  0 J  T£° 

s 

TaD 

K0  1  77 

/GET  STAPTI  J  3  ADD^E  c  S 

r 

DCA 

1  1 

/SET  STAPTI  JG  ADDRESS 

3  I 

DO  32 

1=  1,  33  0 

32 

J  G= 1*  1 

S 

CL  A  CM  A 

s 

DEC  J 

/CLEAR  OUTP'.JT  BUFFER 

s 

CLA 

s 

TAD  I 

1  1 

/GET  DATA  POI.JT 

s 

TAD 

K  2  0  0  0 

/SET  OUTPUT  FLAG 

s 

DO  SO 

/SET  OUTPUT  BUFFER 

c 

o 

CL  A 

s 

I  ST 

1  0 

/CHECK  POI  JT  CO  U  J  TER, 

s 

JMP 

\  3  1 

I  FC.J  S- 

1  532)  1 6,  74,  1  6 

SK  G  77  7, 

0777 

SK20Q0, 

20  0  0 

SPOI.J, 

0  0  0  0 

3P0I.JT, 

0  0  0  0 

SKO 177, 

0  1  77 

S  SC  AJ, 

0  0  0  0 

700 

DO  9  99 

1  =  1,  4  P 

F  J  C I  )  = 

FJ  C  I  )  *P 

999 

CO  JTI . J 

1JE 

CALL  CALI  (  FJ,  I  CHE,  5) 
CO J TIM US 


SK  I  P  I 


EYO 


SKIP  IF  ZERO 


GO  TO  1  6 

50  7  FJRMATCAl) 

A  CaLa  EXIT 

EJD 


APPENDIX  C 


This  program  is  used  to  read  256  point  lateral  data 
from  magnetic  tape,  perform  an  Abel  inversion  on  the  data, 
and  store  the  radial  data  back  on  tape.  Written  in  FORTRAN  II 
it  runs  on  a  DECTAPE  based  PDP8/e  under  an  OS/8  operating 
system.  The  program  is  written  in  two  parts;  a  main  pro¬ 
gram  and  an  Abel  inversion  subroutine.  It  was  written  by 
the  author. 
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CPROGRAtf  TO  PERFORM  ABEL  I D  VERS 1 00  00  DATA  FROM  256  ARRAY 
COMMO.J  FR/FL/JT/A 

DIME  J  SI  U  J  JC256)/FJC  25  6)/FRC256)/FLC256) 

/-i  —  1  ♦  C 

PROGRAM  TO  READ  FROM  DECT APE 


70  0  VRI  TEC  1  /  7  31) 

731  FORMAT  (  'EJTER  DEV., FILE  JAME/DO.  OF  PO I  I  TS/  C  2A6I  4)  ') 
REaDC  1/  7  32)  F1/F2/DP 
7  32  FORMaTC  2A6/  I  4) 

J  =  1 

CALL  I  0  °EJ  C  F  1  /  F2 ) 

READ (4/ 5C5)CJCI)/I=1/  JP) 

505  FORMAT C  1  0  (  I  5/  IX)  ) 

P0RTIOJ  TO  SCALE  THE  SPECTRUM  READ  FROM  TaPE 
URI  TEC  1/  733) 

733  FORMATC  'LITER  SCALIJG  FACTO RC F6.  3)  ') 

READC  1/734)  S 

7  34  FORMATC  F6.  3) 

DO  9  9  :J  =  1/  .1  P 

99  FJ  C  D  )  =  C  FLOAT  CJC.J)))/S 

DO  5  6  1  =  1/ DP 

5  6  FL  C I  )  =  F  J  C  I  ) 

C 

25  VRI TEC  1/ 60) 

60  FORMATC  'A7ERAGEC  1 )/  STOPEC  2)  /  ABEL C  3) /  EJDC  4) / RESTART C 5)  ') 
READC 1/ 61 )  J 

61  FORMATC  II) 

GO  TO C  60  0/  183/  1  04/  29  1/  70  0)1 
C  EJTER  CEJTRE  OF  SYMMETRY  OF  AXIAL  PLOT 

60  0  VRI TEC  1/ 601) 

60  1  FORMATC  'EJTER  START ID  G  DIODE') 

READC 1/ 602) IMI D 

60  2  FORMATC  I  3) 

VRI  TEC  1/  60  3) 

603  FORMATC  'EJTER  EJD  DIODE/  CO  OR  256)') 

READC  1/  60  2)  I  EJD 

JT=I A3SC IMI D-I EJD) 

C  AVERAGE  BOTH  SIDES  OF  AXIAL  TO  GIVE  123  EL EM ED T  AXIAL 

DO  100  1  =  1/  JT 
J  =  I M  I  D+  1  -  I 
K= IMI D+I -  1 

10  0  FLCI  )  =  CFJCJ)  +  FJCX)  )/2. 

C  ?ERO  U’JT  THE  EJD  EL  EXE  ITS 

D  =  1  23 -DT 
DO  10  1  I =  1  / :  J 
X=DT+I 

10  1  FL  C  X )  =  0 . 


. 

. 
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D 0  1  j  j  1=1^1  2o 

1  0  3  FJ  C I  )  =  F_  C I  ) 

GO  TO  05 


1  Do 

0  T=  1 20 

DO  ajLL  IJ7£r?IJ.J  JJ  LATERAL  ARP 
CaLL  A3LL 

DO  20  0  ,J=bl  23 

20  0 

FJ  (  1 )  =  FPC.J ) 

GO  TO  25 

SCALE  ARRAY  AJD 

STORE  OJ  MaG  TA° 

1  3  3 

JT=  1  23 

AH AX*  FJ  (  1  ) 

DO  17  1  =  2,  JT 

I  FC  FJ  C  I  ) -AH AX)  1  7 

,  17,  13 

18 

iJlAX=  FJ  C  I  ) 

1  7 

CJ  JTIO’JE 

AX  1 J  =  FJ  (  1  ) 

DJ  19  1  =  2,  JT 

I  FC  FJ  C  I  ) -AMI  J  )  20 

,  19,  19 

2  0 

AH  I  J  =  FJ  (  I  ) 

19 

CiUTIOOL 

S  =  1  0  U  0  •  /  C  i-i: 'J  A  ’  —  J 

V  “  I  T  E  (  1 ,  5  1  4 )  S 

I  J) 

5  1  4 

FOPLATC  *  SCiVL I  JO 

JO  23  1  I =  1 ,  : J  T 

FACTOR®  ',  FI  2.  4) 

£5  1 

FJ  C  I  )  *FJ  (  I  )  -  AM IO 
J J  2  0  2  I =  1 , O  T 

202 

J  C  I )  =  FJ  C  I  )  *  S 

R I  TEC  1,510) 

5  1  0 

FORMAT  (  'EO  TE”>  FI 
READC  1,511)  FIL 

L  E  J  AM  E,  A  5  '  ) 

5  1  1 

FORMAT  CD. 5) 

Call  0 o J C  '  DTA 1  ',  FILL) 
VRITEC4,  5  0  5) (J< I ), 1= 1,  JT) 
CaLL  0  CLOSE 
GO  TO  25 
29  1  STOP 

; :  j  ■> 


' 


. 
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S’JJFO’JT I  J  X  ABEL 

c  y-v.-u  j  ff*  jt*  a 

DI  M  L-J  SI  0.  J  FP <  25  S )  ->  FL  (  25  S> 

C=-2.  /(  3.  1  5  155*  A) 

JTT=  jT_  ] 

jj  id  x=i,.;t 

j  =  x 

X  J  =  F^OaT(  j  ) 

xx=fl:i  atcx> 

xx2=  cxx-  i.  )*  cxx-  i.  ) 

A  1  =  S  3P.T  (  C XX* *2)  -XX 2) 
a2=S  C  (XJ-  1.  )  *  *  2)  -XX  2) 

AX  J=  C  A  1  -  a2>  /  (  <  2.  *X  J  )  -  1 .  ) 

FP.C  C)  =-C*FL  (  J )  *AX  I 

IX  =  X  +  1 

JJ  1  1  1=30  ITT 

X  J=  -LOaTC  j  ) 

«1=  SOPTC  (X  J*X )  )  -XX2) 

i-.2=  5  3~''r  C  C  (  X  J  -  1 .  )  *  (  X  J  -  1  •  )  )  -  X  C  2 ) 

XX = FLOAT (J- 1 ) 

A  1  1  =  5  ?  FT  (  (  Xtt  * X!  1 )  -  X  C  2 ) 

A 2 2= SOFT C ( C  XX- 1 .  ) * C XX- 1 .  ) ) -”X£) 

A  ;  J  =  Ca1-A2)/  (  (2.  *XJ>-  1.  ) 

AX JM=(Al  1 -a22) / (  (2.*XF)-  1.  ) 

FF ( X )  =  FFCC)  +  (  C*  ~L  C  J  )  *  ( aX  J M-  AX  J )  ) 
1  1  CO  J  T I  J'TL 

F^CX)  =FF(X)  *  FLOAT  C.J  TT) 

VP. I  TL(  U  53)  X*  FL CC)  *  FRCX) 

6  6  FO  PM  aT  (  I  2r  1  2«  A) 

1C  GOJTIXUL 

PLT'JP  J 
L.J  J 


- 


